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PREFACE 

The  work  reported  herein  was  performed  by  Arnold  Engineering  Development  Center 
(AEDC),  Air  Force  Systems  Command  (AFSC),  under  Program  Element  65807F.  The  Air 
Force  Program  Manager  was  Dr.  K.  L.  Kushman.  This  work  was  performed  by  Calspan 
Corporation/ AEDC  Operations,  operating  contractor  for  the  Aerospace  Flight  Dynamics 
testing  effort  at  the  AEDC,  AFSC,  Arnold  Air  Force  Base,  Tennessee.  The  work  was 
performed  in  the  Aerospace  Systems  Facility  (ASF)  under  AEDC  Project  Number  DC71PW, 
The  task  in  the  DC71  project  addressed  here  was  .to  document  the  current  version  of  the 
CUNIFLOW  analysis  code.  This  report  was  submitted  for  publication  on  January  23,  1989. 
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1.0  INTRODUCTION 

Current  interest  in  hypersonic  vehicles  such  as  the  National  Aerospace  Plane  (NASP) 
and  the  Aeroassisted  Orbital  Transfer  Vehicle  (AOTV)  has  prompted  renewed  attention  to 
high-temperature  thermochemical  effects,  including  a  radiative  heat  load  to  the  body  and 
the  formation  of  a  plasma  sheath  which  will  affect  communication  with  the  vehicle.  At 
somewhat  lower  altitudes,  the  partially  dissociated,  vibrationally  excited  air  ingested  by  the 
NASP  engines  may  well  affect  the  efficiency  of  the  combustion  process  within  the  engines. 
Since  reproduction  of  the  hypersonic  operating  environment  within  a  ground-based  facility 
appears  to  be  impractical,  there  exists  a  powerful  impetus  for  numerical  simulation  of  the 
flight  conditions. 

The  Computational  Fluid  Dynamics  (CFD)  community  is  actively  working  to  develop 
the  capability  for  making  such  computations.  Gnoffo,  et  al.  (Ref.  1),  Li  (Ref.  2),  Prabhu, 
et  al.  (Ref.  3),  and  Candler,  et  al.  (Ref.  4)  have  recently  reported  Navier-Stokes  codes 
incorporating  nonequilibrium  thermochemistry,  but  these  codes  have  not  yet  reached  a 
production  stage,  and,  as  Barber  and  Cox  (Ref.  5)  point  out,  the  use  of  such  codes  for  design 
purposes  is  prohibitively  expensive. 

The  current  report  describes  a  mature  and  well-developed  code  which  provides  a  large 
part  of  the  capability  for  analyzing  the  aerophysical  problems  described  above  and  may  be 
as  close  to  production  status  as  any  code  of  similar  capabilities  will  ever  be.  This  code, 
CUNIFLOW,  is  now  operational  on  the  AEDC  Cray  X/MP  and  on  several  other  machines 
around  the  country  (Refs.  6  -  9).  The  original  version  of  CUNIFLOW  and  its  first  extension 
were  thoroughly  documented  in  Refs,  10  and  11,  but  these  documents  are  getting  hard  to 
find  and  there  have  been  many  changes  to  the  code  since  1967,  cf.,  Appendix  C. 

CUNIFLOW  calculates  the  inviscid,  hypersonic,  chemically  reacting  flow  field  about  a 
blunt  axisymmetric  or  two-dimensional  body  like  that  in  Fig.  1.  The  chemical  composition, 
pressure,  density,  and  velocity  components  of  the  free  stream  may  be  specified  as  functions 
of  position,  but  it  is  also  possible  to  input  only  the  vehicle  velocity  and  altitude  and  let  the 
code  calculate  the  other  free-stream  variables  from  fits  to  the  1953  U.  S.  Standard  Atmosphere. 
Embodied  in  the  code  are  thermochemical  models  for  treating  either  equilibrium  air  or  a 
chemically  relaxing  arbitrary  gas  mixture.  Additionally,  nonequilibrium  vibrational  relaxation 
for  two  diatomic  species  and  radiant  energy  loss  from  the  gas  may  be  modeled.  Since  the 
displacement  effect  of  the  boundary  layer  is  approximated  by  a  semi-empirical  method, 
CUNIFLOW  is  well  adapted  for  use  with  a  reacting  boundary-layer  code  such  as  BLIMP 
(Ref.  12). 
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The  advantage  of  the  interface  with  BLIMP  is  especially  apparent  at  low  altitudes  { =23,000 
ft)  where  more  sophisticated  programs  fail  to  converge.  OtheT  advantages  of  CUNIFLOW 
include  moderate  storage  requirements  and  high  speed.  As  currently  dimensioned, 
CUNIFLOW  requires  only  250,000  words  of  memory.  Run  times,  of  course,  will  depend 
on  the  body  size,  the  complexity  of  the  thermochemical  model,  and  the  specified  free-stream 
conditions.  However,  a  fairly  typical  shock  tunnel  simulation  (a  7-deg  blunt  cone  with  a  1-in. 
nose  radius  and  an  overall  length  of  10  in.)  was  calculated  in  69  sec  of  CPU  time  on  a  Cray 
X/MP. 

This  report  guides  the  user  through  the  current  version  of  the  code  and  its  usage  without 
repeating  all  the  details  contained  in  Refs.  10  and  11.  Sections  2.0  and  3.0,  together  with 
Appendix  A,  present  without  formal  derivation  the  mathematical  models  and  numerical 
methods  employed,  with  some  comment  on  their  validity.  Section  4.0  and  Appendix  B  describe 
the  input  procedures  and  suggest  some  cures  for  common  difficulties  which  may  arise.  The 
overall  code  structure  is  presented  by  the  tables  and  charts  of  Appendix  D. 

2.0  OPTIONS,  ASSUMPTIONS,  AND  RESTRICTIONS 
2,1  BODY  SHAPE 

A  typical  configuration  which  may  be  analyzed  with  CUNIFLOW  is  shown  in  Fig.  1 . 
The  body  must  be  convex  and  continuous;  that  is,  steps  are  not  allowed,  but  sharp  expansion 
corners  and  one  sharp  compression  corner  are  permissible. 

The  desired  body  may  be  specified  by  up  to  six  equations  of  the  form 

(r  -  rj)P  =  (z  -  Zj)Q  [C0  +  C,(z  -  Zj)  +  C2(z  -  zj)2...C4(z  -  Zj)<]  (1) 

where  z  and  Zj  are  measured  from  the  body  apex.  Since  the  stagnation  point  is  determined 
under  the  assumption  that  the  body  is  a  sphere,  it  is  necessary  that  n  and  Z|  be  zero. 
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As  indicated  by  Fig.  1,  the  inviscid  calculation  may  be  extended  beyond  the  end  of  the 
vehicle  by  the  artifice  of  introducing  a  “virtual  body”  and  replacing  the  tangency  condition 
by  a  specified  pressure  distribution.  The  surface  pressure  is  determined  as  part  of  the 
computation.  At  present,  only  a  constant  pressure  may  be  specified. 

The  physical  significance  of  this  constant  pressure  curve  is  minimal;  it  should  be  considered 
as  a  numerical  boundary  condition  which  is  required  by  the  numerics  of  the  reference-plane 
characteristics  method  (Ref.  13).  Although  the  computed  flow  field  in  the  base  flow  region 
is  of  dubious  validity,  the  flow  field  above  the  left-running  characteristic  from  the  corner 
will  be  valid. 

2.2  THERMOCHEMICAL  MODELS 

Five  separate  thermochemical  models  are  embedded  in  CUNIFLOW  and  may  be  selected 
through  the  input  constant  MODEL  (Section  4.2.4).  If  the  constant  is  not  specified,  the  default 
model  described  in  Section  2.2. 1  will  be  used. 

2.2.1  MODEL  =  0  (Default) 

The  gas  in  the  shock  layer  consists  of  Ns  monatomic  and  diatomic  species.  The  specific 
molar  concentration  of  the  jth  species  is  moles/gm.*  Because  the  finite-rate  relaxation  of 
the  vibrational  degrees  of  freedom  of  Nvi|>  species  is  to  be  modeled,  the  caloric  equation  is 
taken  in  the  form 


hj  =  Ij  +  h°  +  cj  ergs/mole  (2) 

where  ej  is  the  vibrational  energy  and  hj°  is  the  heat  of  formation.  The  energies  of  translation, 
of  rotation,  and  of  electronic  excitation  are  included  in  the  term  Ij,  which  may 
be  written  (Ref.  14) 


5  +  2(nj  -  1) 

RT  + 

Eif  Sjf  e  Ejf/T 

2 

m 

where  nj  denotes  the  number  of  atoms  per  molecule,  Ejj  is  the  energy  of  the  f1*1  electronic 
level  of  the  jlh  species,  and  gjf  is  the  degeneracy  of  this  level. 


*The  bar  (  )  signifies  a  dimensional  quantity. 
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The  chemical  potential  of  the  jth  species  is  given  by 


The  constant,  bj,  includes  the  translational  and  rotational  contributions,  as  well  as  the 
logarithm  of  the  ground-level  degeneracy.  The  characteristic  vibrational  temperature  is  By 


For  the  species  in  vibrational  equilibrium,  the  simple  harmonic  oscillator  approximation, 
(Ref.  14) 


(n,  -  1) 

[exp(0j/T)  -  1.0] 


(5) 


is  used  for  the  vibrational  energy.  For  the  nonequilibrium  diatomic  species,  the  vibrational 
energy  is  obtained  by  integration  of  a  differential  equation  relating  tj  to  the  temperature, 
T,  and  to  the  rates  of  certain  reactions. 


Each  of  the  component  gases  may  be  assumed  to  be  thermally  perfect,  so  that  one  may 
write  for  the  equation  of  state 


Ns 

p  =  q  T  y?  £  7j  (6) 

j  =  i 


where  y?  is  the  universal  gas  constant,  8.313405  x  107  ergs/mole  °K. 

The  default  thermodynamic  model  (MODEL  =  0)  is  based  on  certain  common 
simplifications  of  molecular  physics.  These  assumptions  are  useful  in  describing  the  thermal 
behavior  of  a  low-temperature,  low-density  gas,  but  become  increasingly  inapplicable  as  the 
temperature  increases.  In  particular,  the  simple-harmonic  oscillator  assumption  of  Eq.  (5) 
must  be  supplemented  by  anhannonicity  corrections  (Ref.  15)  for  temperatures  above 
approximately  1000°  K.  Furthermore,  the  separability  of  the  vibrational,  rotational,  and 
electronic  modes  implied  by  Eqs.  (2)  and  (3),  although  commonly  assumed,  is  not  in  fact 
true  (Refs.  15  and  16)  and  may  become  significant  at  temperatures  greater  than  2500°  K 
(Ref.  16).  Fortunately,  at  the  very  high  temperatures  (T  **  50,000°K)  predicted  for  the  AOTV, 
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MODEL  =  0  is  adequate  because  dissociation  is  so  rapid  that  the  vibrational  and  rotational 
modes  contribute  very  little  to  the  mixture  energy. 

2.2.2  MODEL  =  1 

The  gas  is  assumed  to  be  perfect  with  a  specific  heat  ratio  of  1 .4.  If  the  user  does  not 
specify  the  free-stream  composition,  air  will  be  assumed.  There  is  very  little  reason  for  selecting 
this  option,  as  the  code  numerics  are  not  well  adapted  to  it. 

2.2.3  MODEL  =  2 

Equilibrium  air  (02  and  N2  only)  is  assumed.  The  well-known  RGAS  and  SERCH 
routines  (Ref.  17)  are  used  to  obtain  the  thermal  properties.  No  composition  data  are  output 
and  the  radiation  option  is  turned  off. 

2.2.4  MODEL  =  3 

The  fluid  is  assumed  to  be  an  arbitrary  mixture  of  thermally  perfect  gases  —  monatomic, 
diatomic,  polyatomic  —  with  fully  equilibrated  internal  states.  The  specific  heat,  enthalpy, 
entropy,  and  chemical  potential  of  each  species  are  calculated  from  polynomial  fits  to  existing 
property  tables  (Refs.  18  and  19). 


t 


5* 


=  ajj  +  a2j  T  +  a3j  T^  +  T3  +  a3j  T4 


hj  a2j  -  a3j  —  —  aJj  —  a$i 

—  =  a,i  +  — T  +  — T2  +  —  T3  +  —  t4+- 

J  2  3  4  5  T 


(7) 

(8) 


—  aij  f n  T  +  a2j  T  + 


MT2 
2 


T2  x 


T4  + 


4 


+  a7j 


(9) 


(10) 


Gordon  and  McBride  (Ref.  20)  have  fitted  the  data  of  Ref.  18  with  polynomials  of  the 
form  of  Eqs.  (7)  -  (9),  and  furnish  an  extensive  tabulation  of  the  fit  coefficients  with  the 
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NASA  LeRC  Chemical  Equilibrium  Code  (Ref.  20)  in  the  format  required  by  CUNIFLOW. 
The  LeRC  Code  fits  are  valid  only  in  the  ranges  300°K  to  1000°K  and  1000°K  to  5000°K. 

MODEL  =  3  requires  that  the  internal  modes  —  rotation,  vibration,  electronic  —  be 
in  equilibrium  with  the  translational  temperature.  This  assumption  is  justified  for  most 
propulsion  problems  to  which  CUNIFLOW  might  be  applied  and  in  many  lower-speed,  higher- 
density  flow  problems,  except  directly  behind  a  thin  shock. 

Both  the  JANNAF  (Ref.  1 8)  and  the  Browne  (Ref.  19)  property  tables  incorporate  a  variety 
of  empirical  corrections  for  anharmonicity  and  vibration-rotation  effects  so  that  they  are 
more  accurate  than  the  MODEL  =  0  formulae  quoted  in  Section  2.2.1.  The  accuracy  of 
the  CUNIFLOW  results  will  depend  on  the  precision  of  the  curve  fits,  which  usually  deteriorate 
rapidly  outside  the  announced  limits.  CUNIFLOW  will  not  issue  any  messages  when 
extrapolation  occurs.  The  effect  of  such  extrapolations  on  the  solution  will  depend  on  the 
species  concentration,  on  the  reactions  in  which  the  species  is  involved,  and  on  the  magnitude 
of  the  extrapolation  error. 

2.2.5  MODEL  =  4 

The  thermodynamic  model  presumed  for  MODEL  =  4  is  very  similar  to  MODEL 
=  3.  The  thermal  properties  are  calculated  from  the  equations: 

Cpj  =  ®3j  +  a4j  T  +  ajj  T2  (11) 

hj  =  a,j  +  a2j  +  a3j  (T  -  3000)  +  —  (T2  -  9.0  x  106)  -  a5j  A  -  (12) 

2  \T  3000/ 


Sj  =  a* j  +  a3j  to/— ^  +  a4j  (T  -  3000)  -  ^  II-  -  —\  (13) 

\3000/  2  \T2  3000 2) 

which  are  those  used  by  the  KBLIMP-A  (Ref.  12)  boundary-layer  code. 

If  this  option  is  selected,  the  pressure,  entropy,  and  composition  at  the  body  surface  or 
boundary-layer  edge  are  output  on  Fortran  logical  Unit  3  in  a  format  acceptable  to  BLIMP.* 


'File  usage  is  described  in  Table  1,  Section  4.1. 
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2.3  CHEMICAL  AND  VIBRATIONAL  RELAXATION 


The  rate  of  change  of  the  concentration  of  the  j01  species  along  a  streamline  is  described 
by  two  equations: 


NS 

E 


i=i 


0 


k  =  1,  ...Nq 


(14) 


Nr  _  _  _ 

£7J  _  ^  Qij  (T»{l'}) 

as  ~  i=i  eU 


j  =  Nc  +  1 . Ns 


(15) 


in  which  Nc  represents  the  number  of  elements  and  ajk  the  number  of  atoms  of  kth  element 
in  species  j.  Equation  (14)  states  the  principle  of  atom  conservation.  Conservation  of  electric 
charge  is  secured  by  treating  the  electron  as  an  atom  and  assigning  ajk  =  - 1  if  the  jth  species 
is  a  positive  ion.^ 


The  source  term  on  the  right  of  Eq.  (15)  is  a  function  of  the  temperature  and  of  all  of 
the  species  concentrations  {7} .  When  vibrational  nonequilibrium  is  considered,  the  chemical 
source  term  Qy  also  becomes  dependent  on  the  rate  of  change  of  the  vibrational  energy  ej 
of  certain  species. 

For  the  Nyib  species  not  in  vibrational  equilibrium,  the  vibrational  energy  is  determined 
from 

flei  _ 

~  =  Wj  j  =  Nc+ 1,  ...,  Nc  +  Nvjb  (16) 


where  wj  is  coupled  to  the  chemical  source  terms  Qy  in  Eq.  (15). 

2.3.1  Reactions 

To  lend  generality  to  the  program,  the  Nr  reactions  are  written  in  the  form 
NS  NS 

X)  Vjj  Mj  V*  E  /ij  Mj  i  =  1,2 . Nr  (17) 

j=i  j=i 


tNote  that  the  convention  for  numbering  the  species  lists  the  monatomic  elements,  including  the 
electron,  First,  followed  by  the  molecular  species  with  nonequilibrium  vibrational  energies. 


11 


AEDC-TR-88  31 


where  Mj  represents  the  j»h  chemical  species,  and  the  stoichiometric  coefficients  are 
represented  by  i»y  and 

The  forward  and  backward  rate  constants  for  the  ith  reaction  are  related  by 


where 


and 


(18) 


(19a) 


(19b) 


By  using  Eq.  (18),  either  the  forward  or  backward  reaction  rate  may  be  specified  and 
the  other  computed  from  it.  The  temperature  dependence  of  the  reaction  rate  constants  is 
assumed  in  the  form 


kfieq  0r  kbi«i  =  Aki  TBki  eXp 


The  source  term  Qy  is  written* 


Qij  -  Wj 


(20) 


(21) 


tThese  formulae  assume  that  Z  and  D  reactions  are  binary  in  the  forward  direction. 
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In  Eq.  (21),  the  parameters  Wif  Dj,  Z;  may  be  either  zero  or  one  and  serve  to  denote  the 
type  of  reaction.  Thus,  the  first  term  alone  (Wj  =  1,  Zj  =  0,  Dj  =  0)  accounts  for  specific 
collision  partners  as  in  the  reaction  02  +  O  **  20  +  O.  The  second  term  is  used  when  all 
collision  partners  are  lumped  together,  as  in  the  reaction  02  +  M  20  +  M.  The  third 
term  is  used  when  some  collision  partners  are  lumped  and  others  specified.  Note  that  for 
a  Z-type  reaction,  the  specified  reactant  may  not  be  included  among  lumped  collision  partners. 
For  instance,  N2  +  M  **  2N  +  M  with  M  including  N2  will  not  permit  N2  to  collide  with 
itself  if  the  reaction  is  specified  as  a  Z  type.  This  must  be  written  as  a  separate  reaction  of 
type,  i.e.. 


N2  +  N  s*  2N  +  N  W  Type 

■N2  +  M  ^2N  +  M  Z  Type 


The  parameter  Xj  appearing  in  Eq.  (21)  is  defined  as 


X; 


=  1  - 


(22) 


where  Kj  is  the  equilibrium  constant.  The  significance  of  X;  is  that  it  represents  the  degree 
of  nonequilibrium 

reverse  rate  (23) 

X  =  1  —  - 

1  forward  rate 

and  is  obviously  equal  to  zero  at  final  equilibrium.  Immediately  behind  the  shock  Xj  =  1**, 
while  in  a  region  of  rapid  recombination,  xt  may  assume  large  negative  values. 

2.3.2  Vibrational  Relaxation 

Since  the  energy  required  for  dissociation  is  diminished  by  the  internal  energy  stored  in 
vibration,  a  finite  rate  of  vibrational  relaxation  affects  the  dissociation  rates.  On  the  other 
hand,  the  upper  vibrational  levels  are  selectively  depopulated  by  dissociation.  In  other  words. 


‘These  formulae  assume  that  Z  and  D  reactions  are  binary  in  the  forward  direction. 
“If  all  reactions  are  input  with  dissociation  on  left  side  of  Eq.  (17). 
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vibrational  and  chemical  relaxation  are  coupled.  The  effect  of  the  nonequilibrium  vibrational 
energy  on  the  reaction  rates  is  described  by 


"Svib 

kfl  =  kf  II  VAti  (24) 

^  i=i  i 

and  the  subscript  “eq"  on  the  reaction  rates  indicates  that  these  are  evaluated  for  equilibrium 
vibration.  The  exponent  Ay  may  be  either  0  or  1  and  is  used  to  specify  which  reaction  will 
be  affected  by  the  vibration-dissociation  coupling  factor, Vj  [Eq.  (28)]. 

In  developing  an  estimate  of  the  so-called  CVDV  coupling  terms  which  describe  these 
opposing  effects,  it  is  assumed  that  the  vibrational  energy  levels  are  those  of  a  truncated 
simple  harmonic  oscillator,  that  the  levels  are  populated  according  to  a  Boltzmann  distribution 
with  a  vibrational  temperature  Tvib,  and  that  the  probability  of  dissociation  depends  only 
on  the  collision  energetics.  For  this  model  (Refs.  21,  22),  the  relaxation  rate  uj  is 


Wi  =  ~~z  — 


d«j  eieq  "  ei  1 


dS 


7]U  eU 


j  \  Njg_i  _  I  g  Ajj 

|exp(0j/Tmj  -  1)  cxp(0jN'j/Tmj  -  1)  J  j=i 


Nr 


Gu 


(25) 


where  N’j  is  the  number  of  vibrational  levels  in  the  truncated  oscillator  and 


T„,  «)  l  <j  /  T 

The  characteristic  relaxation  time  r,  is  assumed  to  be  of  the  form 


(26) 


raj  T^i 
P 


exp 


(27) 


where  Taj,  rc.,  are  user-supplied  constants.  For  the  nonpreferential  model  (Ref.  21), 
the  coupling  factor  is 
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1 

'l  -  exp(  -Njflj/Tjnj)' 

'exp(0j/Tvj)  -  1 

N7] 

exp  (Bj/Tm.)  -  1 

exp  (flj/T)  —  1 

(28) 


where  the  vibrational  temperature  Tvj  is  just  the  reciprocal  of  the  first  term  on  the  right- 
hand  side  of  Eq.  (26),  and  T  is  the  translational  temperature. 

At  very  high  velocities  (greater  than  30,000  ft/sec),  the  vibrational  coupling  model  outlined 
above  breaks  down,  resulting  in  the  unstable  growth  of  vibrational  energy  and  temperature. 
Since  the  rate  of  dissociation  is  high,  only  a  small  fraction  of  the  total  enthalpy  is  involved 
in  the  vibrational  modes,  so  that  the  use  of  equilibrium  vibration  is  a  reasonable  approximation. 

Note  also  that  the  coupling  model  is  intended  only  as  a  correction  to  empirical  reaction 
rates  measured  under  conditions  that  assure  equilibrium  vibration.  The  CVDV  coupling  model 
has  recently  been  discussed  and  modified  by  Park,  et  al.,  (Refs.  23,  24)  and  by  Carlson,  et. 
al  (Ref.  25)  and  is  currently  in  a  state  of  flux. 

2.4  RADIATION 

If  the  nonequilibrium  option  (MODEL  =  0)  is  selected,  it  is  also  possible  to  consider 
a  class  of  reactions  which  result  in  the  emission  of  radiation.  Since  the  shock  layer  is  assumed 
to  be  optically  thin,  this  radiant  energy  is  lost  to  the  gas. 

The  radiation  from  the  excited  particles  is  computed  as  if  it  were  in  equilibrium  with  the 
local  translational  temperature;  i.e.,  the  energy  states  which  are  involved  in  radiative  transitions 
are  assumed  to  be  populated  according  to  a  Boltzmann  distribution  with  the  translational 
temperature  as  an  independent  variable.  The  overall  rates  of  radiation  from  the  various 
radiating  systems  were  taken  from  Ref.  26.  These  rates  can  be  fitted  over  a  considerable 
range  by  expressions  of  the  type 

i' 

AkNj  exp  (-  Bk/f)  (29) 

where  Ak  and  Bk  are  input  constants  and  Nj  is  the  nonequilibrium  number  density  of  the 
j*  species. 

2.5  FREE-STREAM  OPTIONS 

To  study  entry  into  a  nonterrestrial  atmosphere  or  phenomena  in,  say,  an  arc-heated  tunnel, 
the  chemical  composition,  pressure,  density,  and  velocity  components  of  the  free  stream  may 
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be  specified  as  functions  of  position.  Most  of  the  time,  however,  it  is  appropriate  to  consider 
a  uniform  stream  of  standard  air  and  input  only  a  velocity  (ft/sec)  and  an  altitude  (feet) 
and  let  the  code  calculate  the  other  free-stream  variables  from  fits  to  the  1953  U.  S.  Standard 
Atmosphere  (subroutine  HATM). 

2.6  BOUNDARY-LAYER  OPTIONS 


One  of  the  strong  features  of  CUNIFLOW  is  its  ability  to  approximate  some  of  the  effects 
of  a  boundary  layer  without  ever  actually  calculating  a  viscous  profile.  This  procedure  needs 
a  wall  enthalpy  ratio  (or  wall  temperature,  Tw)  which  must  be  input. 

While  the  boundary-layer  procedure  increases  the  realism  of  the  simulation  and  allows 
the  program  to  avoid  the  difficult  “reaction  layer”  near  the  wall,  it  slows  down  the  calculation 
and  sometimes  introduces  numerical  instability.  Because  the  numerical  methods  used  at  slope 
discontinuities  are  incompatible  with  the  presence  of  a  boundary  layer,  the  boundary-layer 
calculation  should  be  omitted  when  discontinuities  are  present. 

2.6.1  Transport  Properties 


The  viscosity  must  be  evaluated  in  connection  with  the  boundary-layer  matching  procedures 
to  be  described  later.  Since  the  boundary  layer  is  calculated  only  approximately,  an  exact 
calculation  of  the  viscosity  is  not  required.  At  the  stagnation  point,  a  program  based  on  the 
work  of  Peng  and  Pindroh  (Ref.  27)  is  used  to  compute  the  stagnation  density  and  stagnation 
viscosity.  In  the  program,  an  equilibrium  air  composition  involving  14  species  is  obtained 
as  a  function  of  temperature  and  pressure.  After  computing  the  density,  the  program  lumps 
the  original  14  species  into  4  (molecules,  atoms,  ions,  electrons)  with  molecular  weights  m] 
=  28.8  gm/mole,  m2  =  14.4  gm/mole,  m3  =  14.4  gm/mole,  and  1114  =  5.5  x  10-4 
gm/mole,  respectively.  Using  the  same  subscripts  for  the  lumped  mole  fractions,  Xj,  one 
may  write  an  approximate  expression  for  the  mixture  viscosity  as  follows: 


e 


1.45 


Xk 


,-1 


i=l  I  Mi  G  XiOlj  k  =  l,k*i  Dik 


(30) 


where  the  viscosity  of  the  pure  species  is 


w  = 


ViPi  T 

11920*22 


(31) 
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1  / - \W2 

—  =  171.1[=?3_I  T-^yP  (32) 

D'j  \mi  +  mj/ 

The  cross  sections  for  viscosity  <j>a  and  binary  diffusion  <j>\}  were  obtained  by  a  least-squares 
curve  fit  of  the  results  of  Peng  and  Pindroh. 

The  stagnation  point  viscosity  evaluated  from  Eqs.  (30)  to  (32)  is  used  as  a  reference  for 
viscosity  calculations  at  all  other  points  of  the  nose-region  boundary  layer  since  it  may  be 
rather  safely  assumed  that  the  flow  close  to  the  body  is  chemically  frozen  at  a  composition 
near  that  of  the  stagnation  point  equilibrium.  Away  from  the  stagnation  point  in  the  nose- 
region,  the  viscosity  is  computed  as 


m  =  poVttx;  03) 

Over  the  afterbody,  temperatures  get  low  enough  (~  3000°  K)  that  the  Sutherland  formula 


(1.4582  x  10- 5  )T1/2 

a  =  - — - gm/cm-sec 

T  +  110.3 

is  applicable.  The  reader  should  note  that  Eq.  (34)  applies  only  to  air. 


2.6.2  Displacement 


(34) 


If  INBLR  =  I  (Section  4.2.4),  the  displacement  interaction  between  the  viscous  and  inviscid 
flow  fields  is  computed  using  a  simplified  boundary-layer  theory.  While  this  boundary  layer 
is  approximate,  the  neglected  items  —  chemical  reactions  and  species  diffusion  —  have  no 
significant  effect  on  the  variables  which  affect  the  inviscid  flow  field,  namely,  the  mass  flux, 
ife.  through  the  viscous  layer  and  the  velocity  thickness,  8.*  The  basic  relations  governing 
these  variables  are 


Ik  =  (2ir)“  \/2  s  F,  (n,  Gwall) 


(35) 


V2I 

8  -  =^r-^ 
PeUerg 


(36) 


'The  displacement  thickness  is  not  used  because  it  tends  to  be  negative  under  hypersonic  flight 
conditions. 
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where 


=i: 


Qt  fit  u*  rg"  dS 


(37) 


is  the  transformed  arc  length  used  in  laminar  compressible  boundary-layer  theory.  Away 
from  the  nose,  it  is  necessary  to  replace  the  customary  pressure  gradient  parameter,  v,  by 
its  average  over  the  length  of  the  boundary  layer. 


For  consistency  this  usage  was  extended  to  the  nose  region  boundary  layer.  The  quantities 
Fj,  Fj,  and  Fj  appearing  in  Eqs.  (35)  and  (36)  are  integrals  over  the  ideal  gas  boundary-layer 
profiles  reported  by  Cohen  and  Reshotko  (Ref.  28).  The  values  of  these  integrals,  taken  from 
Ref.  28,  were  fitted  by  the  third-degree  polynomials: 


Fi  =  1  f  dij 

0 

=  ajo(Gwa]i)  +  aj|(Gwan)  it  +  a,2(Gwaii)xJ  +  a13T3 

(39) 

F2=  1 

• 

rcu 

L  ^  =  a20  +  •  • 

•  +  ^(Gwqi)*3 

(40) 

F3  =  ' 

P  di>  =  a3o(Gwau)  +  .  . 

.  +  a33(GwaU)7r3 

(41) 

where 


G  =  H(ij)/He;  f  =  -^-  =  U(,);  ,  = 


dy 


For  a  given  value  of  Gwan,  the  coefficients  ay  are  obtained  by  table  look-up  and  linear 
interpolation.  The  evaluation  of  the  Fj  is  performed  by  subroutine  INTGRL. 


3.0  NUMERICAL  METHODS 


A  complete  presentation  of  all  the  numerical  details  would  make  this  report  unwieldy 
and  contribute  little  to  the  would-be  user's  understanding  of  what  is  going  on  in  the  calculation. 
Since  most  of  this  information  is  available  in  Refs.  10  and  11,  this  section  will  attempt  to 
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present  only  the  information  necessary  for  understanding  the  input  procedures  and  for 
evaluating  the  output  of  the  code.  Additional  information  will  be  found  in  Appendixes  A 
and  D. 

Since  the  shock  layer  flow  is  both  subsonic  and  supersonic,  distinct  numerical  methods 
are  required  for  the  flow  field  over  the  nose,  which  is  described  by  elliptic  partial  differential 
equations,  and  for  the  remainder  of  the  body,  which  is  described  by  hyperbolic  flow  equations 
which  permit  the  use  of  the  method  of  characteristics.  The  evaluation  of  the  thermochemical 
properties  and  of  the  reaction  rates  is  common  to  both  procedures,  as  is  the  evaluation  of 
the  boundary  conditions  at  the  thin  shock. 

This  distinction  in  numerical  methods  is  reflected  in  the  structure  of  CUNIFLOW  (cf., 
Appendix  D)  which  consists  of  a  small  MAIN  program  which  calls  BREGN  to  compute  the 
blunt  body  flow,  FRMAIN  to  calculate  the  supersonic  flow,  and  BRIDGE  as  a  transition 
between  the  two.  The  program  is  arranged  so  that  the  named  segments  can  be  loaded  as 
overlays,  although  this  is  not  required  for  most  modern  computers  such  as  the  Cray  X/MP. 

3.1  THE  NOSE  REGION  PROGRAM 

3.1.1  General 

The  elliptic  partial  differential  equations  which  describe  the  flow  in  the  subsonic-transonic 
region  between  the  body  and  the  detached  shock  wave  are  solved  by  the  “Method  of  Lines” 
(Ref.  29).  That  is,  derivatives  parallel  to  the  shock  are  evaluated  numerically  so  that  the  partial 
differential  equations  may  be  converted  to  a  set  of  ordinary  differential  equations  which  are 
then  integrated  by  an  explicit  marching  procedure  starting  from  the  shock  and  marching  inward 
toward  the  body. 

The  body  surface  is  located  by  mass  conservation.  An  iterative  procedure  is  used  to  vary 
the  shock  shape  until  the  resulting  body  shape  corresponds  to  that  specified.  When  the 
boundary-layer  displacement  option  is  elected,  the  mass  balancing  takes  place  at  the  edge 
of  the  boundary  layer. 

3.1.2  The  Nose-Region  Mesh 

The  geometry  for  the  nose-region  calculation  is  shown  in  Fig.  2.  The  mesh  consists  of 
data  surfaces  parallel  to  the  shock  and  rays  normal  to  the  shock.  Since  the  data  surfaces 
are  parallel  to  the  shock,  there  is  a  discrepancy  between  grid  and  body  which  increases  as 
one  moves  away  from  the  axis.  The  calculation  on  the  six  outer  rays  is  not  carried  completely 
to  the  body  surface  (Section  3.1.4).  A  minimum  of  seven  equally  spaced  rays  is  required; 


19 


AEDC-TR-88-31 


25  rays  are  allowed.  The  ray  spacing,  Ax,  and  the  position  of  the  first  off-axis  ray,  xo,  may 
be  adjusted  to  ensure  a  satisfactory  supersonic  ray  for  starting  the  M.O.C.  (hyperbolic  region) 
program.  If  x<>  <  Ax,  BREGN  will  make  use  of  symmetry  to  create  three  virtual  rays  so 
that  x-derivatives  on  the  axis  will  be  handled  properly  by  the  standard  derivative  routines. 
The  axis  will  then  appear  in  the  printout  as  ray  number  4. 


(Note  the  "virtual"  grid  points  below  the  symetry  axis.) 

Figure  2.  Physical  and  computational  geometries  for  the  blunt  body  code. 

The  spacing  in  the  y-direction  (normal  to  the  shock)  is  determined  dynamically  (within' 
user-specified  limits)  by  the  program.  Spacing  will  be  fine  near  the  shock  and  near  the  body. 


To  relate  the  physical  coordinates  z  and  r  to  the  computational  coordinates  x  and  y,  it 
is  necessary  to  specify  a  shock  shape  in  some  parametric  form.  The  shape  recommended  by 
Lomax  and  Inouye  (Ref.  17) 


z 

RSH 


(42) 


is  used  for  this  purpose.  This  shock  has  a  radius  of  curvature  (RSH)  at  the  axis. 
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The  metric  for  the  nose-region  mesh  is 


dS2  =  (1  +  xy)2  dx2  +  dy2  +  rM^  (43) 

where  x(x)  is  the  local  shock  curvature,  da/dx  (Fig.  2). 

3.1.3  Boundary  Conditions 

When  the  free  stream  is  nonuniform,  the  local  pre-shock  conditions  are  calculated  by 
linear  interpolation  in  the  two-dimensional  tables  input  via  Unit  4  (cf..  Table  1). 

Since  the  shock  is  assumed  to  be  thin,  the  composition  of  the  gas  does  not  change  during 
its  passage  through  the  shock  unless  equilibrium  has  been  assumed  (MODEL  =  2).  The 
vibrational  energies  of  NVIB  (Section  4.2.4)  molecules  are  held  at  their  free-stream  values 
during  the  shock  passage.  The  vibrational  energies  of  the  remaining  molecules  are  equilibrated 
to  the  post-shock  translational  temperature. 

Because  of  the  variety  of  possible  thermochemical  options  involved,  the  shock  conservation 
conditions  are  solved  by  an  iterative  process. 

3.1.4  Marching 

Sections  3.1.4  and  3.1.S  describe  the  inviscid  procedure.  The  optional  viscous  corrections 
are  discussed  in  Section  3.1.6.  In  the  following  discussion,  “step"  means  advancing  by  the 
same  increment  Ay  on  each  ray. 

A  marching  cycle  begins  by  using  the  most  recent  data  to  evaluate  the  derivatives  in  the 
x-direction.  A  seven-point  algorithm  (Appendix  A)  derived  from  filter  theory  (Ref.  30)  is  used. 

Since  the  x-derivatives  are  now  known,  one  has  on  each  ray  a  coupled  set  of  ordinary 
differential  equations  which  are  integrated  for  a  single  y-step  by  an  algorithm  attributed  to 
Lomax  (Ref.  31).  The  algorithm,  which  consists  of  an  explicit  predictor  and  4,  6,  8,  or  10 
implicit  correctors,  is  presented  in  Appendix  A. 

The  Lomax  algorithm  was  devised  to  overcome  the  step  size  limitations  associated  with 
the  “stiffness"  of  the  differential  equations.  A  set  of  differential  equations  is  said  to  be  stiff 
when  the  set  possesses  large  negative  “parasitic"  eigenvalues  which  limit  the  step  size  but 
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are  uninteresting  physically,  and  some  very  much  smaller  “driving”  eigenvalues  which  convey 
the  interesting  physics. 

Failure  of  an  integration  step  is  indicated  by  a  negative  value  for  a  dependent  variable 
or  by  an  excessive  change  (15  percent)  in  the  temperature.  When  a  failure  is  detected,  the 
previous  (successful)  step  size  is  halved  and  the  entire  step  is  repeated.  The  possible  repetition 
of  a  successful  step  requires  that  storage  be  provided  for  two  complete  data  surfaces  (BO 
and  B1  in  Fig.  2)  as  well  as  the  surface  under  construction,  B2. 

If  a  step  is  successful,  the  data  in  array  BO  are  written  to  Unit  9,  B1  is  copied  to  BO, 
B2  is  copied  to  Bl,  and  the  size  of  the  next  step  is  increased  according  to  the  rule 

STEP  =  STEP  X  (1.0  +  1.0/FLOP)  (44) 

where  FLOP  is  initialized  to  1.0  and  incremented  by  1.0  at  each  failure.  In  addition,  the 
estimated  distance  to  the  body,  calculated  at  each  step,  is  used  as  an  upper  bound  on  the 
step  size. 

i 

The  body  is  located  by  mass  balancing.  At  a  general  point  on  ray  M,  the  mass  flux  across 
the  ray  between  the  point  and  the  body  is 


where  V'sh  is  the  nondimensional  mass  flux  into  the  shock  between  the  axis  and  the  point 
(x,0).  The  inward  integration  on  a  ray  is  continued  until  <,  some  specified  fraction  of 

00-  When  a  ray  “reaches  the  body,”  the  results  on  that  ray  are  “frozen”  and  the  ray 
is  dropped  from  the  computation.  Since  seven  rays  are  required  for  evaluation  of  the  x 
derivatives,  the  integration  stops  when  the  number  of  active  rays  becomes  less  than  seven. 
It  is  desirable  that  the  rays  drop  one  at  a  time  from  the  axis  outward,  but  frequently  they 
will  be  dropped  in  clumps  depending  on  step  sizes  and  the  degree  of  parallelism  between 
shock  and  body. 

3.1.5  Shock  Iteration  Procedure 

When  the  inviscid  flow  calculation  stops,  data  for  the  last  three  points  on  each  ray  are 
left  in  the  arrays  B0,  Bl,  and  B2. 

The  First  step  in  the  shock  iteration  procedure  is  to  linearly  extrapolate  the  remaining 
mass  flux  across  each  ray  to  zero.  This  defines  the  points  (zed(m),  r^fm))  at  the  surface  of 
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the  computed  body.  The  desired  body  is  then  translated  along  the  axis  of  symmetry  to  bring 
it  into  coincidence  with  the  point  z^  (MATCH),  red  (MATCH),  where  MATCH  is  the  first 
ray  off  the  axis.  The  intersection  of  each  ray  with  the  desired  body  is  designated  as  zB(m), 
rB(m). 


The  scheme  for  correcting  the  shock  shape  involves  translating  each  ray  along  its  length 
to  bring  the  point  (z^,  r^)  into  coincidence  with  zB,rB.  In  practice,  this  idea  is  expressed 
in  the  equations. 


z'Sh(m)  =  zSh(m)  +  x[l  -  Xx(m)][zB(m)  -  zed(m)] 
r'Sh(m)  =  rSh(m)  +  x[l  -  Xx(m)][zB(m)  -  z^m)] 


(46) 


Only  a  fraction  x[1.0  -  Xx]  of  the  z  and  r  errors  is  removed  at  each  iteration.  Appropriate 
values  for  x  and  X  have  been  found  to  be  x  -  0.35,  X  =  -0.5. 

Two  tests  for  convergence  of  the  shock  iteration  must  be  simultaneously  satisfied: 

1 .  The  RMS  deviation  between  extrapolated  points  and  body  points  must  be  less 
than  an  input  constant,  RMSTST,  which  is  defaulted  to  5.0  x  10~4. 

2.  The  extrapolated  flow  angle  error  at  the  junction  point  with  the  M.O.C. 
calculation  must  be  <,  10.0  x  RMSTST. 

If  these  criteria  have  not  been  met,  the  shock  shape  parameters  RSH,  SIP,  and  SI  of 
Eq.  (42)  for  the  next  iteration  are  determined  by  a  least-squares  fit  of  the  revised  shock 
coordinates  x^,,  rBh- 

3.1.6  Viscous  Interaction  Option 

When  the  viscous  interaction  option  is  invoked,  the  calculation  proceeds  as  above  except 
that  Zed.  r«d  are  determined  by  interpolating  to  the  boundary-layer  mass  flux  from  the 
previous  iteration  instead  of  to  zero.  After  zB,rB  have  been  determined,  the  boundary-layer 
edge  and  mass  flux  are  calculated.  The  mass  flux  is  used  to  re-evaluate  Zed.r«i  and  the  edge 
coordinates  are  substituted  for  the  body  coordinates  in  Eq.  (46). 
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3.2  TRANSITION  TO  M.O.C.  CODE 

When  BREGN  returns  to  MAIN  with  a  converged  shock  shape,  the  amount  of  data  on 
Unit  9  is  highly  variable  —  anywhere  from  20  to  1 ,000  y-steps.  Subroutine  BRIDGE  reads 
Unit  9  and  selects  NOUT  records  as  nearly  equally  spaced  in  mass  flux  as  possible  and  calls 
Subroutine  OUT  to  print  the  data  in  dimensional  form. 

Starting  the  M.O.C.  program  is  a  little  “touchy,”  especially  when  a  boundary  layer  is 
present.  The  supersonic  program  is  sensitive  to  any  inconsistency  between  the  body  shape 
and  the  flow  angle  determined  from  the  flow-field  calculation.  Such  inconsistencies  arise 
because  the  shock  iteration  process  satisfies  the  flow  tangency  condition  only  on  the  average. 
To  remove  or  reduce  these  inconsistencies,  BRIDGE  calls  Subroutine  BLAYST,  which  adjusts 
the  axial  position  and  the  scale  (nose  radius)  of  the  body  to  make  the  body  slope  match  the 
flow  angle.  The  changes  in  scale  and  position  are  small  if  a  reasonable  shock  fit  has  been 
obtained.  (Section  4.2.7) 

When  the  body  adjustments  are  completed,  the  flow-field  data  on  the  junction  ray  at 
the  selected  y-stations  are  written  to  Unit  2  for  transfer  to  the  M.O.C.  code. 

Finally,  the  flow-field  data  are  extrapolated  or  interpolated  to  the  boundary-layer  edge 
or  the  body  surface.  These  data  are  also  printed  by  Subroutine  OUT. 

3.3  THE  HYPERBOLIC  REGION 
3.3.1  General 

The  flow-field  calculation  in  the  supersonic  region  is  performed  by  the  method  of 
characteristics  using  the  so-called  reference-plane  form  which  is  sometimes  attributed  to  Hartree 
(Ref.  13).  This  procedure,  described  in  Section  3.3.3  below,  converts  the  hyperbolic  partial 
differential  equations  into  ordinary  differential  equations  —  the  “compatibility  equations” 
—  which  apply  along  the  Mach  lines  and  the  streamlines  in  the  flow. 

Chemical  kinetic  effects  appear  directly  in  the  compatibility  equations  (Appendix  A)  and 
also  indirectly  through  their  influence  on  temperature  and  molecular  weight.  Since  the  latter 
effects  are  secondary,  it  is  possible  to  converge  the  fluid  dynamic  properties  at  the  downstream 
point  while  holding  the  chemistry  frozen.  After  convergence  is  obtained,  the  downstream 
composition  is  calculated  by  integrating  the  species  conservation  equations  (A-27)  -  (A-29) 
along  the  streamline  through  the  point. 
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The  bow  shock  is  ‘Titled,”  i.e.,  the  infinitely  thin  wave  is  located  by  matching  the  pressure 
calculated  with  the  Rankine-Hugoniot  equations  with  the  pressure  calculated  by  applying 
a  compatibility  equation  along  a  characteristic  from  the  interior  of  the  shock  layer. 

3.3.2  The  Mesh 

3.3.2.1  General 

As  shown  in  Fig.  3,  the  M.O.C.  data  surfaces  —  often  referred  to  as  rays  in  the  text 
—  are  cones  (or  planes)  which  intersect  the  body  surface  in  a  right  angle.  When  the  input 
data  surface  is  prepared  by  subroutine  BRIDGE  as  is  usually  the  case,  the  first  “ray”  is  normal 
to  the  shock  rather  than  to  the  body  surface.  Then,  as  indicated  in  Fig.  3,  the  next  NROTAT 
rays  are  pivoted  about  the  foot  of  the  input  ray  to  adjust  to  the  body-normal  mesh.  The 
default  value  NROTAT  =  4  (Section  4.2.8)  is  usually  satisfactory. 


Figure  3.  Overall  geometry  for  the  afterbody  program. 

The  mesh  consists  of  points  at  the  intersection  of  each  ray  with  the  body  surface  (or  the 
boundary-layer  edge),  each  of  the  selected  streamlines,  and  the  shock  wave.  Ray  spacing  is 
discussed  in  Section  3.3.3.4.  Following  a  given  streamline  from  ray  to  ray  is  facilitated  by 
assigning  an  index  NUM  which  never  changes.  Since  the  boundary-layer  edge  is  not  a 
streamline,  it  is  assigned  the  special  value  NUM  =  0.  The  true  streamlines  are  numbered 
consecutively  outward,  beginning  with  NUM  =  2. 

Occasionally,  the  growing  boundary  layer  will  swallow  a  numbered  streamline,  and  the 
corresponding  value  for  NUM  will  simply  disappear  from  the  output. 

In  the  natural  order  of  events,  a  new  point  would  be  added  to  the  mesh  at  each  ray-shock 
intersection.  Because  of  the  continually  decreasing  curvature  of  the  shock,  these  added 
streamlines  would  be  tightly  clustered  near  the  shock  and  would  soon  overwhelm  the  available 
storage.  To  avoid  this  situation,  a  streamline  is  started  from  the  shock  only  if  the  increment 
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is  more  than  an  input  fraction  CSST  (Section  4.2.7)  of  the  spacing  between  the  two  existing 
streamlines  closest  to  the  shock. 

When  the  number  of  streamlines  threatens  to  overflow  the  existing  storage,  the  program 
drops  odd-numbered  streamlines. 

3.3 .2.2  Exceptions 

The  mesh  just  described  is  sometimes  interrupted  to  accommodate  body  slope 
discontinuities  such  as  expansion  and  compression  comers.  These  exceptions  will  be  explained 
in  Section  3.3.5. 

3.3.3  Marching 

3.3.3. 1  Interior  Point 

In  Fig.  4,  the  flow  field  is  completely  known  at  the  points  designated  by  solid  symbols 
on  rays  A  and  B.  The  flow  field  at  point  IT  on  ray  B  is  to  be  calculated  by  an  iterative  process. 
The  calculation  is  initialized  by  extending  the  streamline  through  IK  to  intersect  ray  B  and 
by  estimating  the  pressure  and  flow  angle  at  IT.  Using  the  composition  at  DC  and  the  pressures 
at  IT  and  IK,  Eqs.  (A-26)  and  (A-30)  are  integrated  along  IK-IT  to  determine  the  “frozen” 
temperature  and  velocity  at  IT.  The  Mach  angle  j2(IT)  is  calculated,  and  the  intersections 
IR  and  IS  of  the  characteristic  lines  with  ray  A  are  found.  The  fluid  dynamic  properties  at 
IR  and  IS  are  interpolated  linearly,  and  the  compatibility  equations  are  integrated  along  the 
characteristics  to  obtain  a  pair  of  algebraic  equations  which  are  then  solved  simultaneously 
for  the  pressure  and  flow  angle  at  IT. 

If  the  relative  pressure  change  from  the  previous  cycle  exceeds  the  allowed  tolerance,  IT 
is  relocated,  f(IT)  is  averaged  with  <A(IK),  and  new  values  for  temperature  and  velocity  at 
IT  are  computed.  In  general,  three  or  four  iterations  are  required  to  converge  the  solution 
at  an  interior  point  to  a  relative  error  of  1  x  10-6  in  pressure. 

When  the  solution  is  converged,  the  species  conservation  equations  are  integrated  along 
the  streamline  from  IK  to  IT  and  the  temperature  is  updated.  This  integration  is  performed 
with  the  Lomax  algorithm  described  in  Section  3.1.4.  If  a  concentration  becomes  negative 
at  any  phase  of  the  procedure,  the  marching  step  is  halved.  The  calculation  is  aborted  after 
eight  unsuccessful  attempts.  Although  several  applications  of  the  Lomax  algorithm  may  be 
required  to  get  from  IK  to  IT,  only  the  point  IT  is  stored. 
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Figure  4.  Schematic  of  geometry  for  calculation  of  an  internal  point. 

3.3.3. 2  Interpolating  Along  the  Boundary 

One  or  the  other  of  the  backward  characteristics  from  points  near  the  extremities  of  ray 
B  may  not  intersect  ray  A.  This  difficulty  is  overcome  by  calculating  the  shock  and  body 
points  before  starting  the  interior  points.  Two  shock  points  and  two  body  points  are  calculated 
to  improve  the  accuracy  of  the  interpolation. 

3.3.3. 3  Data  Storage 

To  economize  on  storage,  data  for  points  on  a  completed  ray  and  also  for  points  on  a 
ray  under  computation  are  stored  in  the  same  arrays.  At  the  beginning  of  a  marching  cycle, 
the  data  for  points  on  the  completed  ray  are  shifted  upwards  by  MDELAY  points,  where 
MDELAY  is  the  number  of  points  saved  for  interpolation.  / 

The  default  value  of  MDELAY  is  5.  A  larger  value  restricts  the  number  of  streamlines 
in  the  mesh  while  a  smaller  value  restricts  the  size  of  the  marching  step.  Since  the  first  nine 
positions  in  the  data  arrays  are  reserved  for  special  uses  such  as  characteristic  base  points 
and  intermediate  body  and  shock  points,  the  effective  number  of  streamlines  is  NMAX  where 
NMAX  =  NY  -  MDELAY  -9.  NY  is  currently  fixed  at  50  by  PARAMETER  statements. 
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3. 3.3 .4  Ray  Spacing 

The  ray  spacing  is  determined  so  that  the  points  IR  and  IS  which  limit  the  zone  of  influence 
of  IT  (Fig.  4)  will  fall  within  the  IK±MDELAY  known  points.  This  is  accomplished  by 
calculating  the  length  of  a  representative  characteristic  segment  using  a  length  on  ray  A  equal 
to  MDELAY/2  times  the  minimum  separation  of  any  two  points  on  the  ray  and  the  maximum 
Mach  and  flow  angles.  The  characteristic  length  is  then  multiplied  by  a  safety  factor 
ZETAMX  <  1. 

The  intersection  of  the  new  ray  with  the  shock  is  located  by  projecting  the  characteristic 
segment  onto  an  extension  of  the  known  shock. 

Finally,  an  iterative  process  is  used  to  drop  the  ray  from  the  known  shock  point 
perpendicular  to  the  body  surface. 

As  the  calculation  moves  downstream,  the  marching  step  size  increases  because  the  Mach 
angle  decreases  and  the  separation  between  streamlines  increases  as  the  density  decreases. 
The  user  can  control  the  increase  in  ray  spacing  through  his  choice  of  ZETAMX  and,  to 
a  lesser  degree,  by  limiting  the  number  of  points  normal  to  the  body,  NMAX,  the  number 
of  points  on  the  initial  ray,  NOUT,  and  the  rate  at  which  points  are  added  at  the  shock,  CSST. 

3.3.4  Viscous  Interaction 

Figure  5  shows  a  typical  situation  at  the  boundary-layer  edge.  At  the  point  marked  0, 
the  external  flow  and  the  boundary-layer  similarity  parameters,  s  and  n  (Section  2.0),  are 
known.  The  problem  is  to  locate  point  IT  and  to  compute  the  flow  quantities  there. 


Figure  5.  Geometry  of  boundary-layer  Interaction  procedure. 
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Since  the  point  B2  is  known,  specifying  a  value  for  8  locates  the  point  IT  on  the  ray. 
Estimate  a  value  for  the  flow  angle  <)>n  on  the  streamline  IR-IT  which  enters  the  growing 
boundary  layer  at  IT.  The  pressure,  composition,  etc.,  at  IR  may  be  determined  by 
interpolation  along  the  previous  ray.  The  Mach  angle,  0,  may  be  determined  by  integrating 
along  the  streamline,  and  one  application  of  the  compatibility  equations  along  IS  -  IT  yields 
the  pressure,  p(IT).  With  the  external  flow  conditions  known  at  IT,  one  can  solve  the  boundary- 
layer  equations  to  evaluate  5 '  and  the  boundary-layer  mass  flow,  0,  between  B2  and  IT. 
Interpolation  for  ft  along  the  initial  ray  locates  a  new  position  for  the  streamline  base  IR. 
Finally, 

4>'  =  tan-,[(riT  -  ^/(zit  -  Z|R]]  (47) 


Thus,  the  interaction  problem  is  reduced  to  the  solution  of  equations  of  the  form 

6'  -  5  =  f(M,ir,s)  -6  =  0  (48a) 

<t>'  -  <t>  =  g (6,0,x,s)  -  <t>  =  0  (48b) 

where  the  functions  f  and  g  are  known  implicitly.  Equations  (48a)  and  (48b)  are  solved 
simultaneously  by  a  Newton-Raphson  procedure. 

3.3.5  Body  Discontinuities 

If  there  is  no  boundary  layer,  the  M.O.C.  code  will  treat  a  body  with  one  or  more  expansion 
corners  and  one  sharp  compression  corner.  The  program  will  not  follow  the  coalesence  of 
Mach  waves  into  a  secondary  shock. 

At  a  discontinuity,  the  fluid  dynamic  variables  are  multivalued.  For  the  finite  reaction 
rate  models,  the  composition  is  frozen  for  calculations  at  the  discontinuity  but  is  allowed 
to  relax  downstream. 

3.3.5. 1  Prandtl-Meyer  Expansion 

Since  most  of  the  calculation  through  an  expansion  fan  is  identical  with  the  usual 
characteristics  procedure,  it  is  necessary  only  to  set  up  the  new  ray  directions,  (subroutine 
RAYB),  and  to  calculate  the  conditions  at  the  singular  point  on  the  body,  (subroutine 
PRANDTL).  An  iterative  procedure  is  used  to  vary  the  temperature  at  the  singularity  until 
the  Prandtl-Meyer  angle  computed  from  the  Mach  number  matches  the  required  turning  angle. 
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3.3.5. 2  Imbedded  Shock 

When  a  compression  corner  is  present,  the  mesh  is  altered  to  that  shown  in  Fig.  6.  Below 
the  secondary  shock,  the  rays  are  normal  to  the  flared  surface;  above  the  secondary  shock, 
the  rays  are  kept  parallel  to  the  corner  ray,  which  is  normal  to  the  original  body  at  that  point. 
In  the  printed  output,  there  is  a  double  point  at  the  secondary  shock.  If  the  identifying  number 
(Section  3.2.2. 1)  assigned  to  the  double  point  is  greater  than  1 ,000,  it  represents  a  streamline 
which  was  interpolated  to  better  define  the  secondary  shock  but  will  not  be  continued  into 
the  secondary  flow  field. 


Figure  6.  Computing  mesh  in  presence  of  secondary  shock. 

Although  the  presence  of  the  secondary  shock  complicates  the  logic  of  the  program,  no 
new  techniques  are  required.  The  imbedded  shock  is  located  first  and  then  treated  as  if  it 
were  the  bow  shock  for  computation  of  the  secondary  flow  field.  In  the  computation  of  the 
primary  flow  field,  the  front  face  of  the  secondary  shock  serves  as  a  virtual  body.  The  run 
will  abort  when  the  secondary  shock  intersects  the  primary  bow  wave. 

The  run  will  also  be  aborted  with  an  explanatory  message  if  an  abrupt  expansion  corner 
(NCORN  =  1)  is  encountered  following  a  flare.  However,  a  sharp  compression  corner  may 
be  followed  by  a  smooth  convex  surface,  as  indicated  in  Fig.  6. 

3. 3. 5. 3  Downstream  Continuation 

If  the  continuation  option  has  been  selected  (I BODY  >  NBODY),  CUNIFLOW  will  add 
a  tangent  sphere  to  the  prescribed  body  (Fig.  1)  and  switch  to  a  constant  pressure  boundary 
condition.  This  pressure,  PLIMIT,  may  be  input.  The  default  value  for  PLIMIT  is  the  last 
body  pressure  calculated. 


4.0  USING  CUNIFLOW 

Sections  2.0  and  3.0  have  presented  the  physical  models  and  some  of  the  numerical 
techniques  employed  in  the  CUNIFLOW  program.  The  application  of  the  code  will  now  be 
discussed.  Sections  4.1  through  4.3  present  the  information  necessary  to  set  up  the  input 
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and  to  execute  different  types  of  CUNIFLOW  calculations.  Interpretation  of  the  printed 
output  is  discussed  in  Section  4.4.  Section  4.S  presents  data  relevant  to  the  problem  of  defining 
the  “operating  envelope”  for  CUNIFLOW. 

4.1  SYSTEM  REQUIREMENTS 

CUNIFLOW  is  written  in  the  Cray  version  of  FORTRAN77.  The  Cray  extensions  to  the 
basic  FORTRAN  language  include  Namelist  input  and  comments  on  the  same  line  as  a 
statement.  The  source  code  contains  5,704  statements  filling  8,250  lines.  The  field  length 
of  the  current  version  on  the  AEDC  Cray  X/MP  12  is  approximately  250,000  64-bit  words. 

CUNIFLOW  will  run  in  single  precision  on  a  32-bit  machine,  although  with  reduced 
reliability;  that  is,  there  are  a  number  of  places  in  the  code  where  the  reduced  precision  might 
occasionally  lead  into  a  wrong  branch  or  cause  a  lack  of  convergence.  However,  the  code 
was  originally  developed  on  IBM  machinery  and  the  default  values  for  the  convergence  criteria 
have  been  set  from  years  of  experience  on  32-bit  machines. 

Table  !  lists  the  assignments  of  the  10  I/O  units  referenced  in  the  code. 

Table  1.  File  Usage 

Unit  Use 


1  Input  for  a  restart 

2  Scratch  file  used  in  BRIDGE 

3  Edge  conditions  for  input  to  BLIMP 

4  Input  nonuniform  free-stream  conditions 

5  System  input 

6  System  output 

7  Output  for  plotting 

8  Input  gas  property  and  reaction  rates 

9  Output  for  a  restart  in  nose  region 

10  Output  for  a  restart  in  M.O.C.  region 


4.2  INPUT 

Figure  7  is  a  schematic  of  typical  input.  The  first  record  of  the  input  file  contains  only 
one  entry  —  the  eight-character  run  number  NR  UN  which  will  appear  on  each  page  of  output. 
Enter  NRUN  with  an  “A”  format.  The  second  record  contains  an  80-character  comment 
or  description  of  the  run,  HEAD.  If  the  word  NONUNIFORM  appears  in  the  first  ten  positions 
(columns),  the  code  will  read  the  free-stream  conditions  from  Unit  4  and  bypass  Namelist 
UPSTREAM. 
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The  problem  —  free-stream  conditions,  body  shape  and  thermochemical  model  —  is 
specified  through  the  Namelists  UPSTREAM,  OPTIONS,  BODY,  SHOCK,  and  CONTROL 
(Tables  2,4  -  7),  which  are  read  in  subroutine  BBIN,  and  SUPER,  which  is  read  in  subroutine 
START. 


Figure  7.  Input  schematic. 

The  SHOCK  and  CONTROL  namelists  supply  initial  estimates  for  the  shock  shape  parameters 
and  values  for  convergence  criteria,  grid  spacing,  and  print  controls.  Namelists  RESTART 
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and  CHANGE  are  read  only  when  a  run  is  to  be  restarted  after  some  interruption.  Namelist 
BODY  is  read  once  for  each  body  segment. 

Thermodynamic  and  reaction  rate  data  are  read  from  a  formatted  file  assigned  to  Unit 
8.  For  equilibrium  cases.  Unit  8  is  referenced  in  subroutine  RGAS;  for  nonequilibrium  cases, 
Unit  8  is  read  in  subroutine  CHEMIN.  Since  the  thermochemical  data  are  rarely  changed, 
it  is  advantageous  to  maintain  a  semi-permanent  file  which  may  be  copied  and  edited  as 
required.  The  equilibrium  data  file  should  never  be  touched. 

4.2.1  Units 

The  user  may  choose  the  English  engineering  system  (feet,  slugs,  seconds),  the  “cgs”  system 
(centimeters,  grams,  seconds)  or  the  SI  system  (meters,  kilograms,  seconds)  for  input. 
Composition  data  may  be  input  as  mole-fractions,  as  mass-fractions,  or  as  moles  per  gram. 
Unless  otherwise  specified,  “cgs”  units  and  moles  per  gram  are  assumed.  Output  is  in  the 
cgs  system. 

Units  for  the  thermochemical  input  data  depend  on  the  value  of  MODEL  (cf . ,  Appendix  B). 

4.2.2  Free-Stream  Conditions 

The  flow  field  upstream  of  the  bow  shock  is  input  through  the  Namelist  UPSTREAM 
and/or  a  formatted  file  on  Unit  4.  The  entries  in  UPSTREAM  are  defined  in  Table  2. 

Table  2.  Variables  in  Namelist  UPSTREAM 


Array  Default 


Item 

Length 

Value 

Definition 

ALT 

... 

Geometric  altitude 

DENS 

— 

Free-stream  density 

PRES 

— 

Free-stream  pressure 

VEL 

— 

Centerline  axial  velocity 

FRACTN 

30 

— 

Chemical  composition 

EPS 

2 

EPS(I)  is  vibrational  energy  of  species 
OSC(I).  Use  erg/gm  or  ft-lb/slug  or  J/kg 

OSC 

2 

Name  of  oscillating  species 

TEMP 

— 

Free-stream  temperature 

SPEC 

30 

... 

Names  of  species  in  free  stream.  Enclose 
in  single  quotes  and  left-adjust. 
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For  a  typical  reentry  problem,  the  only  required  entries  are  the  velocity,  VEL,  and  the 
altitude,  ALT  —  if  cgs  units  are  acceptable. 

If  the  airflow  about  a  wind  tunnel  or  ballistic  range  model  is  to  be  investigated,  one  should 
specify  IUNIT  =  1  and  input  values  for  the  velocity  VEL  and  for  two  of  the  three  state 
variables  PRES,  TEMP,  DENS.  Also  specify  the  units,  the  composition,  and  vibrational 
energies  EPS  if  not  in  equilibrium.  If  units  other  than  moles/gm  are  used,  specify  by  selecting 
either  ‘MOLE’  for  mole-fraction  or  ‘MASS’  for  mass-fraction  as  the  value  for  FRAC.  The 
species  names  are  entered  into  the  SPEC  array,  and  the  identifiers  for  species  with 
nonequilibrium  vibration  are  entered  into  the  OSC  array.  These  names  are  to  be  left-adjusted 
and  enclosed  in  single  quotes.  The  run  will  stop  with  a  diagnostic  message  if  the  input  species 
name  has  no  match  in  the  user-supplied  thermochemical  data. 

4.2.3  Nonuniform  Free  Stream 

If  the  free  stream  is  spatially  nonuniform,  the  comment  card  HEAD  must  begin  with 
the  keyword  NONUNIFORM  in  column  1,  IUNIT  in  OPTIONS  namelist  must  equal  2,  and 
a  table  of  aerodynamic  and  composition  variables  as  functions  of  position  must  be  assigned 
to  Unit  4.  The  position  is  specified  in  *'cgs”  units  using  the  shock -based  cylindrical  coordinate 
system  of  Fig.  2.  Composition  is  assumed  to  be  in  mole  fractions. 

The  format  of  the  free-stream  file  is  shown  in  Table  3.  The  first  few  records  define  the 
size  and  contents  of  the  file,  including  the  names  of  the  free-stream  species  and  the  two  species 
which  may  be  in  vibrational  nonequilibrium.  The  flow  field  is  to  be  input  with  the  radius 
increasing  at  each  axial  station.  Centerline  data  must  always  be  included. 

Table  3.  Input  Format  for  Nonuniform  Free-Stream  Flow 


Card 

Column 

Format 

Item 

Description 

I 

5 

15 

NZZ 

Number  of  axial  stations 

10 

15 

NRR 

Number  of  radial  stations 

15 

15 

NSP 

Number  of  species  <,  30 

20 

15 

NVIB 

Number  of  nonequilibrium 
vibrators  <  2 

2 

1-4 

A4 

SPEC(l) 

Name  of  free-stream  species 

5-8 

A4 

SPEC(2) 

Name  of  free-stream  species 

9 

— 

— 

— 

— 

A4 

SPEC(NSP) 

Name  of  free-stream  species 
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Table  3.  Concluded 

Card 

Column 

Format 

Item 

Description 

2  (Contd.) 

— 

A4 

OSC(l) 

Name  of  vibrator 

— 

— 

OSC(NVIB) 

Name  of  vibrator 

3 

1-10 

E10.4 

ZSPOT 

Z  coordinate  from  shock  vertex 

11-21 

E10.4 

RSPOT 

R  coordinate  from  centerline 

21-30 

E10.4 

VINF 

Radial  velocity  component 

31-40 

E10.4 

PRES 

Pressure,  dyne/cm2 

41-50 

E10.4 

DENS 

Density,  gm/cm3 

51-60 

E10.4 

UINF 

Axial  velocity  component,  cm/sec 

61-70 

E10.4 

FRACTN(l) 

Mole-fraction  of  species  1 

71-80 

E10.4 

FRACTN(2) 

Mole-fraction  of  species  2 

E10.4 

FRACTN(NSP) 

Mole-fraction  of  species  NSP 

E10.4 

EPSINF(l) 

Vibrational  energy,  erg/gm 

E10.4 

EPSINF(NVIB) 

Vibrational  energy,  erg/gm 

Notes: 

1 .  Card  2  is  read  only  if  NSP  >  0. 

2.  NSP  is  number  of  species  in  free  stream,  not  number  in  model. 

3.  VINF  must  be  zero  on  the  centerline. 

4.  For  each  axial  station,  NRR  sets  of  card  3  are  needed.  The  centerline 
(RSPOT  =  0)  must  be  included  at  each  axial  station. 

4.2.4  Namelist  OPTIONS 

Table  4  defines  the  variables  in  Namelist  OPTIONS.  The  options  for  including  sharp 
corners  and  for  the  downstream  continuation  will  be  discussed  in  connection  with  Namelist 
SUPER  (Table  8). 


Table  4.  Variables  in  Namelist  OPTIONS 


Item 

Type 

Default 

Function 

FRAC 

Character 

‘M/GM’ 

Composition  unit.  Choose  ‘MASS’  or  ‘MOLE’  or 

‘M/GM’ 

1NAFT 

Integer 

1 

0  Do  only  nose  region 

1  Do  both  regions 

2  Do  only  aft  region 


35 


AEDC-TR-88-31 


Table  4.  Concluded 


Item 

Type 

Default 

Function 

INBLR 

Integer 

0 

0  No  viscous  interaction 

1  Include  viscous  interaction 

1NDIM 

.  Integer 

1 

0  Planar  (2-D)  flow 

1  Axisymmetric  flow 

IUNIT 

Integer 

0 

0  standard  atmosphere;  supply  only  VEL  and 

ALT 

1  Uniform  free  stream;  supply  composition 
and  any  2  state  variables 

2  Nonuniform  free  stream.  Use  Unit  4 

MODEL 

Integer 

0 

0  Nonequilibrium  model  Sec.  2.2.1 

1  Perfect  gas  Sec.  2.2.2 

2  Equilibrium  air  Sec.  2.2.3 

3  Nonequilibrium  Sec.  2.2.4 

4  Nonequilibrium  Sec.  2.2.5 

NEW 

Logical 

li'J’D 

Set  NEW  =  “F”  for  a  restart 

•NVIB 

Integer 

0 

Number  of  nonequilibrium  vibrators,  NVIB  2 

UNIT 

Character 

*CGS’ 

Input  unit  system.  Use  ‘ENG’  for  ft-slug-sec 
system;  ‘CGS’  for  cm-gm-sec  system;  ‘MKS*  for 

SI  system 

NVIB  =  0  unless  MODEL 

=  0 

4.2.5  Namelist  SHOCK 

Table  5  defines  the  three  initial  shock  parameters  input  through  Namelist  SHOCK.  No 
default  values  are  provided.  The  initial  estimate  for  the  shock  radius  of  curvature  RSHP 
must  be  oon-zero  and  in  the  units  specified  in  UPSTREAM.  Initial  non-zero  values  for  the 
other  two  parameters,  SI  and  SIP,  are  useful  but  not  essential,  A  plot  of  shock  radii  RSH 
normalized  by  the  nose  radius  for  a  number  of  converged  nonequilibrium  solutions  is  shown 
in  Fig.  8.  The  significance  of  the  parameter  is  discussed  in  Section  4.5.  Table  E-l  of 
Appendix  E  presents  more  detail  as  well  as  the  corresponding  values  of  SI  and  SIP. 

For  equilibrium  air  (MODEL  =  2),  the  tabulated  results  of  Inouye  (Ref.  32)  are  helpful. 
(SIP  corresponds  to  Inouye’s  A7  and  SI  to  2A7/-^(M2«,  —  1).  The  equilibrium  air  shock 
parameter  values  are  also  reasonable  estimates  for  nonequilibrium.  For  flow  conditions  not 
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encompassed  by  the  tables  of  Ref.  32,  use  RSHP  =  1.1  x  RN,  SIP  *  0.1  and  SI  *  0.01 
where  Rn  is  the  centerline  radius  of  curvature  of  the  nose. 

Table  5.  Parameters  in  Namelist  SHOCK 
Item  Default  Definition 

RSHP  —  Dimensional  radius  of  curvature  at  centerline. 

SIP  ---  Coefficient  of  r+  in  Eq.  (42) 

SI  —  Coefficient  of  r3  in  denominator  of  Eq.  (42) 

An  initial  shock  which  is  too  blunt  (too  open)  for  the  flow  conditions  will  occasionally 
cause  a  run  to  abort  with  diagnostic  messages  indicating  negative  species  or  temperatures 
on  the  higher  numbered  rays.  When  this  is  the  case,  either  make  the  initial  shock  “tighter” 
(by  decreasing  RSHP  or  increasing  SIP)  or  else  decrease  the  ray  spacing  DELX. 


Figure  8.  Converged  values  for  radius  of  curvature  of  nonequilibrium  shock. 
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4.2.6  Namelist  BODY 

Table  6  defines  the  entries  in  Namelist  BODY: 

Table  6.  Variables  in  Namelist  BODY 

Description 


C0  in  Eq.  (49) 


C4  in  Eq.  (49) 

If  “T”  another  segment  follows.  Read  BODY  again 
Exponent  *‘P”  of  radius  in  Eq.  (49) 

Radius  at  junction  point  “rj”,  Eq.  (49) 

Exponent  of  axial  coordinate  “q”,  Eq.  (49) 

Axial  coordinate  of  junction  point  between  body  segments. 

“zj”  in  Eq.  (49) 

The  body  is  assumed  to  consist  of  NBDY  segments  (NBDY  <  6).  Each  segment  is  described 
by  an  equation  of  the  form 

(r  -  rj)p  -  (z  -  Zj)<i  [C0  +  C|(z  -  zj)  +  C2(z  -  Zj)2...  +  C4(z  -  Zj)4] 

(49) 

where  z  and  Zj  are  measured  from  the  nose  of  the  body.  The  code  will  adjust  to  the  shock- 
based  coordinates  used  internally.  The  first  equation  must  represent  a  sphere,  i.e., 

r2  -  2zR  -  z2  (50) 

At  the  junction  points  zj.tj  the  radius,  r,  must  be  continuous. 

BODY  is  read  once  for  each  segment.  The  logical  flag  MORE  is  used  to  indicate  that 
another  segment  follows.  The  segment  equations  are  input  in  dimensional  form  using  the 
system  specified  by  UNIT  (Table  2). 

4.2.7  Namelist  CONTROL 

Table  7  defines  the  entries  in  Namelist  CONTROL: 


Variable  Default 


CO 

Cl 

C2 

C3 

C4 

MORE 

REXP 

RJUNC 

ZEXP 

ZJUNC 
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Table  7.  Entries  in  Namelist  CONTROL 


Item 

Default 

Definition 

DELMAX 

5.0  x  10-3 

Maximum  step  in  y-direction 

DELMIN 

1.0  X  10-6 

Minimum  step  in  y-direction 

DELX 

0.04 

Ray  spacing 

ERROR 

1.0  x  lO"6 

Convergence  criterion  for  all  iterations  except  shock 

GAMTST 

1.0  x  10~7 

Minimum  mole-fraction  for  inclusion  in  source  term 

IPRNT 

0 

If  non-zero,  print  status  line  every  IPRNT  steps 

ISW6 

0 

Print  diagnostic  information  if  non-zero 

JOIN 

NXDIV-7 

Junction  ray  to  M.O.C.  code 

METHOD 

6 

Use  4,  6,  8,  10  to  select  integration  algorithm 

NBSIT 

1 

Number  of  bow  shock  iterations 

NOUT 

15 

Number  of  points  saved  for  M.O.C.  code;  also  number 
of  full  printout  stations 

NXDIV 

20 

Number  of  DELX  increments  <21 

PSIRAT 

0.10 

Flux  ratio  to  stop  integration 

RMSTST 

5.0  x  10“4 

Convergence  criterion  for  shock  iteration 

SKAP 

0.35 

Damping  factor  in  Eq.  (46) 

SLAM 

-0.50 

Damping  factor  in  Eq.  (46) 

TESTM 

1.10 

Minimum  Mach  number  on  junction  ray 

TIME 

28 

Maximum  run  time  in  seconds 

TTEST 

0.15 

Maximum  relative  change  in  temperature  over  a  y-step 

XO 

0.04 

Position  of  first  non-axis  ray 

The  default  values  which  are  provided  for  the  parameters  in  the  CONTROL  namelist 
are  adequate  for  nonequilibrium  problems  in  the  central  part  of  the  operating  envelope.  Section 
4.5.  The  value  for  TIME  is  based  on  experience  with  the  Cray  X/MP12  and  is  consistent 
with  NBSIT  =  1,  a  chemical  model  with  10  reactions  and  a  body  length  of  approximately 
15  nose  radii.  A  single  iteration  of  the  bow  shock  shape  is  rarely  sufficient  unless  the  problem 
amounts  to  a  small  perturbation  of  an  existing  solution  with  accurately  known  shock 
parameters. 

The  shock  iteration  process  is  regulated  by  the  controls  RMSTST,  SKAP,  SLAM,  and 
PSIRAT,  as  well  as  NBSIT.  The  default  value  for  RMSTST  is  probably  too  lax  since  root- 
mean-square  deviations  between  desired  and  calculated  body  points  <,  10  ~ 4  are  readily 
obtained.  The  control  for  ending  the  integration  of  the  flow  equations  is  PSIRAT,  which 
is  defaulted  to  0.1  to  avoid  the  region  of  “stiff’  chemistry  near  the  body  surface.  If  inspection 
of  the  status  reports  indicates  few  problems  near  the  surface,  PSIRAT  might  well  be  decreased 
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to  0.05  or  even  0.03  to  increase  the  precision  of  the  shock  iteration  process.  The  convergence 
rate  of  the  shock  iteration  process  is  affected  by  the  damping  factor  SKAP  which  has  been 
set  conservatively.  If  RMSDEV  increases  or  oscillates  from  iteration  to  iteration,  decrease 
SKAP.  Shock  convergence  has  been  obtained  with  SKAP  as  large  as  0.6.  There  is  no  reason 
to  change  SLAM  unless  MAXDEV  begins  to  increase. 

The  forebody  computational  grid  is  regulated  by  the  controls  DELMAX,  DELMIN, 
DELX,  NXDIV,  and  X0.  The  grid  consists  of  NXDIV  evenly  spaced  rays  normal  to  the  shock. 
The  spacing  between  the  rays  is  DELX  and  the  first  ray  is  located  at  X0.  If  X0  ^  DELX, 
the  axis  will  be  included  in  the  mesh  and  NXDIV  will  be  increased  by  4  to  account  for  the 
axis  and  the  three  “virtual”  rays  which  allow  the  use  of  the  seven-point  differentiation 
algorithm  at  the  axis.  One  might  choose  X0  >  DELX  in  a  situation  where  the  species 
conservation  equations  were  very  stiff  near  the  stagnation  point  and  interest  was  in  the 
afterbody  flow  field. 

The  marching  step  in  the  “y”  direction  is  internally  controlled  between  the  limits  DELMIN 
and  DELMAX.  At  the  shock,  CUNIFLOW  will  try  an  initial  integration  step 

y  =  100  X  DELMIN 

and  will  halve  this  value  four  times  or  until  the  integration  is  successful.  The  default  value 
for  DELMIN  is  adequate  for  most  equilibrium  and  many  nonequilibrium  problems,  but  a 
smaller  value  will  be  required  for  exceptionally  strong  shocks  or  stiff  chemistry.  Too  small 
a  value  for  DELMIN  may  cause  unacceptable  run  times.  Since  the  nondimensional  shock 
layer  thickness  is  *  0.05,  the  default  value  for  DELMAX  is  probably  the  largest  step  which 
will  give  adequate  flow-field  resolution. 

The  integration  algorithm  is  selected  by  METHOD  which  may  assume  the  values  4,  6, 
8, 10.  The  stability  of  the  Lomax  procedure  and  the  number  of  derivative  evaluations  increase 
with  METHOD.  For  equilibrium  air  problems,  the  standard  Runge-Kutta  algorithm, 
METHOD  =  4,  works  well.  For  many  nonequilibrium  problems  it  is  necessary  to  use 
METHOD  =  10. 

GAMTST  represents  the  smallest  mole-fraction  which  will  be  used  in  computing  the  source 
terms  in  the  species  conservation  equations.  Also,  if  y\  <  GAMTST,  the  code  sets  dy/dy 
=  max  (dyj/dy.0).  Note  that  y,  is  left  unchanged.  The  foregoing  GAMTST  artifice  is  useful 
where  the  stream  contains  small  amounts  of  some  species  such  as  02,  which  will  continue 
to  dissociate  at  the  prevailing  temperatures  such  as  those  behind  a  strong  shock  or  near  the 
stagnation  point.  Since  GAMTST  perturbs  the  composition,  the  minimum  value  of  GAMTST 
which  will  allow  the  run  to  finish  should  be  selected.  Trouble  will  also  be  experienced  if  0 
<  GAMTST  <  machine  zero  (—  10-14  for  Cray). 
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Output  from  the  nose  region  part  of  the  code  is  controlled  by  IPRNT,  NOUT,  and  ISW6. 
If  IPRNT  is  non-zero,  a  one-line  status  report  is  printed  every  IPRNT  successful  integration 
steps.  A  non-zero  value  for  ISW6  produces  diagnostic  information  from  certain  subroutines. 
The  amount  of  data  can  be  overwhelming  so  that  care  should  be  taken  when  this  option 
is  used.  NOUT  is  the  number  of  records  to  be  printed  and  also  the  number  of  points  on 
the  starting  line  for  the  Method-of-Characteristics  (afterbody)  code.  NOUT  <  NMAX,  where 
NMAX  (see  Table  8)  is  the  maximum  number  of  points  per  afterbody  ray. 

If  the  minimum  Mach  number  on  the  junction  ray  (designated  by  ray  index  M  =  JOIN) 
<  TESTM,  the  nose  region  code  will  increase  DELX  and  XO  on  the  next  shock  iteration. 
The  default  value  for  TESTM  is  as  low  as  is  safe,  typically  1.10.  The  maximum  allowable 
value  for  JOIN  is  NXDIV-6.  The  default  values  for  JOIN,  XO,  and  DELX  put  the  junction 
line  at  40  to  45  deg  from  the  axis. 

4.2.8  Input  to  the  Hyperbolic  Flow  Code 

The  controls  for  the  Method-of-Characteristics  code  are  input  through  Namelist  SUPER, 
which  is  read  by  subroutine  START.  The  initial  data  surface  is  usually  read  from  a  file  on 
Unit  2,  but  if  INAFT  =  2  (Table  4),  the  starting  line  is  read  from  the  system  input  (Unit 
5).  The  entries  in  Namelist  SUPER  are  defined  in  Table  8. 

'  Table  8.  Controls  in  Namelist  SUPER 


Item 

Default 

Definition 

CSST 

1.0 

Controls  streamline  spacing  at  shock 

IBODY 

0 

Downstream  continuation  (Sec.  2.1)  if  IBODY  >  number 
of  input  body  segments 

.  IPOLY 

1 

Current  body  equation 

ION 

1 

Plasma  properties  calculated  if  ION  >  0 

MAXRAY 

0 

Maximum  number  of  rays  to  be  calculated 

METHOD 

6 

Algorithm  order.  Will  equal  blunt  body  value  unless  otherwise 
specified 

MODRAY 

1 

Print  ray  if  MOD  (NRAY,MODRAY)  =  NRAOUT 

NRAOUT 

0 

Print  ray  if  MOD  (NRAY.MODRAY)  =  NRAOUT 

NCHI 

0 

ERROR  =  10NCH1  if  NCHI  not  =  0 

NCORN(6) 

6x0 

0  This  junction  is  smooth 

1  This  junction  is  expansion  corner 

2  This  junction  is  compression  corner 

NFIL1 

1 

Write  restart  record  every  NFIL1  rays 

NITER 

10 

Maximum  number  of  iterations 

NMAX 

35 

Maximum  number  of  points  on  ray.  NMAX  <  35 
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Table  8.  Concluded 


Item 

Default 

Definition 

NMOD4 

1 

Print  diagnostics  from  certain  subroutines  if 

NREM4 

2 

MOD(NRAY,NMOD4)  =  NREM4 

NMOD6 

1 

Print  diagnostics  from  certain  routines  if 

NREM6 

2 

MOD(NRA  Y ,  NMOD6)  =  NREM6 

NPRROT 

1 

Number  of  rays  in  Prandtl-Meyer  expansion 

NROTAT 

4 

Number  of  rays  used  for  initialization 

PLIMIT 

— 

Fixed  boundary  pressure  for  continuation  option 

TIME 

— 

Time  limit  default  =  value  from  blunt  body  code 

TSTGAM 

GAMTST 

Same  role  as  GAMTST 

ZMAX 

— 

Axial  station  to  end  calculation 

ZETAMAX 

0.8 

■  Safety  factor  in  ray  spacing  algorithm.  ZETAMX  <  1 .0 

The  spacing  between  data  surfaces  grows  as  the  calculation  progresses  downstream  because 
the  local  Mach  number  increases.  The  rate  of  growth  is  regulated  by  ZETAMX  but  is  also 
affected  by  the  streamline  spacing.  A  streamline  is  added  at  the  shock  when  the  distance 
between  the  shock  point  and  the  first  streamline  below  the  shock  exceeds  CSST  times  the 
separation  between  that  streamline  and  the  next  lower.  When  the  number  of  streamlines  exceeds 
NMAX,  approximately  half  of  the  streamlines  are  removed  from  the  grid  with  the  effect 
that  the  ray  spacing  is  approximately  doubled  from  that  station  downstream.  The  initial  ray 
spacing  is  also  affected  by  NOUT,  the  number  of  points  on  the  initial  line.  Halving  ZETAMX 
triples,  at  least,  the  number  of  rays  and  the  CPU  time  for  a  run.  The  practical  limits  for 
ZETAMX  are  roughly  0.2  5  ZETAMX  <  0.9.  The  default  value  CSST  =  1.0  has  proven 
to  be  optimum  for  80  percent  of  the  calculations  made  with  CUNIFLOW.  Experience  has 
shown  that  any  adjustments  should  be  kept  in  the  range  0.9  <  CSST  s  1.2. 

The  maximum  value  for  NMAX  is  NY  -  15;  NY  is  set  by  PARAMETER  statements  in 
the  source  program.  The  user  may  specify  either  MAXRAY,  the  maximum  number  of  rays 
to  be  calculated,  or  the  final  axial  station  ZMAX,  but  not  both.  ZMAX  is  entered  as  a 
dimensional  number  in  the  same  units  as  RSHP.  If  both  ZMAX  and  MAXRAY  are  entered, 
the  code  will  use  MAXRAY. 

The  default  value  for  MODRAY,  the  number  of  rays  to  skip  between  printouts,  is  less 
than  desirable  in  most  cases;  MODRAY  =  10  is  frequently  fine  enough.  The  data  on  the 
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initial  ray,  on  ray  NROTAT,  on  rays  at  each  comer,  and  on  rays  where  streamlines  are  dropped 
are  always  printed.  Default  values  for  NMOD6,  NREM6,  NMOD4,  and  NREM4  have  been 
selected  to  prevent  the  printing  of  the  auxiliary  output  from  each  iteration  at  each  point, 
since  the  amount  of  this  output  would  be  overwhelming. 

To  invoke  the  downstream  continuation  option,  the  user  must  input  a  nonzero  value  for 
PLIMIT  and  set  IBODY  equal  to  the  number  of  input  body  equations.  If  IBODY  is  greater 
than  the  number  of  equations,  the  code  will  construct  an  aft-facing  tangent  sphere  as  shown 
in  Fig.  1.  Whether  input  or  internally  created,  segment  IBODY  serves  only  to  define  the  ray 
geometry. 

4.2.9  Restarting 

The  blunt  body  code  writes  a  restart  record  on  Unit  9  after  each  successful  step.  The 
current  number  of  records  is  printed  every  IPRNT  steps  as  part  of  the  status  report  and  also 
when  the  body  is  reached.  Unit  9  is  rewound  at  the  beginning  of  each  shock  iteration.  The 
M.O.C.  code  writes  a  restart  record  to  Unit  10  and  prints  the  number  of  records  after  each 
NFIL1  rays. 

To  restart  in  either  code  segment,  assign  the  appropriate  file  to  Unit  1,  set  NEW  -  .FALSE, 
in  Namelist  OPTIONS,  and  add  Namelists  RESTART  and  CHANGE  as  shown  in  Fig.  7. 
Table  9  defines  the  entries  in  RESTART: 


Table  9.  Entries  in  Namelist  RESTART 


Entry 

Default 

Definition 

NRUN 

Old  run  ident. 

Run  identifier  can  be  changed 

ITAPE 

0 

1 - Restarts  in  blunt  body  code 

2- Restarts  in  M.O.C.  code 

NRCORD 

0 

Resume  computation  after 
NRCORD 

There  is  no  record  on  either  Unit  9  or  10  with  NRCORD  =  0.  It  is  always  wise  to  set 
NRCORD  less  than  the  maximum  value  printed  in  the  status  report. 

Most  of  the  blunt  body  controls  and  all  of  the  M.O.C.  controls  may  be  changed  during 
a  restart.  Also,  any  part  of  the  body  definition  which  has  not  been  used  may  be  altered  safely. 
Table  10  lists  the  blunt-body  controls  which  may  be  changed. 
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Table  10.  Namelist  CHANGE 


Entry 

Default 

Definition 

DELMAX 

5.0  X  10-3 

Maximum  step  in  y-direction 

DELMIN 

1.0  x  10'6 

Minimum  step  in  y-direction 

ERROR 

1.0  X  lO  6 

Convergence  criterion  for  all  iterations  except  shock 

GAMTST 

1.0  x  10“7 

Minimum  mole-fraction  for  inclusion  in  source  term 

INAFT 

1 

Change  1  to  0  or  0  to  1  if  desired;  do  not  use 

INAFT  —  2  on  restart 

IPRNT 

0 

If  non-zero,  print  status  line  every  IPRNT  steps 

ISW6 

0 

Print  diagnostic  info  if  non-zero 

JOIN 

NXDIV-7 

Junction  ray  to  M.O.C.  code 

METHOD 

6 

Use  4,  6,  8,  10  to  select  integration  algorithm 

NBSIT 

1 

Number  of  bow  shock  iterations 

NOUT 

15 

Number  of  points  saved  for  M.O.C.  code;  also  number  of 
full  printout  stations 

PSIRAT 

0.10 

Flux  ratio  to  stop  integration 

RMSTST 

5.0  x  10-4 

Convergence  criterion  for  shock  iteration 

SKAP 

0.35 

Damping  factor  in  Eq.  (46) 

SLAM 

-0.50 

Damping  factor  in  Eq.  (46) 

TESTM 

1.10 

Minimum  Mach  number  on  junction  ray 

TIME 

28 

Maximum  run  time  in  seconds 

TTEST 

0.15 

Maximum  relative  change  in  temperature  over  a  y-step 

4.3  OUTPUT 

This  section  describes  the  output  produced  during  the  course  of  a  completely  successful 
run.  Examples  are  contained  in  Appendix  F.  Diagnostic  messages  are  interpreted  in  Section  4.4. 

The  first  page  begins  with  a  largely  self-explanatory  "echo”  of  the  input.  Units  have  been 
converted  to  the  cgs  system  and  concentrations  are  expressed  in  moles  per  gram.  Only  the 
centerline  values_of  the  free-stream  parameters  are  printed.  XLMBDA  is  the  dimensionless 
combination  U2a><^(x/<^Tm  which  serves  as  an  inverse  gas  constant  in  the  nondimensionalized 
equations.  HINF  is  the  dimensionless  total  enthalpy  and  BIGAM  is  the  isentropic  exponent 
in  the  free  stream. 

The  quantities  AK(I),  BK(1),  CK(I),  and  DK(I)  printed  in  connection  with  each  reaction 
are  the  rate  parameters  defined  in  Eq.  (21).  The  reaction  type  is  defined  by  Eq.  (22).  The 
matrix  elements  couple  the  ith  reaction  to  the  oscillator  when  nonequilibrium  vibration 
is  specified. 
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4.3.1  Blunt  Body  Output 

The  blunt  body  output  begins  with  a  tabulation  of  the  current  shock  parameters  and 
geometry.  Since  this  table  is  printed  before  the  optional  rearrangement  for  including  the  axis 
of  symmetry  is  performed,  the  first  ray  (M  =  1)  will  be  off  the  centerline  by  XO.  The  shock 
angle,  SIGMA,  is  in  radians,  but  the  local  curvature  CURV  has  been  normalized  by  its 
centerline  value  1 .0/RSHP. 

As  the  computation  marches  inward,  a  progress  message  is  printed  every  IPRNT  steps. 
The  information  presented  consists  of  the  number  of  records  on  Unit  9  (NRCORD);  the 
number  of  failures  (FLOP  -  1);  the  index  of  the  lowest  numbered  ray  still  in  the  computing 
grid,  M INRAY;  the  last  three  y-values,  and  the  last  y-increment.  Since  the  dimensionless 
nonequilibrium  shock  standoff  distance  is  of  the  order  of  0.05  on  the  stagnation  line,  the 
last  y  printed  gives  a  crude  estimate  of  relative  position  within  the  shock  layer. 

If  ISW6  =  2,  the  quantities  involved  in  the  shock  iteration  process  (Section  3.1.5)  are 
printed  for  each  ray;  otherwise,  only  the  rms  deviation  between  extrapolated  and  desired 
body  points,  the  flow-angle  error  at  the  junction  point,  the  normalized  stagnation  line  standoff 
distance  (ZSTAG),  and  the  current  and  proposed  shock  parameters  are  output. 

When  convergence  has  been  obtained  or  the  allowable  number  of  global  iterations  has 
been  performed,  the  nondimensional  body  coordinates  ZB,  RB,  the  body  angle  in  radians, 
the  boundary-layer  thickness,  and  boundary-layer  mass  flux  are  printed. 

Each  page  of  the  blunt  body  flow-field  data  corresponds  to  a  constant  y  value  which 
is  printed  in  the  upper  right  corner.  If  the  stagnation  line  was  included  in  the  computation, 
it  is  printed  as  ray  number  4.  All  data  are  dimensional  except  for  the  Mach  number,  the 
equilibrium  indicator,  labeled  AVER. CHI,  and  the  ratio  of  the  mass  flux  across  the  ray  between 
the  given  y-station  and  the  body  to  the  mass  flux  into  the  shock  between  the  ray  and  the 
stagnation  line.  The  equilibrium  indicator  AVER. CHI  is  defined  as  [E^I/Nr  where  Xj  is 
defined  by  Eq.  (24).  The  flow  angle  is  defined  as  the  angle  (in  degrees)  between  the  z-axis 
and  the  resultant  velocity  VTOT. 

The  flow-field  printout  includes  only  rays  MINRAY  through  JOIN.  The  data  on  rays 
JOIN  +  1  through  NXDIV  are  of  lesser  quality  than  the  data  on  the  lower  numbered  rays 
because  of  the  inadequacies  of  the  x-differentiation  algorithm,  and  will  be  replaced  with 
M.O.C.  data.  As  the  integration  on  a  given  ray  reaches  the  body  surface  or  boundary-layer 
edge,  the  data  on  that  ray  stop  changing  and  the  ray  is  dropped  from  the  printout  to  prevent 
confusion. 

Following  the  printout  of  the  integrated  flow  field,  the  dimensionless  body  shape  parameters 
corrected  for  the  latest  shock  radius  are  printed  before  and  after  the  adjustment  process 
described  in  Section  3.2  is  performed.  The  last  page  of  the  blunt-body  printout  represents 
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the  results  of  a  linear  extrapolation  of  all  flow-field  data  to  the  body  surface  or  boundary- 
layer  edge. 

4.3.2  Plotting  File 

At  each  y-station,  the  dimensional  flow-field  properties  are  written  in  unformatted  form 
to  a  file  assigned  to  Unit  7.  The  first  record  contains,  in  order,  MINRAY,  JOIN,  NYIB, 
and  CSIG,  where  CSIG  (N)  =  cos  (>„)■  Each  of  the  following  JOIN  -  MINRAY  +  1 
records  contains  the  data  for  a  ray.  The  data  entries  are  z,r,P,T,RHO,XMW,VTOT,  flow 
angle,  total  enthalpy,  entropy,  average  chi,  electron  density,  Mach  number,  flux  ratio,  specific 
radiation  rate,  mass  flux,  and  species  mole-mass  ratios  (mole/gm).  Note  that  the  entries  are 
arranged  in  the  same  order  as  on  the  printout. 

4.3.3  M.O.C.  Output 

Printout  from  the  Method-of-Characteristics  code  begins  with  an  echo  of  the  input  Namelist 
SUPER.  ZMAX  and  PLIMIT  will  be  dimensional  (cgs  units).  The  principal  control  parameters 
are  then  tabulated  in  nondimensional  form. 

The  following  output  data  are  arranged  so  that  each  of  the  selected  rays  starts  a  new 
page.  The  first  block  labeled  with  NPT  =  0  is  either  at  the  body  surface  or  at  the  boundary- 
layer  edge,  which  is  not,  of  course,  a  streamline.  The  last  block  presents  the  flow  data 
immediately  behind  the  shock.  The  other  blocks  present  data  on  successive  streamlines 
numbered  from  the  body  outward.  A  streamline  keeps  the  same  identification  number  from 
ray  to  ray.  Occasionally,  the  growing  boundary  layer  will  “swallow”  a  streamline,  and  the 
corresponding  identification  number  will  disappear.  From  time  to  time,  alternate  streamlines 
are  removed  from  the  grid  to  keep  the  number  of  points  within  bounds;  their  identification 
numbers  are  not  reused.  Two  adjacent  blocks  with  the  same  identification  and  the  same 
coordinates  represent  points  on  an  imbedded  shock.  An  “I.D.”  number  greater  than  1,000 
represents  an  auxiliary  shock  point  which  will  not  start  a  new  streamline.  Note  that  the  angle 
of  the  imbedded  shock  is  printed  as  the  last  item  on  the  previous  page,  Appendix  F. 

With  the  exception  of  the  Mach  number  and  the  equilibrium  indicator,  the  printed  data 
are  dimensional.  The  plasma  frequency  is  defined  as  tip  =  (56.4104  x  lO^Vn^  radians/sec. 
The  electrical  conductivity  of  the  plasma  is  stated  in  mho/m  and  is  calculated  under  the 
assumption  that  the  electron  temperature  is  equal  to  Tvib  (N2)  if  NVIB  >  0,  and  to  T 
otherwise. 

4.3.4  Edge  Conditions  File 

For  interfacing  with  BLIMP,  the  blunt  body  and  the  M.O.C.  parts  of  the  code  write  an 
unformatted  file  of  dimensional  boundary-layer  edge  data  to  Unit  3  (Table  1).  The  entries 
in  each  record  are,  in  order  arc,  radius,  pressure,  entropy,  and  concentrations,  where  arc 
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is  the  surface  distance  from  the  nose  vertex.  Both  the  surface  distance  and  the  body  radius 
are  stated  in  feet.  Species  concentrations  are  stated  in  mass  fractions  and  the  entropy  in 
calories/gm°K.  Pressure  is  in  dyne/cm2. 

4.4  DIAGNOSTICS 

Diagnostic  messages,  which  are  sometimes  helpful,  are  printed  by  CUNIFLOW  whenever 
the  computation  departs  from  the  usual  course. 

4.4.1  Non-Fatal  Conditions 

In  the  blunt  body  part  of  the  code,  an  integration  step  is  repeated  with  a  smaller  y-increment 
if  any  dependent  variable  goes  negative  or  if  the  relative  change  in  temperature  over  the  interval 
exceeds  the  input  constant  TTEST.  When  this  occurs,  either  of  the  following  messages  is 
printed: 

B2( — )  NEGATIVE  ON  RAY--,  RUNGE-KUTTA  STEP--,  Y  =  — ,  DELY - 

or 

T  TEST  FAILED  ON  RAY—,  RUNGE-KUTTA  STEP—,  Y  =  — ,  DELY  =  — 

The  elements  of  the  B2  array  are  defined  in  Table  1 1  below. 

Table  11.  Interpretation  of  the  B2  Array 

J  1  234567  .  36 

VARIABLE  ei  v  p  g  u  71  .  730 

The  information  in  this  output  is  sometimes  useful  in  understanding  the  chemical  processes 
taking  place  in  the  shock  layer  or  in  assessing  the  effects  of  a  tentative  change  in  the 
thermochemical  model.  As  the  step-size  will  be  cut  several  times  before  the  run  finally  aborts, 
a  study  of  the  dynamic  output  might  reveal  insidious  input  errors  such  as  misplaced  decimals. 
It  must  be  emphasized  that  step  cutting  is  a  normal  part  of  a  nonequilibrium  run. 

Another  nonfatal  diagnostic  is: 

ITERATION  IN  SUBROUTINE _ IS  NOT  CONVERGING. 

To  interpret  the  accompanying  data,  it  is  necessary  to  find  the  call  to  CNVRGE  in  the  named 
subroutine.  This  message  is  printed  whenever  the  local  iteration  counter  is  2:  NITER-3. 
CNVRGE  is  an  entry  in  subroutine  DIGNOS. 

A  message  which  may  or  may  not  indicate  a  physically  real  situation  is 
“SHOCK  LAYER  FULLY  VISCOUS” 
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which  means  that  the  calculated  mass  flux  through  the  boundary  layer  at  some  station  is 
greater  than  the  flux  into  the  shock.  When  this  condition  is  detected,  the  shock  shape  iteration 
is  abandoned  and  the  calculation  proceeds  as  if  the  shock  shape  were  converged.  The  user 
should  review  the  input  data  carefully,  paying  special  attention  to  ALT,  DENS,  and  TWALL, 
and  the  body  nose  radius.  An  RSHP  which  is  much  too  small  for  the  body  could  be  the 
culprit.  If  TWALL  is  much  too  low,  the  density  within  the  boundary  layer  will  be  too  high. 

Programmed  stops  are  produced  by  a  call  to  either  TRUBEL  or  PAUSE.  TRUBEL  prints 
the  location  of  the  call  before  stopping.  PAUSE  prints  the  location  of  the  call  and  then  calls 
POINT  to  print  the  data  at  all  locations  on  the  ray,  including  those  which  are  reserved  for 
special  storage  and,  consequently,  do  not  appear  in  the  usual  output.  The  first  two  blocks 
give  interpolated  data  at  the  bases  IR  and  IS  of  the  two  characteristics  (Fig.  4);  the  third 
block  contains  data  at  the  starting  point  of  the  streamline  integration.  This  will  differ  from 
the  point  on  the  previous  ray  only  if  the  integration  step  size  has  been  cut.  The  next  blocks 
present  data  at  the  three  successive  shock  points  and  the  three  successive  boundary-layer  or 
body  points.  The  following  blocks,  up  to  the  first  out-of-sequence  change  in  the  point  number, 
present  the  data  at  the  new  points.  The  remaining  blocks  present  the  data  on  the  previous 
ray.  When  the  trouble  occurs  while  calculating  the  flow  over  a  flare,  the  diagnostic  output 
becomes  somewhat  confusing  because  of  the  repeated  transfer  of  data  in  and  out  of  special 
storage.  The  key  to  interpretation  in  this  case  is  to  remember  that  data  on  the  downstream 
side  of  a  shock  are  always  stored  in  rows  4-6  of  the  data  arrays.  Information  on  the  upstream 
side  of  the  imbedded  shock  will  be  stored  in  rows  5-9. 

4.4.2  Fatal  Conditions 

Frequently  occurring  error  messages  are  discussed  below. 

The  following  messages  indicate  that  the  time  limit  has  been  exceeded: 

ERROR  AT  STATEMENT  1201  OF  ROUTINE  BREGN 

or 

ERROR  AT  STATEMENT  1601  OF  ROUTINE  BREGN 

'Timing  out”  in  BSHCK  on  the  first  iteration  means  that  something  is  terribly  wrong  with 
the  input  and/or  program.  Excessive  time  in  BREGN  may  be  attributable  to  bad  input;  check 
free-stream  density,  velocity,  temperature,  and  composition.  Excessive  time  may  also  be  caused 
by  free-stream  compositions  containing  small  amounts  of  some  substance  which  would  tend 
to  further  dissociate  at  the  high  temperatures  behind  the  shock.  For  example,  if  the  free  stream 
were  partially  dissociated  air,  the  step  sizes  needed  to  keep  the  O2  mole  fraction  non-negative 
behind  the  shock  may  be  very  small.  A  possible  cure  for  this  condition  would  be  to  increase 
GAMTST.  Sometimes  it  is  necessary  to  neglect  some  insignificant  free-stream  species.  If  the 
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input  is  correct  and  the  calculation  has  made  significant  progress,  use  the  restart  procedure 
and  possibly  change  TIME,  DELMIN. 

ERROR  AT  STATEMENT  1208  SUBROUTINE  LINK  2: 


“LINK  2”  is  the  former  name  of  BREGN.  The  significance  of  this  statement  is  that  the 
marching  procedure  could  not  get  started.  Check  the  input,  reduce  the  initial  step  size  by 
reducing  DELMIN,  and/or  increase  METHOD. 

ERROR  AT  STATEMENT  1425  OF  SUBROUTINE  LINK  2: 


Significance  of  this  statement  is  that  continuation  of  the  marching  procedure  requires  a  step-size 
less  than  DELMIN.  A  clue  to  the  cause  may  possibly  be  obtained  from  the  dynamic  dumps. 
If  the  input  is  satisfactory,  change  METHOD  and/or  DELMIN  and  restart  at  a  station  well 
upstream  of  the  point  where  the  catastrophic  instability  occurred. 

ERROR  AT  STATEMENT  30  OF  SUBROUTINE  BDYFCT: 

A  value  of  ZJUNC  <  0.0  was  input  for  first  body  segment,  or  PSIRAT  ~  1.0  was  input, 
or  RSHP  <  0  was  input. 

ERROR  AT  STATEMENT  90  OF  SUBROUTINE  INTR: 

The  significance  of  this  message  is  that  the  intersection  of  a  characteristic  with  the  previous 
ray  is  not  within  the  range  of  points  saved  for  interpolation.  The  coordinates  of  the  offending 
point  will  be  in  either  the  second  or  first  block  of  the  following  output.  The  third  block  is 
the  base  of  the  streamline  from  which  one  can  determine  the  “troubled”  process.  Most 
appearances  of  this  message  are  caused  by  a  discontinuity  at  the  junction  between  two  body 
equations.  If  the  message  appears  while  the  nose  is  still  being  calculated,  it  is  possible  that 
the  junction  Mach  number,  TESTM,  is  too  small,  that  the  computational  mesh  is  becoming 
too  skewed  because  the  flow  direction  is  steeply  inclined  with  respect  to  the  ray  direction, 
or  that  there  are  not  enough  points  on  the  starting  line.  To  fix  these  conditions,  rerun  the 
nose  portion  of  the  code  with  a  larger  TESTM  and  smaller  DELMAX.  Also  try  a  smaller 
value  for  ZETAMX.  If  the  error  occurs  well  back  on  the  afterbody,  it  is  probably  attributable 
to  a  divergent  iteration  process  somewhere  and  possibly  may  be  fixed  by  decreasing  the  step- 
size  (adjusting  ZETAMX,  CSST,  and/or  NMAX). 

ERROR  AT  STATEMENT  11  OF  SUBROUTINE  CFFCNT: 

The  significance  of  this  error  statement  is  that  the  Mach  number  <  1.0.  If  the  message  occurs 
early  in  the  M.O.C.  calculation,  it  is  probably  because  of  an  inadequate  starting  line  from 
the  nose  code.  The  statement  may,  of  course,  correspond  to  reality  if  the  flare  angle  is  too  steep. 
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4.5  OPERATING  ENVELOPE 

Figure  9  is  a  plot  of  the  CPU  time  required  to  execute  three  shock  iterations,  including 
the  boundary-layer  displacement  procedure.  The  shock  iteration  time  is  controlled  by  the 
number  of  integration  failures  with  the  accompanying  step  repetitions  and  step-size  cutting 
and,  therefore,  is  affected  somewhat  by  the  step-size  controls,  by  the  input  shock  parameters, 
the  integration  method,  and  the  thermochemical  model.  The  same  thermochemical  model 
(MODEL  =  0,  eight  species,  ten  reactions,  and  equilibrium  vibration)  and  control  settings 
(METHOD  =  10,  DELMIN  =  10-*,  and  GAMTST  =  10~7)  were  used  for  all  the 
calculations.  The  shock  damping  factor,  SKAP,  was  increased  from  0.35  to  0.50  for  the  high- 
velocity  cases  to  secure  better  convergence.  The  procedure  for  constructing  the  graph  was 
to  start  at  15,000  ft/sec  and  build  outwards  in  both  directions  using  as  initial  shock  parameters 
for  each  case  the  final  (converged)  values  from  the  previous  calculation. 


Figure  9.  Effect  of  free-stream  conditions  on  CPU-time  for  nonequilibrium 
blunt  body  solutions. 

The  increase  in  CPU  time  at  low  speed  was'  caused  by  the  presence  of  a  minute  amount 
of  nitrogen  atom  (y  =  1.8  x  10-9)  which  was  being  consumed  by  the  reaction,  N  +  O 
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-  NO+  +  e,  faster  than  it  was  being  produced  and  could  have  been  forestalled  by  choosing 
a  larger  value  for  GAMTST,  such  as  10  _6. 

The  increase  at  high  speed  arose  from  the  virtually  complete  dissociation  of  O2  near  the 
stagnation  point,  where  the  flow  is  nearly  in  equilibrium  (AVER.CHI  =  3.3  x  10-5).  In 
contrast,  at  the  lowest  speed,  the  flow  is  far  from  equilibrium  (AVER.CHI  =  0.98)  at  the 
boundary-layer  edge.  This  increase  illustrates  the  difficulty  in  using  a  nonequilibrium  code 
in  near  equilibrium  conditions. 

The  curves  of  Fig.  9  are  labeled  with  the  parameter,  p^Rn-  It  has  been  shown  (Ref. 
33)  that  the  product  of  the  free-stream  density  with  a  characteristic  body  length  constitutes 
a  simple  and  exact  similitude  for  nonequilibrium  processes  which  involve  two-body  collisions 
only.  Since  the  full  chemistry  of  high-temperature  air  involves  the  three-body  recombinations 
of  oxygen,  nitrogen,  and  nitric  oxide,  binary  scaling  cannot  apply  exactly  to  air,  especially 
at  lower  altitudes  where  the  flow  is  close  to  equilibrium,  or  over  the  afterbody  where  the 
flow  is  recombining. 

Figure  10  is  a  plot  of  the  scaling  factor  against  velocity  for  a  number  of  nonequilibrium 
solutions  which  have  been  obtained  over  the  years  using  various  versions  of  CUNIFLOW 
and  various  reaction  sets.  It  is  probably  safe  to  assume  that  CUNIFLOW  will  run  for  free- 
stream  conditions  which  can  be  plotted  on  the  figure  below  the  indicated  upper  boundary. 
Near  the  upper  boundary  it  will  probably  be  necessary  to  include  both  radiation  and  thermal 
nonequilibrium  options.  For  very  blunt  bodies  or  bodies  with  large  cone  angles,  the  base 
radius  might  be  a  more  appropriate  characteristic  length  than  the  radius  of  curvature  of  the 
nose. 


5.0  CONCLUDING  REMARKS 

The  structure  and  use  of  the  CUNIFLOW  code  has  been  discussed.  CUNIFLOW  is  a 
very  fast,  robust,  reacting  flow  code  which  has  become  the  standard  method  at  AEDC  for 
computing  the  nonequilibrium  flow  about  axisymmetric  or  two-dimensional  bodies.  By 
combining  CUNIFLOW  with  the  well-known  BLIMP  boundary-layer  code,  it  has  been  possible 
to  obtain  nonequilibrium  solutions  at  low-altitude,  high-speed  conditions  where  other  external 
flow  codes  will  not  converge.  Furthermore,  it  is  possible  to  make  computations  for  flow  media 
other  than  air  and  for  nonuniform  upstream  conditions.  The  economy  and  robustness  of 
the  code  are  illustrated  by  plots  (Section  4.5)  of  CPU  time  as  a  function  of  velocity  and  binary 
scaling  parameters. 
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Figure  10.  Nonequilibrium  air  solutions  on  missile-like  bodies. 

The  available  thermochemical  models,  the  flow  and  geometry  options,  and  the  overall 
code  structure  are  presented  in  Sections  2,0  and  3.0.  The  procedures  for  determining  the 
location  of  the  bow  shock,  as  well  as  an  imbedded  shock  and  for  approximating  the 
displacement  effect  of  the  boundary  layer  have  been  explained  in  Section  3.0. 

Instructions  for  using  CUNIFLOW  are  presented  in  Section  4.0  and  in  Appendixes  B, 
E,  and  F.  Appendixes  A  and  D  present  mathematical  details  and  the  code  structure.  Differences 
between  the  present  and  previous  versions  of  CUNIFLOW  are  listed  in  Appendix  C. 
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APPENDIX  A 

MATHEMATICAL  DETAILS 
1.0  BLUNT  BODY  EQUATIONS 

The  nondimensional  conservation  equations  in  the  curvilinear  coordinate  system  of  the 
blunt  body  code  are 


Mass: 


a-i  +  l£ - S -  1  —  (pu)  +  -  — 

dy  v  dy  v(l  +  xy )  g  dx  r  3x 
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Vibrational  Energy: 
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In  these  variables,  the  partial  pressure  of  the  i*  species  may  be  written 


Pi  = 


(A-8) 


and  Rad  is  the  radiant  energy  lost  by  the  gas. 

The  nondimensional  quantities  used  above  are  defined  according  to  the  following  system: 
Types  of  Quantity  Reference 


Length 

Velocity 

Density 

Pressure 


e» 
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Types  of  Quantity 

Reference 

Energy  of  mixture 

V*. 

Energy  of  particular  species 

Molecular  weight 

5. 

Concentration  (yj) 

NS 

£  D.  = 

2.0  DERIVATIVE  ALGORITHM 

Derivatives  in  the  x-direction  (parallel  to  the  bow  shock)  are  evaluated  in  subroutine  PART 
using  the  formula 


dB|  =  SBi-3  -  12B)_;  -  39Bj_i  +  39Bi+L  +  12Bi+2  -  5Bi+3 
dx  96Ax 

where  B  represents  any  variable  and  i,  its  position  in  the  data  array. 

Close  to  the  ends  of  the  array,  special  formulae  are  necessary: 

For  i  =  1, 


dB|  _  -  133B|  +  84B2  +  207B3  -  192B4  -  39B5  +  108B6  -  3SB7 
dx  96  Ax 


(A- 10) 


For  i  =  2, 

dB2  -35B]  -  28Bj  +  49B3  +  32B4  -  17B;  -  4B6  +  3B? 
dx  96 Ax 


(A-ll) 


For i  =  3, 

dB3  3B3  -  44B;  -  25B3  +  64B4  +  17B5  -  20B6  +  5B7 
dx  96 Ax 


(A- 12) 
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For  i  =  N— 2  or  N  - 1,  where  N  is  the  number  of  elements  in  the  array, 

—  _  **‘+1  ~  Bj-i 
dx  2  Ax 


(A- 13) 


For  i  =  N, 


^Bn  _  Bn —2  -  4Bn_i  +  3Bn 
dx  2  Ax 


(A-14) 


3.0  INTEGRATION  ALGORITHM 

The  explicit  algorithm  devised  by  Lomax  for  stiff  equations  is 

B?t,  -  B,  +  h«l,M)B„ 

Bf>,  =  B„  +  h«k,M)B*-,1) 


where 
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dB 

s  ds 


S-Sn 


h  =  Sji  Sn  _  | 

and  B„  represents  any  variable  evaluated  at  point  s„. 
4.0  EQUATIONS  FOR  THE  M.O.C.  CODE 


(A- 15) 
(A-16) 

(A-17) 

(A*  IS) 
(A- 19) 


The  afterbody  program  uses  the  “reference-plane”  method  of  characteristics.  The 
nondimensional  compatibility  and  characteristic  equations  used  are  listed  below: 


CA  dp  ±  d<£  =  CBdJ- 

(A-20) 

CA  =  (cot  p)/ gu2 

(A-21) 

<t>  sin  ^  sin  0 

—  +  F  sin  /3 

r 

(A-22) 
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which  must  be  satisfied  along  the  curves 

£  =  tan (£  ±  ,8) 

In  addition,  the  thermodynamic  and  reaction  equations 
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are  applicable  along  the  streamlines  specified  by 

dr 


dz 


=  tan  <f> 


The  density  is  calculated  from  the  equation  of  state 


(A-29) 

(A-30) 


(A-31) 


(A-32) 
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APPENDIX  B 

INPUTTING  THERMOCHEMICAL  PROPERTIES 

If  a  problem  involves  chemical  or  vibrational  nonequilibrium  (MODEL  =  0,3,4), 
subroutine  CHEMIN  is  called  to  read  a  properties  file  from  Unit  8.  The  format  and  contents 
of  this  file  will  depend  on  the  model,  but,  in  general,  reaction  data  are  input  first  —  two 
logical  records  for  each  reaction  —  followed  by  thermodynamic  data  for  each  species.  All 
input  is  formatted,  and  each  block  of  data  is  terminated  by  a  record  with  the  characters  “END” 
in  the  first  three  columns.  The  code  will  count  the  entries  in  each  category  and  will  select 
the  data  for  the  species  specified  in  Namelist  UPSTREAM  (Section  4.2.2)  or  named  by  the 
reaction  records  from  the  data  base  on  Unit  8.  At  present,  dimension  statements  allow  72 
reactions  involving  not  more  than  30  species. 

A  diagram  of  the  input  procedure  is  shown  in  Fig.  B-l .  The  various  record  types  are 
identified  by  the  numbers  in  the  upper  right  corners  of  the  record  symbol.  The  input  file 
“NEWCHEM”  which  was  used  for  the  exploration  of  the  operating  envelope  in  Section  4.5 
is  shown  in  Fig.  B-2. 

Table  B-2  from  Ref.  8  lists  a  number  of  reactions  among  the  components  of  heated  air. 

Other  tabulations  will  be  found  in  Refs.  15,  30,  and  34-36. 
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Figure  B-l.  Thermodata  input  schematic. 
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02  26.0  I . 6 1SB  €-03  1.01444  E+02  3.333333E-01  2256.34  E+OO 

N2  33. D  1.1153  E-05  0.5  1.5400  +02  3.333333E-01  3374.03  E+OO 

□  1  O  1  O  O  0  1.0 

2.46855  E+12  2.1032  E+OO  0.0  E+OO 

1 

0.0  E+OO  7 .548  E-22  1.077  E-21  7.532  E-20  1.604  E-19 

5.0  E+OO  3.0  E+OO  1.0  E+OO  5.0  E+OO  1.0  E+OO 

N  1  O  1  0  0  N  1.0 

4.712342  E+12  1.67997  E+OO  0.0  E+OO 

i 

0.0  E+OO  9.127  E-20  9.152  E-20  1.369  E-19  4.1819E-19 

4.0  E+OO  6.0  E+OO  4.0  E+OO  6.0  E+OO  1.2  E+01 

E-  1  0  0  O  1  .  E-1.0 

0.0  E+OO  -1.423517  E+01  0.0  E+OO 

1  1 

2.0  €+00 

AR  I  O  1  O  O  AH  I. 

0.0  E+OO  1.86557  E+OO  0.0  E+OO 

0.0  .  E+OO  4.4215E- 19  4.4503E-19 
1.0  E+OO  5.0  E+OO  3.0  E+OO 
02  2  1  0  0  0  2.0 

0.0  E+OO  1.21618  E+OO  2.23897  E+03 

I 

0.0  E+OO  3. 75B3E-20  6.26B5E-20  1.713  E-19  2.3641E-19 

3.0  E+OO  2.0  E+OO  1.0  E+OO  3.0  E+OO  3.0  E+OO 

N2  2  1  0  O  0  N  2.0 

0.0  E+OO  -4.21630  E-OI  3.35324  £+03 

1  1 

0.0  E+OO  2 . 3866E -  19  2.830  E-19  2.84B  E-19 

1.0  E+OO  3.0  E+OO  6.0  E+OO  1.0  E+OO 

NO  2  1  O  O  O  N  1.  0  1.0 

8.98898  E+11  1.2325&  E+OO  2.69918  E+03 

1  1 

0.0  E+OO  5.744  E-22  2.087  E-19  2.180-E19 

2.0  E+OO  2.0  E+OO  2.0  E+OO  4.  E+OO  0.0  E+OO 


Figure  B-2.  Thennochcmlcal  input  file  “NEWCHEM”. 
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1 -O  E+OO  6.0  E+OO  3 

END 
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Figure  B-2.  Concluded 
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Table  B-l.  Formats  for  Thermochemical  Input 


Record  Type  I 


Format  <6A4,3F3.0,I3,2F3.0) 


(SPECIK(1,K,K  =  1,3) 
(SPECIK(J,K),K  =  4,6) 
CZI,CWI,CDI 

KFI1ND 

(A(I,K),K=  1,2) 


Reactant  names  left  justified 

Product  names  left  justified  beginning  in  Column  13. 
Reaction  type  [Eq.  (22),  Sec  2.3.2]. 

Set  just  one  to  1  and  the  others  =  0. 

Set  =  1  if  reaction  rate  constants  apply  to  reverse  of  the 
reaction  as  written. 

Set  =  1  if  this  reaction  is  coupled  to  the  kth  nonequilibrium 
vibrator.  [Eq.  (20)] 


Record  Type  2 


Format  (5E15.7) 


AK 


Constant  Ak  of  rate  coefficient  [Eq.  (21)](cm3/mole-sec  or 
cm6/mole-sec) 


BK 

CK 

DK 

Record  Type  3 

osc 

o 


Exponent  of  T  in  Eq.  (21),  Bkj 
Activation  temperature  [Eq.  (21)],  °K 
Pki  in  Eq.  (21) 

Format  (A3,  F4.0,  5E14.7) 
Vibrator  name  —  left  adjusted 
Number  of  vibrational  levels  for  this  vibrator 


TAUA 

TAUB 

TAUC 

TAUD 

THETA 


Coefficient  in  Eq.  (27)  (gm/cm/sec) 

Exponent  of  T0,  Eq.  (27) 

Argument  of  exponential,  Eq.  (27) 
Temperature  exponent  in  argument  of  Eq.  (27) 
Characteristic  vibrational  temperature  (°K) 
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Table  B-l.  Continued 

Record  Type  4  Format  (7E11.6) 

EPSCNT  Reserved  for  coefficients  of  possible  vibrational  energy  Fit. 

Not  used  at  present. 

Record  Type  5  Format  (A4,4I3,8F3.0) 


SUB 

Name  of  species  —  left  adjusted 

NUMA 

Number  of  atoms  in  species 

LU(1) 

1  if  species  is  neutral  molecule 

0  otherwise 

LU<2) 

1  if  species  is  neutral  atom 

LU(3) 

1  if  species  is  ionized  atom 

LU(4) 

1  if  species  is  electron 

EL(1) 

Name  of  element 

B(I) 

Number  of  atoms  of  EL(1)  in  molecule 

EL(2) 

Name  of  second  element  in  molecule 

B(2) 

Number  of  atoms  of  EL(2)  in  molecule 

EL(3) 

Name  of  third  element  if  present 

B(3) 

Number  of  atoms  of  EL(3)  in  molecule 

EL(4) 

Name  of  fourth  element,  if  present 

B{4) 

Number  of  atoms  of  EL(4)  in  molecule 
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Record  Type  6 

HJOP 

BJ 

THETV 

Record  Type  7 
(LIT(I,J),J  =  1,  NRAD) 

Record  Type  8 
HCFNT 

Record  Type  9 
FCFNT 

Record  Type  10 

(SUB(I),I  =  1.4) 

EL(1) 

B(l) 

EL(2) 

B(2) 

EL(3) 

B(3) 


Table  B-l.  Continued 

Format  (3E15.7) 

Heat  of  formation  at  1  atm,  ergs/mole 
Constant  b  in  Eq.  (4),  erg/ mole 
Characteristic  vibrational  temperature  °K 

Format  (1713) 

Matrix  elements  connecting  this  species  to  Jth  radiating 
system 

Format  (7E11.6) 

Array  of  electron  energy  levels  for  this  species  — 
cal/molecule 

Format  (7E11.6) 

Degeneracies  of  electron-energy  levels 

Format  (4A4,8X14(A2,F3.0),A1,2F10.3,F13.5) 
Species  name  —  up  to  16  characters 
Name  of  element 

Number  of  atoms  of  EL(1)  in  molecule 
Name  of  second  dement  in  molecule 
Number  of  atoms  of  EL(2)  in  molecule 
Name  of  third  element,  if  present 
Number  of  atoms  of  EL(3)  in  molecule 
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Table  B-l.  Continued 

EL(4) 

Name  of  fourth  element,  if  present 

B(4) 

Number  of  atoms  of  EL(4)  in  molecule 

PH 

Enter  “G” 

T1 

Boundary  between  lower  and  middle  temperature  ranges 

T2 

Boundary  between  middle  and  upper  temperature  ranges* 

MW 

Molecular  weight 

’Note:  The  thermal  properties  tables  from  Lewis  Research  Center  have  only  two  ranges: 
100  -  1000  °K  and  1000  -  5000  °K. 

Record  Type  11 

Format  (5E15.8) 

(ARRAY(N,1),N  =  1,5) 

First  5  coefficients  of  mid-range  polynomial  fit,  Eqs.  (7)-(9) 

(ARRAY(N,1),N  =  6,7), 
(ARRAY(N,2),N=  1,3) 

Rest  of  mid-range  and  first  3  low-range  coefficients 

Record  Type  12 

Format  (5E15.8) 

(ARRAY(N,2),N  =  4,7) 

Remainder  of  low-range  thermofit  coefficients 

HF298 

Heat  of  formation  at  298.15  °K,  cal/mole 

Record  Type  13 

Format  <7(F3.0,I3),30X,2A4) 

(B(N),JAT(N),N  =  1,7) 

B(N)  number  of  atoms  of  atomic  number  JAT  in  this 
molecule 

SUB(l),  SUB(2) 

Species  name  —  up  to  8  characters 
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Table  B-l.  Concluded 


Record  Type  14 


Format  (6E9.6,2E6.0,I1) 


HF298 


Heat  of  formation  at  298.15  °K,  cal/mole 


(ARRAY(J,1),J  =  1,5)  Thermofit  coefficients  in  Eq.  (12);  a2  =  ARRAY  (1,1) 

Lower  temperature  range,  cal/mole  or  cal  mole  °K 


T1  Upper  limit  of  lower  temperature  range 

T2  Reserved  for  upper  limit  of  middle  temperature  range.  Not 

currently  implemented. 


KPH 


Enter  1  if  this  species  is  a  gas 
Enter  0  if  otherwise 


Record  Type  15 


HF298 

(ARRAY(J,2),J  =  1,5) 
EMI 


Format  (6E9.6,2E6.0,I1) 

Heat  of  formation  at  298.15  °K.  This  must  be  the  same  as 
on  previous  record. 

Thermofit  coefficients  for  upper  range* 

Molecular  weight 


*At  present.  If  a  third  temperature  range  is  necessary,  alter  subroutine  CHEMIN  to  read 
a  third  record  in  record  type  15  format;  redimension  EPSCNT  for  the  number  of  species 
involved  (30  at  present),  and  enter  appropriate  value  for  T2.  Use  EPSCNT  to  carry  the 
coefficients  to  subroutine  THERMO,  which  will  also  need  modification. 


70 


AEDC-TR-8S-31 


Table  B-2.  Reactions  and  Rates* 


Refection 

Type 

Direction 

Ak 

O 

Aq 

Units 

02  +  O  =  20  +  O 

W 

FWD 

2.1  X  10“ 

-0.5 

59380 

1 

cmVmale-sec 

Oj  +  O2  =  20  +  Oj 

W 

FWD 

3.56  X  IO21 

-1.5 

59380 

1 

cmVmole-sec 

O2  +  M.  =  20  +  M 

z 

FWD 

1.190  x  1021 

-1.5 

59380 

1 

cmVmole-SK 

N2  4  N  =  2N  4  N 

w 

REV 

7.50  X  10“ 

-1.5 

0 

craVmole2-sec 

N2  +  N2  =*  2N  4  N2 

w 

REV 

1.50  X  IO20 

-1.5 

0 

— 

1 

cm6/mole2-sec 

N2  +  M  -  2N  +  M 

z 

REV 

5.00  x  JO19 

-1.5 

0 

— 

1 

cmVmole-sec 

N0  +  M  =  N  +  0  +  M 

D 

FWD 

5.18  X  IO11 

-1.5 

75490 

... 

~ 

cmVmolc-sec 

N  +  Oj  =  NO  +  0 

w 

FWD 

1.0  X  1012 

4  0.5 

3120 

— 

— 

cm3 /mole-scc 

O  4  N2  =  NO  4  N 

w 

FWD 

5.0  x  10“ 

0 

38016 

— 

— 

cm2 /mole-sec 

N  +  O  -  NO+  +  e“ 

w 

REV 

1.8  x  IO2' 

-1.5 

0 

— 

— 

cm2 /mole-sec 

20  =  e“  4  02* 

w 

REV 

1.9  x  IO21 

-1.5 

0 

— 

cmVmole-sec 

2N  =  e~  +  Na+ 

w 

REV 

9.6  x  IO21 

-1.5 

0 

— 

— 

cm2 /mole-sec 

O  4  N  =  N  4  e_  4  0+ 

w 

REV 

6.0  x  1024 

-2.5 

— 

— 

— 

cm4/ mole-sec 

20  =0  +  0+  +  e“ 

w 

REV 

6.0  x  IO24 

-2.5 

... 

_ 

— 

cmVmole-sec 

O  +  N2  =  N2  +  e"  +  0+ 

w 

REV 

2.0  x  IO24 

-2.5 

— 

_ 

— 

cmVmole-sec 

O  +  02  =  0+  +  e"  +  02 

w 

REV 

9.0  x  IO24 

-2.5 

— 

— 

— 

cm6/ mole-sec 

O  4  NO  =  NO  4  e~  4  0* 

w 

REV 

1.0  X  102* 

-2.5 

... 

— 

— 

cm6/mole-sec 

2N  =  N  4  e~  4  N+ 

w 

REV 

6.0  x  IO24 

-2.5 

... 

— 

— 

cm4/rnole-sec 

N  4  O  =  O  4  e~  4  N+ 

w 

REV 

6.0  X  1024 

-2.5 

— 

cm6/mole-scc 

N  4  N2  =  N2  4  e~  4  N* 

w 

REV 

2.0  X  IO26 

-2.5 

— 

— 

-- 

craVraole-sec 

N  +  02  =  02  +  e_  +  N+ 

w 

REV 

9.0  X  IO24 

-2.5 

— 

— 

— 

cmVmole-sec 

N  4  NO  =  NO  4  e"  4  N  + 

w 

REV 

1.0  X  IO28 

-2.5 

— 

— 

— 

cmVmole-sec 

N2  +  0+  =  0  +  N2+ 

w 

REV 

7.8  x  IO11 

0.5 

— 

— 

— 

cmVmole-sec 

N  +  NO+  =  N2  +  0+ 

w 

REV 

7.8  x  10“ 

0.5 

— 

— 

— 

cmVmole-sec 

0  4  0^  =  02  4  0* 

w 

REV 

7.8  x  10" 

0.5 

cmVmole-sec 

O  +  NO+  =  NO  4  0+ 

w 

REV 

7.8  x  10" 

0.5 

— 

— 

— 

cmVmole-sec 

N  4  =  NO  +  O* 

w 

REV 

7.8  x  10" 

0.5 

— 

... 

— 

cmV  mole-sec 

N2  4  N+  =  N  4  Nj 

w 

REV 

7.8  x  10" 

0.5 

— 

... 

~ 

cmVmole-sec 

M+0j*=02  +  N+ 

w 

REV 

7.8  x  10" 

0.5 

— 

— 

cmVmole-sec 

0  4  N0+  =  02  4  N+ 

w 

REV 

7.8  x  10" 

0.5 

— 

— 

~ 

cmVmale-sec 

N  4  N0+  =  NO  4  N+ 

w 

REV 

7.8  x  10" 

0.5 

— 

— 

— 

cmVmole-sec 

0  4  N2f  =  NO  4  N+ 

w 

REV 

7.8  x  IO11 

0.5 

— 

— 

... 

cmVmole-sec 

N  4  0+  =  0  4  N  + 

w 

REV 

7.8  x  IO11 

0.5 

— 

... 

— 

cmVmole-sec 

Nj  4  Oj  =  2  NO 

w 

FWD 

9.1  X  IO24 

-2.5 

cmVmole-sec 

Ar+  4  2e"  =  Ar  +  e' 

w 

— 

2.064  x  I0,T 

0.663 

148000 

— 

— 

cmVmole-sec 

*  Rales  are  of  ihe  form  k  =  Ak  TBk  exp  (-Ci/T0*) 
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APPENDIX  C 
CHANGES  SINCE  1967 


Appendix  C  lists  those  features  of  CUNIFLOW  which  have  been  added,  deleted,  or 
substantially  changed  since  the  publication  of  Ref.  10. 

Magneto-gas-dynamics  —  All  references  to  magneto-gas-dynamic  properties  and  related 
calculations  have  been  removed  from  the  blunt-body  case. 

Equilibrium  wall  temperature  —  Reference  10  describes  an  elaborate  procedure  for  balancing 
the  radiant  and  convective  heat  transfer  to  the  body  surface  with  the  radiant  energy  emitted 
from  the  surface  to  determine  the  wall  temperature.  Since  the  wall  temperature  option  was 
inaccurate  and  slow,  it  has  been  removed. 

Integration  Algorithm  —  The  Lomax  algorithm  described  in  the  text  has  replaced  the  Runge- 
Kutta  and  Treanor  algorithms  discussed  in  Ref.  10. 

Nonuniform  Free  Stream  —  Inputting  free-stream  data  in  a  tabular  format  and  using  linear 
interpolation  for  pre-shock  properties  has  replaced  the  cumbersome  polynomial  fits  described 
in  Ref.  10. 

Centerline  —  The  original  version  of  CUNIFLOW  made  no  provision  for  including  the 
centerline  in  the  calculation. 

Comers  —  The  original  version  of  CUNIFLOW  (Ref.  10)  made  no  provision  for  sharp  corners. 
The  facility  for  treating  embedded  shocks  and  Prandtl-Meyer  corners  was  soon  added  and 
was  described  in  Ref.  1 1 .  The  current  implementation  differs  significantly  from  that  of  the 
early  version. 

Input-Output  —  The  current  input  procedure  is  entirely  different  from  that  of  Refs.  10  and 
11.  The  current  printout  is  organized  in  roughly  the  same  way  as  that  of  the  original,  but 
a  number  of  quantities  which  were  of  interest  only  for  code  development  have  been  omitted 
and  the  printout  has  been  made  dimensional.  The  binary  plotting  file.  Unit  7,  and  the  BLIMP 
communication  file,  Unit  3,  have  been  added. 

Equilibrium  —  Reference  1 1  described  an  elaborate  equilibrium  chemistry  package  which 
was  supposed  to  furnish  the  equilibrium  composition  as  well  as  the  fluid  dynamic  quantities. 
Because  the  composition  package  was  slow  and  temperamental,  it  was  soon  abandoned;  only 
the  RGAS  and  SERCH  routines  have  been  retained. 
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APPENDIX  D 
CODE  STRUCTURE 

The  overall  structure  of  CUNIFLOW  is  displayed  schematically  in  Figs.  D-l  through  D-4. 
A  matrix  of  subroutine  references  is  shown  in  Table  D-l.  SMALLH  is  actually  an  ENTRY 
point  in  subroutine  THERMO;  CNVRGE  is  an  ENTRY  in  DIGNOS;  FELDA  is  an  ENTRY 
in  subroutine  FELD;  STRIP  is  an  ENTRY  in  ISTRIP;  and  QUELLE  is  an  ENTRY  point 
in  subroutine  SOURCE.  A  directory  of  unit  operations  is  presented  in  Table  D-2. 


Figure  D-2.  Transition  from  elliptic  to  hyperbolic  region. 
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Figure  D-3.  Elliptic  region  calculation. 
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Figure  D-4.  Hyperbolic  region  calculation. 


75 


Table  D-l.  Subroutine  Cross-Reference  Chart 
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Key  to  Subroutines 


1. 

AMBEN 

36. 

LEAST 

2. 

BACK 

37. 

MOCSOL 

3. 

BBIN 

38. 

OUT 

4. 

BDYFCT 

39. 

OUTPUT 

5. 

BLAYER 

40. 

PART 

6. 

BLAYST 

41. 

PAUSE 

7. 

BODYB 

42. 

POINT 

8. 

BPTS 

43. 

POLATE 

9. 

BREGN 

44. 

PREPAR 

10. 

BRIDGE 

45. 

PRNDTL 

11. 

BSHCK 

46. 

RADAB 

12. 

CFFCNT 

47. 

RADION 

13. 

CHEMIN 

48. 

RAY 

14. 

CHEMOUT 

49. 

RAYB 

15. 

CNVRG 

50. 

RAYC 

16. 

COUPLE 

51. 

REACTB 

17. 

DAUX 

52. 

RGAS 

18. 

DERAB 

53. 

RNGKT1 

19. 

DERBS 

54. 

SERCH 

20. 

DIGNOS 

55. 

SHIF 

21. 

DIVIDE 

56. 

SHIFTIT 

22. 

ELDENT 

57. 

SHMAIN 

23. 

EQILIB 

58. 

SHOCK 

24. 

EQUIB2 

59. 

SIMEQ 

25. 

ERRMSG 

60. 

SIMSOL 

26. 

FELD 

61. 

SMALLH 

27. 

FLARE 

62. 

SOURCE 

28. 

FLO 

63. 

START 

29. 

FLOW 

64. 

STIP 

30. 

FRMAIN 

65. 

STRMB 

31. 

HATM 

66. 

THERMO 

32. 

IMBED 

67. 

TRANS 

33. 

INTGRL 

68. 

TRUBEL 

34. 

INTR 

69. 

WRITE2 

35. 

ISTRIP 
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Table  D-2.  Directory  of  Unit  Operations 


Operation 


Input  free-stream  conditions 
Input  options  and  controls 
Input  body  shape  equations 
Input  initial  shock  shape  parameters 

Free-stream  conditions  at  z,r 

Evaluate  shock  equation  for  blunt  body  shock 

coordinates 

Iterate  for  M.O.C.  shock  points 

Post-shock  flow  —  elliptic  region 
Post-shock  flow  conditions  —  hyperbolic  region 

Calculate  pressure  at  trial  shock  location  from 
interior  flow-field  data 

Copy  dependent  variable  data  from  one  array  to 
another  —  elliptic  region 
Copy  data  from  one  cell  to  another 
—  hyperbolic  region 

Evaluate  x-derivatives  —  elliptic 

Evaluate  y-derivatives  —  elliptic 

Evaluate  species  source  terms 
Evaluate  vibrational  relaxation  rates 

Evaluate  radiant  energy  loss  rate  —  elliptic 
Evaluate  radiant  energy  loss  rate  —  M.O.C. 

Input  nonequilibrium  thermo  properties 
Input  equilibrium  thermodynamic  tables 
Evaluate  nonequilibrium  thermo  properties 
Evaluate  equilibrium  thermo  properties 


Location 

BBIN 

BBIN 

BBIN 

BBIN 

AMBEN 

BSHCK 

SHOCK 

FLOW 

FLO 

BACK 


SHIFTIT 

SHIF 


PART 

DERBS 

SOURCE 

SOURCE 

RADION 

RADAB 

CHEMIN 

RGAS 

SMALLH 

RGAS 
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Operation 

Location 

Arguments  for  RGAS  —  elliptic 

Arguments  for  RGAS  —  hyperbolic 

EQUIB2 

EQILIB 

Integrate  from  shock  to  body  —  elliptic 

BREGN 

Intersect  ray  with  body  —  elliptic 

Intersect  ray  with  body  —  hyperbolic 

BPTS 

RAY 

Evaluate  polynomials  for  radius  and  slope 

Input  body  polynomials 

BDYFCT 

BBIN 

Evaluate  viscosity  at  stagnation  point 

Look  up  boundary-layer  integrals 

Calculate  mass  flux  and  velocity  thickness 
—  elliptic  region 

TRANS 

INTGRL 

BLAYER 

Iterate  for  point  at  edge  of  boundary-layer 
—  hyperbolic  region 

COUPLE 

Adjust  shock-point  coordinates  to  match  mass 
flux  —  elliptic  region 

BREGN 

New  shock  shape  parameters  by  least-squares  fit 

LEAST 

Select  points  for  hyperbolic  starting  line  and 
write  to  Unit  2 

BRIDGE 

Adjust  body  shape  for  best  match  at  initial  line 

BLAYST 

Extrapolate  elliptic  solution  to  body  or 
boundary-layer  edge 

BRIDGE 

Print  elliptic  region  solution 

OUT 

Read  M.O.C.  starting  line  from  Unit  2 

Read  M.O.C.  starting  line  from  cards 

START 

START 

Manage  M.O.C.  calculation  of  hyperbolic  region 

FRMAIN 
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Operation  Location 

Locate  next  ray  RAY 

M.O.C.  calculation  at  point  on  body  surface  BODYB 

Manage  calculation  of  Prandtl-Meyer  expansion  RAYB 

Calculate  flow  conditions  behind  Prandtl-Meyer  PRNDTL 

expansion  fan 

M.O.C.  calculation  at  point  in  interior  of  shock  FELD 

layer 

Perform  one  cycle  of  interior  point  M.O.C.  MOCSOL 

calculation 

Integrate  for  frozen  temperature  and  velocity  STRMB 

along  streamline  segment 

Integrate  species  conservation  equations  along  REACTB 

streamline  segment 

Derivatives  of  composition  and  vibrational  DERAB 

energy  in  streamwise  direction 

Interpolate  for  dependent  variables  INTR 

Test  for  completion  (time  or  position)  STIP 

Write  restart  file  (Unit  10)  STIP 

Rearrange  storage  omitting  alternate  points  DIVIDE 

Print  M.O.C.  solution  on  ray  OUTPUT 

Write  communication  file  for  BLIMP  OUTPUT 

Calculate  collision  and  plasma  frequencies  ELDENT 

M.O.C.  solution  at  constant  pressure  boundary  BODYB 
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Operation 

Location 

Manage  M.O.C.  calculation  with  embedded 
shock 

SHMAIN 

Calculate  initial  angle  for  embedded  shock 

FLARE 

Locate  embedded  shock  and  calculate  flow 

IMBED 

Print  diagnostic  message  and  stop  —  elliptic 

TRUBEL 

Print  diagnostic  message  and  stop  —  MOC 

PAUSE 

Print  diagnostic  and  continue 

ERRMSG 

Print  diagnostic  message  and  continue 

CNVRGE 

Read  restart  file  (Unit  1),  free-stream  conditions, 
thermochemistry,  body  shape 

Read  restart  file  and  restart  elliptic  region 
calculation 

BBIN 

BREGN 

Read  restart  file  and  locate  desired  restart 

station 

START 
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The  nonequilibrium  shock  shape  parameters  RSHP,  SIP,  and  SI  which  have  been  obtained 
for  various  free-stream  conditions  and  thermochemical  options  are  tabulated  in  Table  E-l . 
Unless  otherwise  indicated,  the  reaction  set  of  Ref.  7  was  used  in  making  the  calculations. 


Table  E-l.  Nonequilibrium  Shock  Shape  Parameters 


UM,  Kfps 

eroRN,  gm/cm2 

RSHP/Rn 

SIP 

SI  x  100 

RMS  Dev 
x  I04 

Comment 

2.76  x 

10-7 

mm 

0.08412 

_ 

3.80 

No  B.L.,.  vib. 

1.95  x 

I0-6 

:IS31l 

0.1473 

1.283 

No  B.L.,  vib. 

23.1 

4.119  x 

10-6 

0.15163 

0.739 

7.455 

B.L.,  no  vib 

14.7 

7.11  x 

10-6 

.  1.2566 

0.0931 

1.372 

4.687 

Vib.,  no  B.L. 

10' 

2.0334 

x  10~5 

1.1967 

0.10437 

1.197 

2.119 

NV1B  -  0,  w/  B.L. 

12 

1.1697 

0.11161 

1.064 

5.101 

15 

1.1154 

0.11879 

1.360 

17.1 

1.1173 

0.13880 

0.927 

IS 

1.1186 

0.14699 

— 

5.61 

19 

1.1199 

0.15474 

7.47 

1.1192 

0.15872 

8.873 

21 

1.1177 

0.16137 

0.877 

10.08 

22 

1.1127 

0.16130 

0.836 

9.18 

24 

1.1023 

0.16430 

0.7B1 

5.68 

26 

1.0976 

0.16614 

0.729 

7.673 

1.0843 

0.16117 

0.656 

4.765 

2.0334  X  10-5 

1.0708 

0.15227 

0.579 

NVIB  =  0,  w/  B.L. 

3.2  x  I0-5 

1.2520 

0.0768 

1.728 

1.584 

Vib.  +  B.L. 

4.0668  x  10“J 

1.1935 

0.10420 

1.195 

5.718 

■fa 

1.1709 

0.11403 

1.088 

3.186 

fj/ 

1.1156 

0.11924 

KflKEGSiZiHM; 

2.446 

fH3 

1.1159 

0.14004 

0.888 

3.048 

1.1161 

0.14164 

0.808 

1.638 

22 

1.1077 

0.14663 

0.760 

4.194 

25 

1.0967 

0.13812 

0.630 

1.920 

28 

1.0832 

0.13519 

0.551 

2.214 

30 

4.0668 

x  10- s 

1.0756 

0.13365 

0.S08 

4.208 

10 

7.4703 

x  10-! 

1.1962 

0.10712 

1.228 

— 

12 

1.1683 

0.10961 

1.045 

15 

1.1068 

0.09B60 

0.751 

3.358 

18 

1.1282 

0.14343 

0.910 

2.233 

20 

1.1225 

0.13727 

0.783 

1.779 

22 

1.1132 

0.13007 

0.675 

1.523 

25 

1.0997 

0.12760 

0.582 

1.767 

28 

7.4703 

x  10~5 

1 .0859 

0.12667 

0.516 

1.898 

12 

1.5225 

x  10“4 

1.1695 

0.10356 

3.825 

15 

1.1263 

0.10732 

0.817 

3.750 

18 

1.1313 

0.14225 

0.905 

2.309 

20 

.1.1264 

0.13611 

0.777 

1.867 

22 

1.5225 

X  10~4 

1.1166 

0.12645 

0.656 

1.521 

Vib.  +  B.L. 

23 

1.9  X  10-4 

1.1239 

0.1611 

0.679 

7.302 

Vib.  +  radiation 

IS 

8.7  x  10-4 

1.1109 

0.09562 

0.622 

4.506 

Vib.  +  radiation 
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APPENDIX  F 
EXAMPLES 

The  run  times  shown  in  Fig.  9  were  obtained  from  calculations  made  with  the  input  file 
“ENVELOPE.”  The  sample  shown  in  Fig.  F-l  is  actually  a  composite  which  includes  the 
controls  necessary  to  restart  the  blunt  body  calculation  and  to  continue  it  into  the  Method- 
of-Characteristics  calculation.  All  of  the  computations  used  the  thermochemical  data  set 
NEWCHEM  displayed  in  Fig.  B-2.  The  reactions  and  rates  in  NEWCHEM  are  those  discussed 
in  Ref.  7. 

Consider  the  following  exercise: 

Pari  1.  With  NEWCHEM  assigned  to  Unit  8  as  in  the  example,  set  NEW  =  .T.  to  indicate 
a  new  run,  and  INAFT  =  0  to  limit  calculation  to  the  blunt  body  region  in  Namelist  OPTIONS. 
Save  the  output  file  assigned  to  Unit  9  and  submit  the  job,  allowing  55  sec  of  CPU  time 
if  Cray  X/MP  is  used. 

Selected  pages  from  the  resulting  printout  are  displayed  in  Figs.  F-2a  through  F-2e.  The 
table  at  the  bottom  of  Fig.  F-2a  displays  for  each  species  the  electron  energy  levels  (HCFNT) 
and  the  degeneracy  of  each  level  (FCFNT).  The  ALPIJ  matrix  indicates  the  number  of  atoms 
of  element  i  in  species  j.  The  input  reaction  data  are  echoed  at  the  top  of  Fig.  F-2b.  The 
shock  geometry  calculated  from  the  initial  shock  parameters  appears  next.  The  coordinates 
X,  SZ,  and  SR  are  dimensionless,  as  is  the  curvature  CURV.  The  shock  angle  SIGMA  is 
in  radians.  The  lower  part  of  the  figure  is  filled  with  status  and  error  messages.  MINRAY 
is  the  index  of  the  lowest  numbered  ray  still  in  the  computing  grid  (the  axis  is  numbered 
4)  and  the  number  of  repeated  steps  is  one  less  than  FLOP.  The  error  messages  may  be 
interpreted  with  the  aid  of  Table  11  and  the  species  listed  in  Fig.  F-la.  For  example,  B2(U> 
is  the  fifth  species,  i.e.,  02. 


The  results  of  the  third  shock  iteration  are  displayed  in  Fig.  F-2c.  The  average  deviation 
between  the  desired  and  calculated  boundary-layer  edge  points  is  quite  acceptable  (1.855 
x  10-4),  but  the  relative  error  in  flow  angle  fit  on  the  junction  ray  (MAX  DEV)  is  more 
than  ten  times  the  specified  value  of  RMSTST  and,  therefore,  the  “DID  NOT  CONVERGE” 
message  is  issued.  Note  also  that  there  are  76  records  on  the  restart  file. 

The  results  of  the  70th  successful  integration  step  are  shown  in  Fig.  F-2d.  Only  rays  16 
and  1 7  are  still  active;  95  percent  of  the  input  mass  flux  on  the  other  rays  has  been  accounted 
for.  The  next  six  lines  are  written  by  subroutine  BLAYST  to  report  the  results  of  the  body 
adjustment  procedure  (Section  3.2)  which  moved  the  body  slightly  forward  (compare  the 
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first  value  of  ZBDYN  to  the  ZSTAG  shown  in  Fig.  F-2c)  and  decreased  the  nose  diameter 
by  0.00005  (compare  BDCN(2)  with  BDC(2)).  The  polynomial  coefficients  are  normalized 
by  the  latest  shock  radius  as  are  the  coordinates  of  the  intersection  between  the  junction 
ray  and  the  boundary-layer  edge  —  BZ,BR. 

Parts  of  the  extrapolated  results  at  the  boundary-layer  edge  are  shown  in  Fig.  F-2e.  Ray 
4  is  the  axis  and  ray  17  is  the  junction  ray. 

Part  2.  As  part  of  the  exercise,  restart  at  step  65  of  the  blunt  body  calculation  and  attempt 
to  continue  onto  the  afterbody.  Assign  the  restart  file,  which  was  written  on  Unit  9,  to  Unit 
1.  Set  NEW  =  .F.  in  OPTION.  Set  ITAPE  =  I,  NRCORD  =  65  in  RESTART.  Also  set 
INAFT  =  l  in  SCHANGE  and  choose  ZETAMX  =  0.6  in  SSUPER.  Note  that  SCHANGE 
must  follow  $  RESTART,  and  that  this  pair  can  go  anywhere  in  the  file  between  SOPTIONS 
and  SSUPER. 

The  first  page  of  the  resulting  output  is  shown  in  Fig.  F-2f  and  the  last  page  in  Fig.  F-2g. 
The  Method-of-Characteristics  calculation  stopped  after  completing  7  rays  with  the  message 
“ERROR  at  STATEMENT  11  SUBROUTINE  CFFCNT,”  which  means  that  the  Mach 
number  became  subsonic.  Sometimes  round-off  errors  will  allow  computation  past  the  test 
in  CFFCNT  but  will  stop  the  run  with  the  message  “ERROR  AT  STATEMENT  90  OF 
SUBROUTINE  INTR,”  which  means  that  the  characteristic  did  not  intersect  either  the 
previous  ray  or  the  known  body  or  shock  points.  In  any  event  the  iteration  was  evidently 
diverging.  The  step  size  was  probably  too  large  for  the  rapid  curvature  of  the  body.  The 
message  “ITERATION  NOT  CONVERGING  IN  SUBROUTINE  FELD,”  and  the  following 
line  of  numbers  were  written  by  subroutine  DIGNOS.  Since  the  message  appears  only  once 
per  ray  at  point  number  27,  the  interior  point  computation  evidently  converged  on  the  eighth 
iteration.  To  interpret  the  other  numbers  in  the  message,  look  at  the  arguments  of  the  call 
to  CNVRGE  in  FELD. 

Part  3.  As  Part  3  of  the  exercise,  resubmit  the  job  files  of  Part  2  with  just  one  change; 
in  SSUPER,  decrease  ZETAMX  to  0.4. 

The  last  two  pages  of  the  resulting  printout  are  reproduced  as  Fig.  F-2h.  This  run  was 
stopped  by  subroutine  STIP  because  the  input  time  limit  of  98  sec  was  exceeded.  If  the  output 
file  from  Unit  10  was  saved,  the  run  may  be  continued  by  setting  ITAPE  =2  and  choosing 
an  appropriate  value  for  NRCORD  somewhat  less  than  the  final  value.  Since  the  MOC 
calculation  has  a  good  start,  although  it  is  still  on  the  spherical  nose,  it  is  probably  possible 
to  accelerate  the  computation  by  increasing  ZETAMX  to  0.6.  Once  the  calculation  is  well 
past  the  sphere-cone  junction,  ZETAMX  =  0.8  would  probably  maintain  stability  and  further 
increase  the  marching  speed. 
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The  heat-transfer  (QDOT)  and  FRICTion  coefficients  which  appear  at  the  top  of  Fig. 
F-2h  are  meaningless  since  the  input  of  TWALL  was  neglected.  The  flow  angle  PHI  at  the 
boundary-layer  edge  (NPT  =  0)  is  slightly  greater  than  the  body  angle  (PHI  BODY).  Note 
that  the  average  equilibrium  indicator  —  AV.CHI  —  is  small  and  negative  close  to  the  body 
but  increases  toward  the  shock  indicating  that  the  gas  is  recombining  near  the  body  and 
dissociating  close  to  the  shock  where  AV.CHI  is  positive. 
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NON -EQU1  LIBRIUM  WITH  STANDARD  RATES 

LB  I, NV18  =0. 1NAF  T  » I . NE W-  F  .  I  UNIT- ' .UNIT-'CGS* 
; THIS  IS  A  NEW  RUN 
NON  EOUILBRIUM  AS  IN  SECTION  2.2.1 
;  EQLJI  LBQI'JM  VIBRATION 
; INPUT  ALL  FREE  STREAM  CONDITIONS 
: USE  CGS  UNITS  DTTAULT 
;  I NCLUOE  VISCOUS  0I5»»«  ACCMEN1  OPTION 
: DO  ONLY  BLUNT  BOOT  CALCULATION 


ENVELOPE 

SPHERE  CONE 
SOP  I  IONS 
NEW- . T.  . 

MOO EL  -O. 

NVI8-0. 

IUNIT-1. 

UNIT- 'CCS* . 

INBLR-1 . 

INAFT'O 
t 

SUPSTREAM  VEL-B  09G2E-05.PPES- I  0I32E05 . TEMP -300. . 

VEL-6.0962E-05,  ;CENTERL INf  VELOCITY  IN  CM/SEC 

PRES- 1  OI325EOS  : FREE  STREAM  PRESSURE  IN  DYNE/XM«t2 

TEMP-300.  : FREE  STREAM  STATIC  TEMPERATURE  DEG  K 

FRACTN-32  19490E-05. ;MOLES/GRAM  OF  AR  OEEAULT  VALUE 
72.43929E04,  ;  02  * 

26  96274E-03.  •  '  N2 

t 

S SHOCK  RSMP-0  71527. S:RADIUS  OF  CURVATURE  OF  SHOCK  AT  AXIS  IN  UNITS 

SIP-0.13611.  SPECITICO  IN  OP  11 ONS  NAMELIST 

SI-0.00777 


f BOO Y  REXP-2  .C1-1.2J.C2-1  .MORE-  T  ♦ 

REXP-2  ;  EXPONENT  OF  < R  RjUNC )  TERM  EON  53 

CI-1.27.  ;  TWICE  THE  NOSE  RAOIUSSFE  EO  54 

C2* -  1 .  :  SEE  EQN  54 

MORE-.T.  ;  ANOTHER  BODY  SEGMENT  r  01  LOWS 

SBOOY  RvJUNC =0.61  3933 .  ZJUNC  ■  472792 . RE XP - !  . C 1 •  2642  I  14 . MORE • . F  S 


RdlJNC-'O.  6 ' 3933  . 
Zi  UNC-  .472792. 
REAP- 1 .0. 

Cl-  26421 »4. 
MORE-  F 


: RADIAL  COORD  OF  SPHERE-CONE  JUNCT I  ON -01  MENS I ONAL 
;  A  X  I A  l  COORDINATE  OF  SPHERE-CONE  JUNCTION 
;F!RS1  POWER  OF  R  FOR  CONE 
; TANGENT  UF  HALF  CONE  ANGLE 
; LAST  SEGMENT 


SCONTROL 
N8SIT-3. 

IPRNT-b. 

SKAP^O  5. 

PSJRAT» .05. 

ERRORS  1  0€-6. 

DE l MIN- I  E-9. 

OElX-  05. 

AO-  05. 

GAMTST-1 .£ -07 

methoo-io. 

TIME -98. 

S 

SRESTART  ITAPE 
I T  APE  -  1  . 

NRC0R0-65. 

% 

SCHANGE 
1NAFT-1. 

% 

S  SUPER  2MAX-10  6  .  NMAX-29  .MOOR  A  V  *  ^  .  2E  T  AMX-0  4  .  ME  1 1 IOD- 10.  1  ON- I 


OC -9.GAMTST- 1  E -7 . 


BSIT-3.DELX-  OS.XO-  05, DEL  MI N- 
;  TAKE  THREE  GLOBAL  ITERATIONS 
;  STATUS  REPORT  EVERY  5TH  STEP 
; DECREASE  DAMPING  -EON  50  DEFAULT  IS  .35 
STOP  WHFN  95  OF  MASS  FLUX  IS  ACCOUNTEO  FOR. 

; DEFAULT  CONVERGENCE  CRITERION  ALL  ITERATIONS 
MINIMUM  STEP  SIZE  INITIAL  SIZE* 100-0ELM1N 
; AVERAGE  RAY  SPACING  DEFAULT  IS  0.04 
: F IR5T  OFF  AXIS  RAY  AT  X •  05 .  INCLUDE  AXIS. 

; MI N1 MUM  MOLE  FRACTION  DEFAULT  VALUE- II -7 
; USE  9  CORRECTORS  LOMAX  ALGORITHM  APPENDIX  A 


•1.NRCORO-65  V 

; RESTART  IN  BLUNT  BODY  COOE  IF  NEW-.F 
: RESUME  CALCULATION  AT  65TM  RECORD 


;  IF 


RESTART.  CHANGE  TO  INClUDt  M . 0 . C . CALCUL AT  I  ON 


Figure  F-l.  Input  file  ‘ENVELOPE’  for  nonequilibrium 


run. 
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XLMBDA  -  4  j 1 50984E  *02  HINF  - 

SPECIES _ 

ALPf  J 

O  0.000006*00  0  762i4E*0O  0  108856*01  O 

0  0  0.0.0.0  .  0 . 50000C *0  1  O.SOOOOEHji  O.  10000E*01 

_ N _  O. 000006*00  O  922426*02  0.924956*02  0 

.0.0. 0.0. 0.0.  0.40000E*01  O  80000E*01  0.4000bF*Ot 

E- _ O.OOOOOE*OQ  0 . 000006*00  0000006*00  0 

.1.0.0.0.0  0  .  6  200006*0 1  0. 00000 E *00 0. OOOOOE* 00 

_  *«  _ 0 ■ OOOOOE+OO  0446856*03  0.449776*03  0. 

0.1  00  do.  0.10000E*01  0  500006*01  0300006*01 

_ _ 02  _  O  OOOOOE  *00  0  379B3E  *02  _0_.  63353E  ♦  02  0 . 


.0.0. 0.0  0  0 


MO 

.0.0. 0.0. 

NO* 


O.  30000E  *01  0  20000E*01  0.  1006o€*dl 

*00  0.241 106*03 0.2860 IE *03  0. 

♦01  0 . 30000E  *0 1  0.600006*01 

OOOOE  *00  0. 58052E»00  0.2 1092 E*03  Q. 

O • 20000E  *01  0  20000E *0 1  0.200006*01 

OOOOE *OQ  0  193386*03  0.383586*03  Q. 

0.I0000E*01  0.60000E401  O  300006*01 


2 .  I924662E  *02  8 1  GAM  •  1.400 
FCFNT 

7i122£»Oa  O. 1f| 1 1 E  *03  0 

0  50000E*01  0  1<XXX>E*01 

1 3836E  *03  0. 4 2 264 E +03  0 

0.60000E*01  O. 120006*03 

OOOOOE  *00  0.  OOOOOE  KH) _ 0 

0.00006c«00  O. OOOOOE 400 
OOOOOE *00  0. OOOOOE *00  O 

O. OOOOOE *00  O  OOOOOE *00 

17312E4Q3  0 ■ 33893E  *03  O 

0 . 30000E  *0 1  6  300006*01 

28783E4Q3 _ 0, OOOOOE »00  O 

O. IOOOOE401  0  OOOOOE 400 

22Q32E4Q3  Q, OOOOOE *00  Q 

6. 40000E*0 I  o.oooooc*oo 
35Q?BE*03  0351 69E*03 0 

'OOOOE+OI  0.200006*01 


■  OOOOOE 4QQ  0 

O. OOOOOE  400 
OOOOOE  ♦  00  0. 

O.OOOOOEiOO 
■OOOOOE ♦ OQ  Q, 

O. OOOOOE 400 

■  OOOOOE 4QQ  __0. 

O  OOOOOE *00' 

OOOOOE *00  Q. 
0  OOOOOE *00" 

■  OOOOOE 4Q0  O. 
O. OOOOOE 400 

oggooE*gp _ o. 

O. OOOOOE 400 
OOOOOE  4Q0  O. 
0.  OOOOOE  400 


000006  *00 
o.oooooc 
030006400 


O  OOOOOE *00 

OOOOOE »OQ _ 

0. OOOOOE 4 X 
OOOOOE 4Q0 
O  OOOOOE 400 

OOOOOE *0O _ 

O.OOOOOC 400 

000006*00 _ 

0.000006400 


a.  Output  from  nonequilibrium  run  ‘ENVELOPE’ 
Figure  F-2.  Sample  calculations. 
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I 

1  J. 

0 

REACTION 

0  0 

0 

At  RK(t> 

2. lOOv  .JEtiS  -0  5000 

CRM) 

0.5938C-05 

OKI  I  ) 

1.000 

AM.  1) 

1  .0 

AM. 2 
.00 

I  T* 

DIR 

FWO 

MU-IOOOIOOO  MUP ■ 

3  0 

0 

00000 

2  02 

02 

0  0 

02 

3  5GOOOOOE  *2 1  1  500 

0.5938E-05 

1  oco 

00 

00 

W 

F  WO 

MU-  00002000  MUP- 

2  0 

O 

01000 

3  02 

M 

0  0 

M 

1  1900000F  *21  1  500 

0 . 5938E  *05 

1.000 

1  .0 

00 

2 

F  WO 

NU-  0  1  1  1  1  1 

1  1  MUP- 

2  1 

1 

1  0  1  n 

4  N2 

N 

M  N 

M 

7  5000000E  *20  -1.500 

0.0000 

1.000 

00 

1.0 

REV 

NU-  0  1  0  0  0  1 

0  0  MUP- 

0  3 

0 

00000 

5  M2 

M2 

N  N 

M2 

1  5000000C *20  1  500 

0.0000 

1.000 

00 

1.0 

w 

REV 

NU"  000003 

0  0  MUP- 

0  2 

0 

0  0  10  0 

6  M2 

* 

N  N 

M 

5  OO0OOOOE  » 1 9  1  500 

0.0000 

1 .000 

00 

1.0 

2 

REV 

NU-  1  0  1  1  1  1 

1  1  MUP- 

1  2 

1 

110  11 

■  ■ 

7  NO 

M 

N  0 

M 

5  1 900000E ♦ 2 1  1.500 

0 . 7549E  *09 

1  000 

00 

.00 

0 

F  WO 

NU-  0  0  0  0  0  0 

1  0  MUP- 

1  1 

0 

OOOOO 

a  n 

02 

NO  0 

1  OOOOOOOE  *  1 2  0.5000 

3120 

1.000 

00 

.00 

W 

TWO 

NU*  0  1  0  0  1  0 

0  0  MUP- 

1  0 

0 

0  0  0  10 

9  0 

M2 

NO  N 

5  OOOOOOOE  *  1 3  0  OOOO 

O  3802E*05 

1.000 

.00 

00 

w 

rwo 

NU-  1  0  0  0  0  1 

0  0  NUP- 

0  1 

0 

0  0  0  10 

10  N 

0 

NO*  E- 

5  4000000E  *21  1 . 500 

0.0000 

1  .oco 

.00 

00 

w 

REV 

NU-  1  10000 

6  0  NUP- 

0  0 

t 

0  0  0  0  1 

SHOCK  GEOMETRY 

M  xTST  sftSl  S»TiO  sioiA  curvTmT 


1 

2 

0  500006-01  O  I2498E  02 
0.10000  0  4996BE02 

0. 49979E -01 

0  99833E  01 

1  5208' 

1.4708 

•1.0003 

•1.0010 

3 

4 

O  15000  6  1 1233E01 

0.20000  0  19945E01 

0. 14944 

0  19867 

1 .4207 

1.3706 

1  0018 

1  0014 

5 

6 

0  25000  0  31110E-01 

0.30000  0  44893E-01 

0  24740 

O. 29551 

1  3206 

1 . 2707 

-0  99927 
-0  99437 

9 

8 

0.35000  0  6085  IE-01 

0.40000  0  7B922C-OI 

0  34289 

0. 38943 

i .  22 12 

1  .  1722 

-0  9B57B 
-0.97280 

9 

10 

0.45000  0  99414F-0I 

0.50000  0  12210 

0  43502 

0  47959 

1.  1240 

1  0768 

-0  95501 

0  93231 

1  1 

12 

0.55000  0  14681 

0.60000  0.17347 

0.52305 

0.56534 

1 .0309 

0  98638 

-O  90494 
-0  87343 

13 

14 

0.65000  0  20194 

0.70000  0  23212 

6  60644 

0  64630 

0.94357 

0.90257 

-0  83851 
-0.80103 

19 

18 

0  75000  0  26386 

0.80000  0  29705 

0  68493 

0.72232 

0.86349 

0.82640 

-078188 

-0.72191 

17 

18 

6  85000  0  33157 

090000  0  38731 

0  75B49 

0. 79346 

0  79130 

0.  75820 

-0  68189 
-0  64244 

19 

20 

0  95000  0.40415 

1.0000  O  44199 

6  82726 

0  85994 

0.72704 

0  69777 

-0  60408 

0  56719 

NRCORD  •  5  FLOP  •  1  MINRAY 
NRCORD  •  10  FLOP  •  1.  MlNRAv 

4  Y*0  31000E-050. 160OOE-O5O  OOOOO  LAST 

•  4  Y-0. 10230E-030.51200E  040  OOOOO  LAST 

STEP  -0  OOOOO 
STEP  -0  OOOOO 

PSI(MMIN). 

PS ! ( MM IN)* 

82 (  11)  NEGATIVE  ON  RAY  4RUNGE - 

KLTTTA  STEP  10  V  •  1  53BE  03  DELY 

-  8  I92E-04 

NRCORD  -  15  FLOP  •  2.  M I MP4 V 
NRCORD  •  20  FLOP  •  2.  MI MR AY 

•  4  v-0  13873E  020  460806  030 

■  4  Y*0  10497E-010. 34992E -020 

OOOOO  LAST 

OOOOO  LAST 

STEP  -0  OOOOO 
STEP  -0  OOOOO 

PSI(UUIN)- 

PSHMMIN)* 

82 (  11)  NEGATIVE  ON  RAY  4PUNGE- 

KUTTA  STEP10  V  •  1.5501-03  DELY 

-  5.000E-03 

B2<  11)  NEGATIVE  ON  RAY  4RUNGF. 

KUTTA  STEP  10  Y  •  1  B34E-02  DELY 

-  4. 147E-03 

NRCORO  •  29  FLOP  •  4.  M|NP/\V 

•  4  V-0.70047F  O10  303751  -020 

OOOOO  EAST 

STEP  -0  OOOOO 

psi(mmEnT^~ 

B2<  13)  NEGATIVE  ON  RAY  4RUNC.E 

KUTTA  STEP  <!  V  -  J  TIRE-03  DEL* 

■  3.1 idi-03 

B2<  1l)  NEGATIVE  ON  RAY  4 BUNGE 
NRCORO  •  30  FLOP  »  6.  MINRAV 

KUTTA  STEP  10  v  ?  5 1 BE -02  DELY 

•  4  Y-0  75**  151  OIO.  173b4E -020 

'  2  624E  03 
OOOOO  LAST 

STEP  -0  OOOOO 

PSI(MMIN)- 

B2<  13)  NEGATIVE  ON  SAY  4RUNGE • 

KUTTA  STFP  9  V  J  OISE-02  DELY 

■  2.  187E-03 

B2(  11)  NEQATtVE  ON  RAY  4RUNGE - 

KUTTA  STEP10  F  •  3  26BE-02  DFLY  ' 

■  1.7281-03 

b.  Output  from  nonequilibrium  run  ‘ENVELOPE’ 
Figure  F-2.  Continued. 
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NBCORO  ■  35  FLOP  ■  9  MINPAY.  4  Y  0  1  K*4nf  0 10  SSOli'Jf  030. OOOOO  LAST  STIP  *0  00000  PSIINMIN)* 

P2I  Bl  NEGATIVE  ON  RAY  4PUNGI  KIJtlA  STEP *0  *  3  5GSF  Or  DELT  ^1,362EH}3  - 

NProso  -  40  riCP  3  MINPAc  4  » »0.  3*  1O4E-0I0  747-1E  030.00OIYO  LAST  STEP  -0  00Q0Q  P5HW4IN)  • 

E3I  Ml  NrGATIVE_ON  Rt>  4RUNGI  Kbit*  SltPIO  V  3  7BII-0J  RELY  •  1. 0781-03 

NRCORD  •  45  FLOP  •  10  M1NPA  y *  *  »  0  3ltd76=  -liin  0  M  15T  030  O00M0  LAST  STCP  >0.00660'  PSlfWlN)- - 

B3i  ni  Mrr.ifTvfoN  ray  aponk  pmii*  step  to  y  .  3  him  or  iyely  -  7  pi  if  <>4  -  - 

NRCORD  50  PLOP  *  11.  MINPAY-  4  V  .p.  4O102E  •  nio  4 74  1  1 i  030.00000  LAST  STEP  •0.00000  P5I(MM1N)= _ 

831  lit  NEGATIVE  ON  RAY  4HLML.I  KU  M  A  Sit  J  •  4  1131-03  IYELY  -  5  693E-04 

NHcr.no  5"¥  FLOP  13.  MINFAYy  4  Y  O  41N43E  OIO.399.3l  -030.00700  LAS1  SIEF  -O  OOOOO  OSTIHHINT* 


831  lit  NI  GAT  1 VF  ON  PAY  4HL4YK -I  Kill  1 A  S1LPH.I  .  4  3't3E-‘V  IHIY  -  3  AS  IF  -04 


NHCURO 
NRCC1P0 
NRCUHO 
NRCORD 
NRCORD 

NHCMRII 


OO  FLOP  13.  MINPAY-  /  V.0.42NG8E  :i  1(1 . 23it(.3t  030 . OOOOO  LAST  STEP  -O  OOOOO  PSI  (MMI N)- 

65  FLOP  I  (  MIIRAY.  10  Y-0  14M9 1 T  010.5300.11-030  140631-0'  LAST  SIFP  >0  1371BE-01  PS  1  ( ifiii  ill- 

70  riot  l.i  PINPA.-  13  Y-0.  I  tier  010. STI17E -030.2 I4B3E-0’  LAST  STEP  -O  33833E  01  PSIlMMINl" 
75  FlOP  11  I  IIP.Y.  15  Y-o  •(  73E  OIO  i0974f  030  5IA10F-0'  LAST  STEP  -0  42974E-01  PSl  tiSfllili- 

80  F  L  DP  11  RtM.AY-  IH  Y-l*  49l|  no  1014  3 1  030  53  I44E -0 1  LAST  SI  EP  -0  57769E-OI  PSI(MMIM). 

B4V  (.3  |W<5  <1073(1  0I0CIYI3)-  - 5 - 


ITERATION  ho.  |  I  UI4  RAPIDS  o  IIS fi  r.  IP  n  l.mil  si  6  oil 77 7 

RMS  0EVIA1 ION 4  3646334E04  MAA  PE V  t  4433B56E  03  2MAG  -  4  43I6737E  02 _ 

_ NIK  SHOCK  RADIUS  ■  7_  14 1551  IE  01  NEK  SHOCK  pahame TERS  SIP*  I  .  383308 IE  -  01 _ SI-  7  8921060E  03 


SHOCK  GEOMETRY 

M 

MMI 

SZIM1 

SIGM4 

CURV(M) 

1 

U. 500001 

CM  0.  124«flE>02 

0  49979E  01 

1  5200 

-  1  0004 

2 

0.  IOOOO 

O  49970E-02 

0  9'IAHF  01 

1  4  707 

1  0013 

3 

0  15000 

o  11234E-OI 

O  14944 

1  4207 

1  0022 

4 

0-20000 

M  19949E  C)  1 

O  19967 

1  3705 

•  0023 

5 

0  25000 

0.31«l7f  01 

O  24740 

1  3205 

1-0006 

6 

0  30000 

M  447091  01 

O  29551 

1  2705 

O  99617 

7 

0.39000 

O  60678E  01 

0  34  289 

1.2209 

-0  98799 

B 

0  40000 

0  789661  01 

O  38942 

1  1 7  IB 

-0  97535 

9 

0  45000 

0  99500E  01 

0.43500 

1.  1235 

-0  95781 

10 

0.50000 

O  12219 

0  47955 

1  0761 

-0  93526 

1  1 

0  55000 

0  14694 

0-53799 

1.0300 

-0.90793 

12 

0  60000 

O  17363 

0.5G52B 

0  98538 

0  8^635 

_ _ _  13 

0  65000 

O  20215 

0  60633 

0.94243 

-0  84128 

0  70000 

O  23238 

0  64615 

0.00130 

•0  80358 

15 

0  75000 

n  26417 

0  68474 

0.86210 

-0-76416 

16 

o  aoooo 

29742 

0  7 2 208 

0  87400 

O  73390 

17 

0  85000 

0.33200 

0  75819 

0.78971 

-0.68358 

n 

0  90000 

O.  36779 

0.79310 

0  7S653 

-0  64383 

19 

0 . 95000 

O- 40468 

0  02684 

0.72531 

-0  60520 

20 

1  0000 

0  44259 

0  B5945 

0  69598 

-0  56805 

-  NKCORO  •  .  5  FLOP  ,  MtNHiY- 

A  VO 

3  IQCOt  *050  I6000C  050 

OOOOO  LAST 

51 EP  -0  OOOOO 

PSI(MMIN)- 

NRCORD  •  10  FLOP  1  MINRAY- 

4  vo 

10230C  030  5 1 201*  040 

OOOOO  LAST 

STEP  -0  OOOOO 

pshmmin)- 

B2(  11?  NFGA I I VE  ON  BA  >  4RUNGE -KUTTA  STEP  10  V  =  I.63BE-03  OELV 

H  I97E-04 

NRCOBO  -  15  FLOP  •  2  MINRAV- 

4  VO. 

13823E  020. 49000E  030 

OOOOO  LAST 

STEP  -0  OOOOO 

PS  I  ( MM  IN)* 

NRCORD  ■  20  FLOP  •  3  MINRAY- 

4  VO 

I0497F  HO  34992E-030  OOOOO  LAST 

STEP  -0  00000 

PSilMPHN). 

B3(  III  NEGATIVE  ON  RAY  4RDNGE  KUTTA  STEPIO  V  •  1  55OF02  DETy '•  5  OODE-03 


B2(  III  NEGATIVE  ON  RAY  4RUNGE  KUTTA  SIFriO  Y  5  I.834C-03  DELyTTT  147E  03 - 

NRCOBO  •  35  FLOP  ■  4.  HIINRAI-  4  V-0  31I047E-010  30375E  030  OOOOO  LAST  STEP  -0. OOOOO  PSI(MIHN)- 

83 (  131  NEGATIVE  ON  RAY  4HUNGE  KUTTA  STEP  9  Y  2. 3165-02  PELY  •  3. 11QE-Q3  _ 

B2  I  NEGATIVE  ON  RAY  ARUNGF  KUI I*  SltPIO  V  •  2.5I6F  03  t  •  J  624E-03 


b.  Concluded 
Figure  F-2.  Continued. 
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B2I  II)  NEGATIVE  ON  RAY  4RIJNCr  KUTTA  STFPIO  Y  *  3  5IRE-02  DELY  -  2.62AE  03 

NRCORO  ■  TO  FLOP  •  «■  HINPAY-  _A  Y-0  259A5E -Old  I736AJ -020  OOOOO  LAS)  STEP  -O  OOOOQ  PSIQMIN)- 

B  3)  NEGATIVE  ON  BAY  A  RUNG  I  -KUTTA  STEP  9  Y  -  30I6E-02  I  .  -  3 .  U7E-03  _ 

B2I  II)  NEGATIVE  OAI  RAY  4RUNGL  KUTTA  ST1PI0  V  •  3.36aE-03  DELY  ■  I  72RF-03 

NRCORO  •  35  FLOP  •  8.  N1NRAV-  4  Y-0  3l045F-r.[0  B5050E -030  OOOOO  LAST  STEP  *6  OOOOO  PSI(MMIN)- 

"557  ■)  NEGATIVE  ON  5AY  ARUNGE  KIJTTA  STEPIO  Y  •  3.565E-02  DELY  *  I  3S3E-03 

NRCORO  •  40  FLOP  - _ 9.  M1NRAY.  4  T-0.3S104E  010  74751E -030.00000  LAST  STEP  -0  00000  PSI1MMINI- 

52(  II)  NEGATIVE  ON  BAY  4RUNGE-KUTTA  STEPIO  Y  •  3  7B8E  02  DELY  •  I.025E-03 

NRCORO  •  45  FLOP  ■  10.  &INRAY-  4  Y-0.  3B070E -0)6  6 14  )6E  -030.  OOOOO  LAST  STEP  -O  OOOOO  pTKHHIN)- 

B2(  II  )  NEGATIVE  ON  RAY  4RUNRF  KUTTA  STEPIO  V  -  3  H9T  CJ  OETy"-  7  43  IE -04 

NRCORO  •  50  FLOP  -  II.  MINRAY »  4  Y -0.40 IQ3E -0 10.4794 It -030. OOOOO  LAS)  STEP  -0  OOOOO  PSIHWIN)- 

63(  IT)  NEGATIVE  ON  RAY  4RUNRE -KUTTA  STEP  B  V  ■  4  IIBE-02  DELY  -  5  S93E-04 

NRCORO  -  55  Flop  -  12  MINRAY-  5  Y-O.4I643E-0IO  1S94 2E -030  OOOOO  list  STEP- ^0  OOOOO  PSI(MNIN)- 

NRCORO  •  60  FLOP  •  13  MINRAY  11  Y-0. 44266E -O 10 . 77875E -030  I6330E-01  LAST  STEP  -0  I7967E  01  PSI(MMIN)- 

NRCORO  ■  85  FLOP  -  12  MINRAY-  14  Y-O  4B054E-010  4235BE  036  34917E-01  LASt  SfEP-0  39B 1 3E  :01PS  t  (KNIM- 

NRCORO  •  70  FLOP  *  12  MINRAY-  16  Y-O. S35S0C -010.  I7339E -030. 66370E -01  LAST  STEP  -0  5510OE  0I  PSIIMMINI* 

NRCORO-  t8V75)-O.S93  ITE  -  lltlELYi  3)"- 

ITERATION  NO  3  SFA1CK  RADIUS  O  71321  SIP  ■  0  13775  O’-  6  007M 

RMS  DEVIATION  -  I  .  BS5 1 394E -04  MAS  of  v  I  096R4I3E-02  2STAG  ■  4  43253435-08 _ 


PROGRAM  DID  NOT  CONVERGE 


VALUES  A I  INTERSECTION  OF  HAY  KITH  BCOV 


*  AY 

ZB 

PR 

PH  IB 

DELTA  _ 

FLUX 

5 

4  56072E-02 

4  77592E  02 

1 .517 13E+00 

2.B2715E-03 

5. 1624 1 E -04 

6 

4  94480E02 

9  53713E  02 

1  46347E*00 

2.B6453E  *03 

2 . 0554 2 E -03 

7 

5  58339E  02 

1  4269IE01 

1  409B4E  *00 

2 . 90300E  *03 

4  603991  03 

B 

6  4  742  IE  02 

1  B9576E'-01 

1  3S623E«00 

2 . 95634E -03 

8  13373E-03 

9 

7  81416E0? 

2.35007E  01 

1  30266E  *00 

3 . 02863E  -03 

1  2607  IE -02 

IO 

0  99945E  02 

?  8I490E-01 

1  249 I2E+00 

3  l?337f-03 

1  79765E-02 

1 1 

1  062S6E  01 

3  262LHE  Ol 

1  1 9562E  *00 

3  24328E-03 

2.4 1B47E-02 

12 

1  24076EO1 

3.70O65E-O1 

1  1 42 1 5E *00 

3 . 39549E -03 

3. 11659E-02 

13 

1  4 5 800 E  01 

4.  12792E  01 

1  -06171  §400 

3  58092E-O3 

3 . 8844  IE -02 

14 

1  6B966E  01 

4  54325E  01 

1  03529E  *00 

3. 8029  IE -03 

4.713S4E-02 

!5 

1  9431  it  01 

4.94B52E01 

9.BIB78E-01 

4  0650 IE -03 

5.59476E-02 

16 

221765E  01 

9.333651-01 

9.29475E-01 

4  367231-03 

6 . 518 1  IE -02 

17 

2  51253E  01 

5.70682E  01 

B  75067E-01 

4.700771  CJ 

7.4727BE-02 

10 

2  63699E  6i 

6 . 06339E  01 

B  2 1646E  *01 

5. 1 13BBE-03 

0  44B27F  02 

c.  Output  from  nonequilibrium  run  ‘ENVELOPE* 
Figure  F-2.  Continued. 
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STEP  SIZE  •  1  236616  03 


NUMM'R  Of  SU<  CLSSTUI  L  STEPS  *  70 


DISTANCE  FROM  SHOCK  *  3 . 8  19236 -02  ~CM 


RAV  Z  CM 

_ HTOT  ERQ/GM 

MOL/QM  0 


“  R  CM 
S  CAL/QM  K 
N 


P  OYNE  /CM • *  7 
AVER  CHI 
gr 


i  crn  k 
ft  liens  */cc 

AR 


PHD  T.M/CC 
MACH  NO 
02 


XMW  CM/MOL 
FLUX  RAT|0 
N2 


VTflt  CM/SIC 
RAU  ERQ/GM  SEC 


7lW1w“ 

PJI  CM/SEC 
NOA 


1  SS64476E-0I  3 . 82069786 -O f  2  2325419E«07  0  39I0102E»03  9.74I3745C  04  2  3»85549E7oT"  3~2834429l«05  S  2895499E+01 

J  88A5574g*  l  »  2  57069306-00  -3  7726495T  01  4  79&2B4SEH5  1  92700281-00  1.01 89992E  +00  O  OOOOOOOE  *00  2  42I2094E-OO 

1.32154591-02  3  97  159  ME -03  8  1702828F  06  3  2194900E  04  8  1 9 1 2302 f  -OS  2  4422657E -62”T  10040406 -0T  8  77028266 -06 


-LZ - L  4  _,09L®99Ll?J  -- ‘°7  ^  J597!07€*03  9.28B205SE-04  7  3423526E  *0  1  3  500993CE-05  5 . QO66035E  *-0  I 

1  8852  1 34E  ♦  1  1  2  55131B1E-00  -5  2232I49F  Ol  3  62753 10E«15  2  10155841*00  1  54689196*00  0  OOOOOOOE *00"  4  2 133TT9FMX) 

- 1  31468536-02  3.  16718BOE-Q3  6  4821975E-06  3,21949006-04  9  84*44581  05  2  48Q7Q5BE-Q2  1 . 1376942E-Q3  § . 482 1975E  Q6 

”  BOD»  POLYNOMIALS  BEFORE  ADJUSTMENT 

_ IBP* _ TBOT _ BOC(L)  _  _  _ _ _ 

- Q_- PPO0000E*O0 _ 0. 00000006*00 _ 0  OOOOOOOE  *00  I  7806877E+00  -  1  OOOOOOOE *00  O. OOOOOOOE +00 _ O.  OOOOOOOE 400 

- MBOMIE-OI _ B,6080545E-0I  _ p  OOOOOOOE  *00  2  642114QE-01  0. OOOOOOOE*OQ  0. OOOOOOOE *00  0.  OOOOOOOE  ♦OO 


BOOT  POLYNOMIALS  AFTER  ADJUSTMENT 

bdcnTL) 

WOOOOC^OO  f.  77985656 *O0  -  f  OOOOt 
000660€*00  26421140601  dTOOOOC 


O .00000006*00 

o.ooooboocioo- 


d.  Output  from  nonequilibrium  run  ‘ENVELOPE’ 
Figure  F-2.  Continued. 


AEDC-TR-88  31 


M 

28 

im 

ARC 

PM1B 

~  *4 

4  450953nr  02 

o.of"  n*«x>r  *00 

O  OOlXURNlEiOO 

1 .570796  JE* 00 

B 

4  5607 195 E  0? 

4  7 7 S9 1651  02 

4  7771777*  02 

1  .5171 293 E  K)0 

6 

4  944B045F-0? 

o  v*vi34r.r  02 

9.  '.538026*  02 

1  463473OE*00 

7 

5.56339221  02 

1  4?'.?10l  >1 

1.43287301  01 

1  409837  IE  *00 

~B 

6. 474 207 ?f  02 

1  7*. 71M  01 

1  ••  1010791  01 

1  35622951*00 

9 

7.614162‘Jt  02 

2  35  Pl*v«l  ot 

?.:0'70428EOI 

1  3026556E  *00 

10 

B  99945301  02 

2.1'  1  1904  l  Ol 

:»  .»nr»4*it  ei 

1 .2491 IBOE+OO 

11 

1  0625620E  Ol 

3  :T.’5769E  Ol 

1  t:**»'»4  f?E  -01 

1  1956 1G2E  +00 

“12“ 

1 . 24B764 1 1  0» 

3  «i'OC4orr-Qi 

1  r  15M47.7C-U1 

1  14214741*00 

13 

1  -4579982E  01 

4  i.  oi 

1  791701GE  01 

1 .088706  IE  *00 

14 

i  _689662“*F  01 

4  54  f240‘1  Ol 

1  707271  -01 

1  0352B55E ♦ OO 

IB 

1  943l095r-0l 

4  9455215101 

5  2427270E  -  01 

9  81B77BBF  01 

fa 

2  .'2 1784501  -01 

5  .73  105  7"E  01 

5  71B1422E  01 

9 . 2847473C  01 

_  17  51079161  01 

5  701952 /E  01 

G. 103596OE  01 

B. 7506723E  01 

ta 

2  8214342E-01 

r  OOJ19B6F  Ol 

r.  089 1 74.1T  of 

8  31121171  0» 

M 

DELTA  _ 

nniBAR 

000  T 

ENRAD 

4 

o.oooooooc+oo 

0  OOXKXXjfOO 

f»  tHJOOOOOC  »C0 

•)  .  OOOOOOOE  *00 

“  5 

7  S271501F03 

•  1‘.  •40,4"F  -04 

;  .1 19805BE*  l  1 

"  OOOOOC'OL  « OO 

6 

2.8645267E  03 

2  054  7530C  03 

7.27466BOE* 11 

0  OOOOOOOE  *00 

*-  7 

2  903000*  H  03 

4  KW3720E  03 

7  19985941  •  1  1 

0  OOOOOOOE »oo 

ft 

2  956339U  03 

P  13  i60G4f  03 

/.096316f-Ef  1  1 

0  OOOOOOOE  *oo 

9 

3.02863331  0? 

1  .’^20041  02 

f  9flfi0ftB9F  *  1  1 

0  OOOOOOOE *oo 

10 

3.  1233  1591  03 

1  «r.8942E  02 

f  fU/M0ii3E»11 

0. OOOOOOOE* 00 

1 1 

3.2432597E03 

2  4  1-40461  07 

r.  0245093EMI 

0  OOOOOOOE +00 

•w.  13 

3  39549161  03 

3  H5H73E02 

i  4  I90I.50E*  '  1 

o . ooooooof  *  oo 

13  ' 

*3  5805186E-03 

3.8n.’4b;OE-02 

“o  I92I279E41I 

6 .  OOOOOOOE  »oo 

14 

3. 6029 14  IE -03 

4. 7UB362E -02 

•»  94646  1  BE  ♦  1  1 

o. oooooooe »oo 

15 

4  06  5006 8  E  03 

5  591I369E  02 

5  6B4tf887E»li 

0  OOOOOOOE »oo 

16 

4  367 2 306 E  03 

6  -i  1 3262 1 E  02 

5  4  1037  IDEM  1 
r>  l2SO.ir.oE4  1  1 

0  00000006*00 

17 

4  70077011-03 

f  46G3995E  02 

0  OOOOOOOE *oo 

IS 

5  1103486E  03 

B  4399288E -02 

4  B34B615EH1 

0. OOOOOOOE *00 

FRIC7 

O  OfXXXXXJt  *oo 
I  96796941 *03 
3  9327244E03 
5  990799  IE -03 
7B3B4Q79E-03 
9  77 10B38E-03 
I  16B3439E  02 
1  356B769E -02 
1.54 196  1BE  02 
1 . 7225950E -02 

1  B977251E-02 
206618171 -02 

2  22664 1 4E  ^02 
2  3776744E  02 
2  5 IB0581E -02 


TV 

GWAll 

0  OOOOOOOE *00 

0  OOOOOOOE *oo 

O  OOOOOOOE. 00 

0  OOOOOOOE  "00 

0. OOOOOOOE *00 

0. OOOOOOOE *00 

0  OOOOOOOE  *oo 

0  OOOOOOOE  +00 

0. OOOOOOOE +oo 

0  OOOOOOOE ♦ 00 

6  OOOOOOOE  *oo 

0  OOOOOOOE *00 

0. OOOOOOOE  400 

0  OOOOOOOE *00 

o.  oooooooe .00 

0  OOOOOOOE +00 

0. OOOOOOOE *00 
0. OOOOOOOE 400 

0  OOOOOOOE *00 

0. OOOOOOOE *00 

0. OOOOOOOE *00 

0. OOOOOOOE *00 

6 . OOOOOOOE  *00 

0  OOOOOOOE *00 

0 . OOOOOOOE *00 

0  OOOOOOOE +00 

6. OOOOOOOE +00 

o  oooooooe. oo 

0  OOOOOOOE *00 

0  OOOOOOOE *00 

d.  Concluded 

Figure  F-2.  Continued.  — * 
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VALUES  At  ROUNOARY  I.Avffl  FDG6  6 1>|  LOU 


step  m  71  i  pnne  is  oi 


NUMBER  or  SUCCESSFUL!.  SIFTS 


DISTANCE  FRIIM  SHOCK  -  3  e  1*3331  f!J  CM 


RAV  *  2  CM 

MTOT  CRG/GM 


B  i.M  P  DVNE/CM*  *2  T  DfC,  K  RMU  GM/CC 

S  CAL/GM  K  AVER.  CHI  IL.OFNS  A/CC  MACH  Ml) 

COMPOS  I  T  I CN  IN  MCI  FS  PFR  GRAM 

N  €  _ _ AR  na 


AMW  GM/MOL  VTOT  CM/SEC  FLOW  A NO 

FLU*  RATIO  RAD  ERG/GM  SEC  PSI  GM/SEC 

N2  _ NO  _  NO  * 


4 

3  1744542E • 02 

D.OOOOOOOE  *oo 

4  *  935464E  *07 

7 . 69630096*03 

1  99 1568  IE*  16 
3.2 1949006  04 

1 . 4G45979E-03 

2.  14328456*01 

1  87751 1QF*  1  1 

1  3462238E  02 

2  69315056*00 
107481816-02 

-3  0623044F  04 

7  353347HE-05 

9  013461  IE  -04 
3.74033796-05 

1 . 00000006 *00 

2.  1 1  1324J6-Q2  _ 

6  6dooo66£*ob  6 

9  27280006  04  2 

.00000006*00 

* 

3  054/773E  05 

1  07759O5F*1» 

3.4 16I276E -02 
2.69213456*00 

4  175 184 1 E  *07 
8.6274630E-04 

7  88760406*03 

1  10O9725E  03 

1 . 547 107 JE -01 

3  78503 »3E  05 

2  1  447  J8  4  f *01 
l.0000000€*00 

2. 1 14325  IE -02 

?  98*901  IF  *04  7 

0  OOOOOOOE*CO  1 
9 . 284516  IE -04  2 

97187706*0* 

1796  14  IE  02 
.  339101  IF  OS 

1  34600 1 SE  OJ 

I.0B87135E-O7 

2  339101  IE -05 

3.7  1949006-64 

6 

3  32342656-0? 

6  82235876-02 

4  12597166*07 

7 . 06? J997E  *03 

2  14895316*01 

1  ooooooof*oo 

2  12 30067 E -02 

5 . 922568 3E *04  8 
6  00000006*00  7 
9  3337639604  2 

D 1  ItffftMl 
*9580436-02 
298C788F  05 

1  8778287011 

I  34S638G6  *0? 

2  G092877E  *00 
I.G508989E-02 

2  33004 I6E  02 

2 . 29807B8E -05 

1  92 7006  IE •  16 
3.2  I9-J900E  04 

3.0802R46E-01 

3. 7557G22E-05 

7 

3 . 77 10863E -02 

1 .  B782237E  Ml 

1  0208758E-01 

2 .  684  4  759E+-00 

4  044  1 78 1  E*07" 

8  533501  IE  03 

7  62016376*03 

1  849J 1 J3E ♦ 16 

1 . 3762065E -03 

4 . €2647 JOF -01 

2  15598  lflE*Ol 

1  00000006*00 

8859553bE*64  j 
0.00000006*00  1 

H492197E*0I 
6789417L  01 

1.34471ROC  O: 

1  C218I32E  0? 

2  2302 11 9E  05 

3.2 I94900F  04 

J  97884  If. t  r.5 

7  1377R29F  -0.* 

9. 30/90081  04 

/  ?"V\'  11lL  *.•> 

a 

4  !7«>76«G9t  **2 

1 .8787H72E  »  I  I 

t.  34  at,  no  it  o? 

1  .^fflOOF  OI 

2  67779726*00 

«  8  1  13C44C  <yr 

3  0r>O5B48E»O7 

1 .  781  .-■/Fine -a? 

2 .  136729CF  05 

7 . 560  ~ lU9E  *n3 

1  /  1  J  23931  1  10 
3.7  I9.:uooF  04 

1  35413076  03 
G  IH62150F  01 
4  0034638F • 05 

7  18574301  *01 

1 .00000006*00 
J  15/ 1G14F -OJ 

t  i :  B73U •  •  ■  >5 

0  G<)CUUU0C*O3 

9  4n»,20«JOC  -O4 

/.OVK45I.T  *i/1 

2  9*;.r.07|ir  oi 

S.  I3G729CH  -05 

9 

f>  202 -VIP  IE  05 

1  87942(1  II  ♦  1  1 

1 . 34?  1 4GM.  0? 

\  aismAE-oi 

2  6G93P02C « OO 

3  10140126  03 

3  78*.  104  16*0* 

3  .  T.  30G4H6  i.'.* 

2  01925  I3C  09 

7  40238456*03 

1  6 1 24  /56L ♦  10 

I  .719490OE  -04 

1  .77539636  OC 
/.  781448  tk  01 
t  70181.471  -05 

7.  178  IF.  . . 

1  OCKX^XIol  «oo 

J  18208566 -07 

1  4LG84?*L*05 

0  00<XXX>OE*CO 

9  OJ  1735.36  04 

7.331  15  . *76  *01 

4. 59725  1. If -OI 
2.01925126-05 

10 

6  1B42v«'Jf  0? 
i  8001075:*  it 

1  340J492E  02 

2-OI48Gt.lt  01 

7  659 1 38  ?f  «UO 
B.69G22GOF -03 

3  6HM0tC*O7 

4  **0762176 -02 ‘ 
1 . BH00469F-05 

7 . 3057275E  *03 
'  46  1  194 tit  *  16 

3 . ? 194 900E  C4 

1  28997B7E03 
9  .153604  IE  01 
4  4278 Jim  05 

9*3*3tt341  i'i 

1  .OOUOOOOf  tM 
7.7116727F-0? 

1  7S15473FIU5 

o  oooonoonoo 

3  77009086  04 

7  04460951 *01 

6  555 12731  fl| 
88O04  G9E  O'j 

0 

7  34U69I3C  03 

1  88099556*11 

1  33R*6R*6  O? 

2. 335502  IF -01 
2.64737126*00 
8  00G2368E  03 

3  4 12990  IE *07 
-C  26G1077F  07 

1  72220176  05 

1  2690  1 101*03 

1  206O766E ♦  16 

3 . 7 194 9006  U4 

1.7-181  14  41-03 

1 . 09632896 »C0 
4 . 76G3786F -05 

2  ; 1039476*01 

1  OOOOOOOF *  00 
2. 7453 1976 -C2 

2  03049336*05 

0  00000<X>6  *00 

9  93G778U6  04 

f  751*11046*01 

8  81676911  -01 

1  7272017F  Ot 

12 

8  .  G602256E  *02 

1  88 I7827E  •  1  1 

1  31558606  02 

7  6493626F-0I 
7.63-1  IG57I  *00 
7.752H8I If  03 

3  193 190OC  *  07 
9.0684/941  02 

1  549/G13F  05 

7  .  134053  IE  *03 

1  12080581 • 16 

3  ?  1949006-04 

1  70035 17C  03 

1  75905186*00 
5  I5HI98AF  05 

7  2297O4t.f*01 

1  .00000006*00 
2  7*22.387 6  0? 

2.  30200266 *05 

0  00000006*00 

1  111232776  03 

G  45  18  1576  *01 

1  1 3b404  71  HK1 

1  54976  IP!  05 

13 

1 .0U>3O3r.r  r>» 
i  8825S09EMI 
1.33286196  Os 

2 . 9554 OBOE  01 

2  61970176*00 
0.44946396  03 

2.R5F27«ME*07 
-  1 .30952721-01 

1 . 307437 li  05 

6  « 7805/96*01 
R. 49480206* 15 
3.2  1949006  04 

1  .  I4  75933E  03 

1  4:349/8E»00 
5  6499'41,96-OS 

7  25006 196*01 

1  00000006  *OL« 

2  32149896-07 

7.  5G46091E  *05 

U  ■  00000006  *00 

1  03756446 -OJ 

6  15118436*01 

1  4  11.3  1446*00 

1  3G7437 1 C -05 

14 

1.13157291  01 

1  .BPrr;  o74r  •  1 1 

1  329834  IE -02 

3  25295556  01 
?  604 17821  *oo 
5  ai91493r-03 

2 . 7 149BB5C  *07 
-1  97I559SE  01 

1 . *800804 E  05 

6  80 1 94  4  3E  *OJ 

7  76  1822  if  •  15 

3  2I94900E  0-1 

1  (19  111896  01 

1  5?96  1  1Hk*CKJ 
6.78 19904C  05 

7  77281666*01 

1  OOOOCXK1F  ♦  OO 
7  3G20174T -02 

2  rt  1 65574  6  (05 

O  .'VloOoOOt  *00 

1  054076 76  03 

Qr. /Ii8  13t  *01 

1  • I9S28C6 *00 

1  19008646  05 

'  15 

1  36313256  Ol 

1  HH3380G* • 1 1 

1 .3259051 E  <>.’ 

3.94I4I37E-01 

2  58780126*00 

4  79b/3Gqt  -03 

2  47052776*07 
-2  6573033601 

9  95nr.48JE  OG 

G. 006*8496 *03 

c  f*2  M09i •  ic 

J  71949001  04 

1 . 03256566 -<’3 
’  75  7364 3E *00 
7  1001114E-OS 

2  79563.136  *0  1 

1  OOOOOOOF  *00 
2 . 40286308 -03 

3 . 0584B3BE  *05 
0.00000006*00 

1 . 07873 1 26 • 03 

5. 5G8935U6  *0* 
2.03966456*00 
9.95164836  OG 

It 

1  55G9  19LE  OI 

1  8B45574I ♦ 1 1 

1  32ISB7QE  02 

3.87037196 -OI 
?  570T50BE*(X1 
3.972G527F-03 

2  23 1G I32E  *07 
-3.79 165YMF  01 

8  I7I2ISGF  OG 

6  39093911*03 
4. 793808  IE*  15 
3. 2194900F -04 

9.  73726.15E  OI 
1.92690G6E»00 
8. 1B7IP89F  OS 

2  31851516*01 

1  00000006*00 
7  4422J94F-02 

3. 70328&3F.  *05 

6  oooooo:ie*oo 

1  10006 /2E  03 

5 . 29025 17E*0* 

2. 376066 IE *00 

8  17I2156F  06 

17 

1  7G49HI4F  01 

1 . 8852  I32FM  t 

1  3IG33T4P  02 

4  0*950106  -0| 

7 . 5537584  E *00 

1  ;00087?r  03 

2  00204636*07 

5  38699436  01 

6  5171897E  06 

6  15/24706*03 
3. 592074BF.  ♦  f  5 
3.21949nor-04 

9. 1550675k  04 
2.09862696*00 
9  64473806  05 

7  14  J99HOI  *01 

1  00000001  00 
2 . 47969076  02 

3  19631316*05 
0.00000006*00 

1  ,  17507676  03 

5 . 02  2BC  7  2  E  ♦  0 1 

2 . 72D7663E  *00 

G  .  5  1?  1H97F  -OFi 

—  - - 

— 

r>RA«;  -  t  Mtiao?c»o7  dynes 


e.  Data  extrapolated  to  boundary-layer  edge 
Figure  F-2.  Continued. 
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AEDC-TR-88-31 


CUN IF LOW 

RUN  NO  ENVELOP 
09/02/88 

SPHERE  CONE  NON  EOUlllBRIUM  WITH  STANDARD  PATES 


CONTINUATION  OF  RUN  ENVELOP  f  ROM  RECORD  '  65  SEGMENT  1 

NPCORD  «  70  FLOP  •  12  MINRAv*  17  v  0- 54066E  010  51623E  030. 58648E  0?  LAST  STEP  -O65I0OE-O1  PSlIMMlNl- _ 

NRCORD ■  75V<3J-0  592 I6E -0  IDE LY I  3  I  • 

~  ITERATION  NO  "5  SHOCK  RADIUS  •  0  71321  SIP  •  0.13775  Si  -  0  00786 

RMS  DEVIATION  -  I  8557528E -04  MAX  DEV  ■  1.0975208F  02  ZSTAG  •  4  4325342E02 


PROGRAM  DID  NOT  CONVERGE 


VALUES  AT  INTERSECTION  OF  RAV  WITH  BDOV 


RAV 

ZB 

PB 

PHIB 

DELTA 

FLUX 

5 

4 

5607 2E  02 

4  77592E  02 

1  51713EKXD 

2.B2715E-03 

5. 1524  IE  04 

6 

4 

944806 -02 

9  53t 131  02 

1  4634  ?E  *00 

2  8645 JE  03 

2  05542E  03 

7 

5 

5B339E -02 

1  .  4269  IE -01 

1.409846*00 

2 . 903O0E -03 

4  60199E-03 

q 

6 

4742  IF  -02 

i  895T6E  Ol 

1  356236  *00 

2  95634 E  0.1 

8  133731 -03 

9 

7 

6  14  16E-02 

2.35B87E-01 

1 . 30266E  «00 

3.02863E  03 

1.2607  IE  02 

»0 

8 

999456 -02 

7  8I490E-O1 

1  249126*00 

3.  *  2332  6  03 

1  797S5E-0S 

1 1 

1 

06  2 56 E -Ol 

3.2625BE-01 

1  1 9562E  * 00 

3.243266  03 

2  4 1847E  02 

-  ■  ‘ij 

1 

24876E -Ol 

3  70O65E-01 

1  142 15E  *00 

3. 39549E  03 

3  11659E  03 

13 

1 

45800E -Ol 

4  12792E-01 

1  08B7 IE  *00 

3.5BU52E-03 

3  8844  IE  02 

14 

1 

689f,6r  Ol 

4  54325E01 

1 .03529E*00 

3.802876-03 

4.?1354{  6} 

15 

1 

9431  IF -Ol 

4  94552E-01 

9.8  I878E  01 

4.O64HBE03 

5  59476E  02 

16 

2 

217656 -01 

5  33366E  01 

9.284T5C-01 

* . 36705 E -03 

6  5  1 8  1 0E  -  02 

17 

2 

51253F-01 

5  70*626 -01 

8  75067E  01 

4  7U052E-C3 

7  47275E02 

18 

2 

826996 -01 

6.06339E  01 

8.216466  01 

5  1I355E-03 

8  4 48  2  76 -02 

f.  First  page  of  output  from  restart  of  ‘ENVELOPE* 
Figure  F-2.  Continued. 


1  2 77 380 1C  02 

3  45G9?b'.1f  03 

5  74  141.836  06 

3 

21949006  04 

1  4f.6  1433E-04 

3 

50738G3F  02 

1  .41501676  03 

5  74146  06 

IS* 

1  5294762E-01 

4  370G  1  TOE  Ol 

3  667  If.  19L  105 

4 

H138«3IE*01 

2  30652636*07 

1 

0776387E  03 

6.  • 32 I534E  *03 

3  .  fl6  149961  *0l 

0  OUOOOOOE  *00 

2  3821)7  1GE  *01 

2. 22448096  MDO 

l 

30234B4E* 12 

3  4271009E* 15 

1. 

B700G47E  *  12 

5  B444404E  *0 1 

9  922407  1C  0? 

1  26464601 *03 

2  2462753E  03 

5. 278334  1C -06 

1 

2 194900F -04 

1  61  30/696  04 

2 

5080? 14F  0? 

1  513A992E  03 

5  278334  IF -06 

14 

1  5O9'i067t  -Oi 

4  3R94?0'»|  oj 

3  65 •6  260F  *Ot 

4 

485Q7G4F  ♦(>  1 

2  4651687F**!/ 

1 

00901816  03 

6  :90?669i*03 

4  1679966F  *0 1 

n  OOOOOOOE*OD 

2  4479 12?f  *0 1 

2  3712289L-00 

2 

61987011*12 

2.  1569527E*  15 

1 

89359R3E • 12 

3  6126?07l*01 

8  3336497E  Ol 

1  1 285557T -02 

I.3527P57F  03 

3  34BB716E  U6 

3 

2 194900F -04 

8  9935179E  04 

2 

55845406  07 

i  40024B6E-03 

3  34887 16E -06 

15 

i  6i wy >?  ;»  •  6  0‘ 

4  395653  IE  01 

4  0 1G 1 767E  *05 

4 

0599772E  *01 

2  19574401  *1)7 

B 

BB27734F  04 

P  6085062E *03 

4  3533444F  *01 

0. OOOOOOOE  »00 

2  896 1 4 R3F  *0 1 

2 . 26  1 30496  *00 

0 

OOOOOOOE  •  00 

0  00000006*00 

0 

OOOOOOOF  *00 

o.  oooooooe  *oo 

1  00000006  *00 

o. oooooooe *oo 

O  OOOOOOOF  ♦Of* 

0  (MXKKICXiE  *00 

3 

2  194900C  -04 

7  74392906  03 

3 

096?74<l|  02 

o.oooooooe»oo 

6  OOOOOOOF ♦ 00 

ELAPSED  time  •  2  09  SIC.  ORAG  *2  33BE*Q7  OTMFS  _  _ 

ITERATION  IN  SUPRI *.> I  I NF  F  f  Lb  IS  NOT  CONVERGING 

27  O  6099033  02 178657  0  6163378  O  2103410  0.6195616  0  2205489  O  8119507  O. 5387694 


NRCORU  •  6 

27  O  6168989  0.2236159  O  G234822  O  2*71849  O  6254600  O  2259747  0.8222203  0  5483861 


NRCORD  *  7 

ERROR  AT  STATEMENT  11  OF  SI  IBROU  T  INFrFFCNT 


g.  Last  page  of  output  from  restart  of  ‘ENVELOPE* 
Figure  F-2.  Continued. 


97 


AEDC-TR  88  31 


PAV 

105 


RUN 

FNVEl 


SIGMA 

4.337  169  I E ♦O 1  b 
SURFACF  DISTANCE 


PI 

168BO36E-01 
O  02  1  I  FEET 

f « 


DEL  IA 
5  58351b2E -03 
ENTROPY 
GWALL 


2  BOUV 
3. 2901  7021  01 


SBAK 

2  6GG789BI -02 
O.  25'j227E  *Ot  CAL/GM/K 

FRICT  OfKH 


R  BOUV 

5  3769437E-0I 


PHI  BOOT 
3. 2096135E*O1 


0  OOOOOOOE  *00  0  OOOOOOOE  *00  1  1 894U32F -03  I  767G375E  OT 


NPT  Z  (CM) 

RAO  ERG/GM  SEC 
MU» /CM  U 


R  I CM  I  U  (CM/SEC  I 

•M«  (GM.'MULE  I  MACH  NO 
W  E 


PHI  IOEO)  P  I OYNF /CM*  *  2  I 

PLASMA  IWEO  EL  UENSL*/CC| 
AR  02 


RHO  IGM/CCI 
COLL  »REli 

M2 


T  (DEG  K )  _  PSI_(GM/CM«»2SCC) 

II  CONO  AV  CHI 

NO  _  NO*  _ 


5.71  I2227E  *  03  3  946IB22E*00 


*  260502 7E  -Ot 
O  00000001*00 
1  30057 33E -02 

3  23038361-61 
0. OOOOOOOE ♦ 00 

•  300443  IE  02 


5  4242440E-O1 
2  3797  T42F  *0 1 

1  920792  IE  03 

5.47226CU  O* 

2  38243321 *0 * 
I  B7830G9E  03 


3  840  1 0631  *  09 
2  40976B0E *  OO 
3.91 32402k  06 

3.R199156E*05 
2  4053180C*IX> 
3.9O77817E-0G 


3  2890M44F10I 

2  05254GOE  ♦  1 2 

3  2 194 BOOT  04 

3.  35878 181*01 
2.  10789  IOC  ♦  12 
3  21949CXX  04 


I  I204538E  *07 
I  3239357E ♦ 1b 
1  27  * BBHbE  04 

1  IB694C5E  *07 
1 . 396295  TE  *  10 
1  33G6756E  04 


5 . G 153000)  04 
9  28279B4L ♦ 1 1 
2  54  134701  02 

5  93048301  C'4 
9  8  I67767L*  » I 
2.5415541E  02 


4 . 5484 2 7 3E  *OI 
1 . 173B338E  03 

5  73501481 *03 


4  4562B20E  01 
3.9132402E  06 

8  24736b7|*00 


4.53609111*01  3  7860747E  01 

1 . 2  17  10 J9E  03  3  9077817E  06 


5  3 .*19075001  01  5.5254570E-OI  3  8447519E*05 

0 . 0000<X)0€  ♦  00  2  386  1  B25E*01  2 . 4059738E  *00 

l  29370596  02  1 . 8 13B5B5E *03  3  867  * 266E -06 


3  44O2660E  *-0 1 

2  16B399 IE •  12 

3  2*94900F  04 


1  27 16890E  *07 
1  4776090E*15 
1  42834R1E  04 


6  34  1B22SE  04 

1  05000  7  2  E  *  12 

2  5423246E  02 


5.75497191*03 
4 . 4R7530GE  *0J_ 
1  . 26  i  2676E  03 


1  .  4  3268  8  51  •  CJ  1 
-3.  1632  1  1 T- F  0  i 
3  867  1266E -06 


7 


3  1 559268 fc  01 
O. OOOOOOOE *00 
1 .207053  IE  02 


5  5909 78 2C  Ol 
2.3B94F50t*01 
1  756639  if -03 


3.85U9592F  *05 
2 . 4099050E *00 
3 . 9  187343F  life 


J  5049811  I r*oi 

2  2 15482  It  *  12 

3  2 194ROOE -04 


I  345295 U *07 
I  54 2473 1 E *  15 
I  5  1  IP4i»F  04 


8  704 1 094 E -04 

1  1096946E • I 2 

2  543094  IF  02 


5. 7670472E  *03 
4  4328946E*0i 
1 . 3031402* -03 


2  00715B7E*0I 
-2  74936761-01 

3  0  I8  /34  if  -06 


9  31 193794k  01 

O. OOOOOOOE «00 
1 . 28 1 7 1 46E  02 


5  6492485E  01 
2 . 39297R0E  *01 
I  6959733F  03 


3  959  16  14E  *05 
2.41 39972E  *00 
3.7570354E  06 


3  5677795E  *01 

2  259  1 8Gr>l  i  12 

3  21949COE-04 


l  4jh?774C407 
1  6039295E *  15 
1  60323B9E  04 


7  0B57O28E-O4 

1  1  /  I8245E  *12 

2  543972  IE  02 


5  775  I OBOE  *03 
4  365I17BE»01 
1  3463068E  03 


2  65  I B 1 37E  *0 1 

2  410O270E0I 

3  7570354E  0€ 


11  3 . 083 101 7C -01  5. 7070RB8F -Ol  3  8688G?7E*G6 

O  OOOOOOOE •  Oil  2.39648851*01  2 . 4209 1 47E *0O 

1  2755649E-02  I.6351820C  03  3  6861751E  06 


3 .  829B096E  H)i 
2.2970B3BE* 12 
3  21949001  04 


1 .4969209C*07 
1  b5H 1 9  I  If •  15 
1  099601  If  04 


7  4662372E04  5  7786997E*03  3 . 335 1788E *01 

1  23308521*1?  4  288594 3E *0 1  2 ■  14B5848E  01 

2  5448904 E  02  }  3B8ft0'28E  03  3  686175 1E08 


13 


3  05O9255E  01 
O  OOOOOOOE *UO 
1  27O0O9OE  02 


5.7583B97EOI 
2 . 3996790E  *01 
1  579927  IE  03 


7  8790633E  *05 

2  42886VOE *00 

3  613824GE  06 


3 , 67B0499C  *0 1 

2  3252370E  *  12 

3  2  194900E  04 


1 . 5623554E  *07 
I  6990867E ♦ 1 5 
1  7BB3435F  04 


7  803535BF  04 

1  2B66635E ♦  1 2 

2  5457727E02 


5  7784893F *03 
4  21 1 3742E*0 I 
1  . 4264850E -03 


3  978477 | E  *0  I 
I  97  I  I 298E -O 1 
3.61 38246E -06 


15  3.0249355E  01  5  793B275E  Ol 

O.OOOOOGOE*00  2.446671CE*01 
1  194  769  IE  02  7 . 36644 19E  04 


4  1 62 2524 1*00 
2. 71 96464 E *00 
1  9Q247 37C  06 


5  7  218 1 3R(  •  >  J  * 

1  795954  7  r « 12 

3  2 I94900F  04 


1  62201111*07 
I  0I36123E*  15 
7  6335909E  l>4 


8  8429295F-04 

1  3624797E *  1 2 

2  5687694F  02 


5  3975565E  *03 
2  37 25485E  *0 1 
I  9  1 154AOE-Q3 


4  572 i5HOE  *0 1 
-1  85 I00B2E -O 1 
1  90247371  06 


’7  3  0C92924F  01 

O. OOOOOOOE *00 
1  1 9  ~ 09  T  4  E  02 


5.82476841  01 
2 . 4  449505E  *01 
7.68904 15E  04 


4.  143790 1C  *05 
0  6984 3G4E  *00 
1  9B21789E  06 


3 . 7443822E  »0* 
1  85 104 90E  *12 
3  2 1949O0E  *04 


I  665209 1E*07 
I  076755  IT* 15 
3  53779871-04 


9  OI8O650E04  &  4J11956£*03  4  9692013E*0I 

1  396007  IE *12  ?  4598 157E *01 _ *  1  6662 38 BE -O I 

2.5673626E  02  1  B073424E-03  \ . 912 1769E  06 


19 


2  99650351  01 
0  OOOOOOOE *00 
1  1 873349E -02 


5. 84*503  1 3E -01 
2.4493737E*01 
.  19180351  -04 


4  i606127E*0b 
2.7i 70068E  *00 
1  8797883E  06 


3 . 7590353E  *01 
1  0?8OO12E_*<2 
3 . 2 194900E  04 


1 .702617 IE*07 
1  05010B2E* 15 
"3. 762566  If  04 


9  271885TF-04 
1  429H49E*  12 

?  S6T4|c4f  02 


5  409 7 44  IE *03 
2  3433663E401 
1  856 JT5TE -03 


5.30066751*01 
1 . 5Q72022E  -Q1 
i  87978B3f-06 


21  2  9B52156E  01  5.86315591-01 

0  OOOOOOOE* OO  2  456G702E  *0 1 
?  17127051-02  6  1816102F  04 


4  I937351F105 
2  75482091  *00 
I  7025933k -06 


3  77 19559E  *0 1 
I  76778291*'. 
3  2194900E-04 


I  736929  IE  *07 
9.82062 I7E*  14 
4  2225B13E-04 


9.5735062E  04 
1.4G2150BE* 12 
2  567809  I E  02 


5  3608 1 32E *03 
2  1 42003 16 *0i 
1  92943391  03 


5  6089824E *01 
1 .4Q51299E  01 
1  7025933E  06 


23  2  9748485E  01 

o  ooooooor*oo 


5  880007  IE -01 
2.4G4397  3F  *01 


4  .?297876r  *05 
?  7P7H74MOO 


3.7B39869E*01 
I  r,0Al379F*  12 


1  7694724E  *07 
o  0620*.  191  *  14 


9  8A6BB6 If  04 
»  494547BT *  1? 


5  304784 1 E  403  5 . 908268 7f *01 

I  933706  IF  *01  -  I . 34 I7890C -Ol 


h.  Typical  output  from  MOC  calculation  run  number  ‘ENVELOPE’ 

Figure  F-2.  Continued. 
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25 


1  .  15368761-02 

2~9«46288E ^01 
0  QOQOOOOE  >00 
1  1361493E-02 


5  58592391  04 

5 . B959790E :0l 

2.4721550E»Q1_ 

4.7968947E-04 


1 .5212945C-06 

4  26640 15E *05 
2  8406055E*00 
1 . 3298G43E  *06 


3  2194900E -04 

3  794857  IE *01 
1  612882  IE* 12 
3  7194900C  04 


4  7528197E  04 

F  80086785 *07 
6  174973  IE *14 

5  4518468E  04 


2  5683235E -02 

1  020286  IE -03 

1  5256746E ♦ 1 2 

2  57048971-02 


2? - 2  “9550922? ^01  5 . 91 1 1839E r6l  4  30556565*05  3  8O22503E*0l 

0  OOOOOOOE *00  2 . 4821791E-*01_  2  B835O9O£*00  i  4540553E*i2 

'  ~  1 . 1 134B9BE-02  3. 8012794E:04  I050B619E-06  32194900E04 


1  8  306503E  *07  I  04939BOE  03 
6  6442037E*  14  1 . 5508502 E M2 

6  9980695E  04  2  5796003E -02 


28  2 . 94704B0E -01  5.9;n  R58I  I  «  36283! 

_ O.OOOOOOOE*OQ  2. 5Q67Q27E *01  2  9326828E*00 

I  04B7464E02  2 . 4023323E -04  5  77597 lOE -07 


3  8034583E*Oi  1  B503962MO7  1  O677164E03 

1 .OB73727E*  12  3  7I56773EM4  1  55223B1EM2 

3  2 194900E -04  1  199926 «E -03  2  60423521  02 


29  2  9400224E-01  5  935210BE-OI 

O  OOOOOOOE *00  2.J493fl96E*OI 

9.24929671-03  1.43181611*04 


4  4053Q32E  *05 
2  94 12299E  *00 
2  1904948E  07 


3. 768508 '1E*01 
6 . 6586 147E* 1 1 
3  2 194900E -04 


1  8G05337E  *07 
I  . 39331361*14 
2.  12321  16E  03 


1  055722  IE  03 
1 .5I929J4EM9 
2 . 6395080E  02 


30 


2  9333463E  01  5  9458550E01  4 . 439 I53UE *05  3 . 7 14B03GE -O 1 

0.00000001*00  2  655»979E*gj!  2  9O7376QEK>0  2.4352428EM1 

fi  iMMWMOST  0998179?  05  3  14612ME  OB  3.2194900E  04 


1 .6369145E«07  98318023104 

1  8636555E  *  13  1.4  103916EM2 

3  99178821-03  26772928E  0? 


31 - Fe 2906171 -Of  5.9526B63E-01  4  46715535*05  3 . 6  1  189B0E  *01  '  7942228F*07 

0  OOOOOOOE *0O  2  09614B3E*O1  2  75430981*00  O  OOOOOOOE *00  0  00000COE»00 

- O  OOOOOOOE Foo  o  bobddOOEKW  O. OOOOOOOE *00  3.2194900E  04  7  24392 90E -03 


8  7329982?  04 
0.0000000F*00 
2  B962740C02 


ELAPSED  TIME  »  2.7.2?  SEC  DOAG  »2.752E*07  OVHES  _ 


NO CORO  -  106 

NR CORO  -  107 

NRCORO  108 


NR CORO  -  109 


h.  Concluded 
Figure  F-2.  Concluded. 


1  99BR969E -03  1.5212845E  06 

5  248 175 IE *03  6 . 2044053? *0? 

I  7Q88254E *01  -1 .0335748E-Q1 
2 ■ 0346662? -03  1  3298843?  06 

5 . 20B0005E  *03  6  49815555*61 

1  3663009E *Q1  _8  6648403E  02 

1  9522949i:03  1  0508619E06 

5 . 22497681*03  67755? 168E*01~ 

7  634Q  I  3  I E  *00  5  30908 96 E  -Q1  _ 

‘l  5999646E-03  5.77597105  07 

5  40377061  *03  6~.  9B3T428E*OT 

2  924698 2 MOO _ 8  .  405 1B04E  *01 

9  91919111  *04  2  19O494iE-0‘> 

"5  9665797E  *03  7  .  F935B25E*01 

4 . 2 14Q452E  ~01  9.723Q34BE  01 

3  0859389E - 04  3. 14812 1  IE  OB 

7  15659695*03  ? 79636562E*0 1 

O  OOOOOOOE *O0  1  OOOOOOOE *00 

6  oooooooe*o6  O. OOOOOOOE *oo 
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Example  2.  The  second  example  to  be  considered  (Fig.  F-3)  illustrates  input  for  a  calculation 
of  equilibrium  air  flow  over  a  flared  body.  The  free-stream  conditions  are  obtained  from 
fits  to  the  1953  Standard  Atmosphere.  Specifying  cgs  units  explicitly  was  unnecessary. 

Figure  F-4a  is  a  reproduction  of  the  first  page  of  the  output  from  run  EQUILI.  Notice 
that  the  altitude  and  velocity  are  printed  out  in  English  units  and  that  the  thermochemical 
properties  are  not  printed.  The  results  of  the  third  marching  pass  are  shown  in  Fig.  F-4b. 
There  were  no  integration  failures  (FLOP  =  1)  and  the  RMS  deviation  between  specified 
and  computed  body  points  is  certainly  acceptable,  but  the  relative  error  between  calculated 
flow  angle  and  the  angle  of  the  body  surface  at  the  junction  ray  is  approximately  1  percent, 
so  the  message  “PROGRAM  DID  NOT  CONVERGE”  is  printed.  DELTA  and  FLUX  are 
zero  because  the  boundary-layer  option  was  not  selected  (INBLR  =  0).  Printout  from  the 
thirteenth  data  surface  is  shown  in  Fig.  F-4c.  AVER. CHI  =  0.0  at  all  points  as  it  should 
for  this  equilibrium  calculation.  No  composition  data  are  obtained  when  MODEL  =  2,  so 
the  electron  density  is  reported  as  zero. 

The  M.O.C.  printout  at  the  flare  discontinuity  is  reproduced  in  Figs.  F-4d  -  F-4g.  The 
message  at  the  bottom  of  Fig.  F-4d  indicates  that  the  next  ray  would  ordinarily  be  located 
at  the  dimensionless  station  ZRAY  =  1 .6772  but  since  this  station  is  downstream  of  the 
segment  end  ZCHECK  =  1.626,  Ray  59  is  placed  at  the  junction  point  and  printed.  Now 
the  code  discovers  that  this  junction  is  a  discontinuity,  so  it  calculates  Ray  60  normal  to  the 
cylinder  (body  segment  2)  at  the  junction  point  (Fig.  F-4e)  and  then  calls  subroutine  FLARE 
to  determine  the  initial  inclination  of  the  imbedded  shock,  Fig.  F-4f.  The  match  error  is  that 
between  the  flow  angle  calculated  from  momentum  and  energy  conservation  and  the  specified 
body  angle.  The  flow  conditions  just  downstream  and  just  upstream  of  the  shock  at  the  body 
surface  are  printed  as  Ray  61 ,  Fig.  F-4g.  Note  that  the  point  identifier  and  the  point  coordinates 
are  repeated.  Figure  F-4h  reproduces  the  first  ray  downstream  of  the  discontinuity. 
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Es»ip'l  input  fop  fuuh  rPRti.w  mb  now  over  sphere -cvunufb  coni 

tOP  T I ONS 

MOnri  .2  :SPECIFICS  EQUILIBRIUM  AIR 

lUNIT-O  ;  r  PE  E  STREAM  CONDI  I  ION'S  1  RDM  *0*3  ST ANOARO  ATMOSPHERE 

INBLR-o!  :  NO  VtSCIM*  IN  II  RAl'T  ION 

INAFT-I  .  INCLUDE  M  0  •  CALCULATIONS 

UNI T" #CG5 #  ;«E?RI»  UNI  IS  rOR  INPUT 


tURSTREAM 
VEL-4  876BE5. 
ALT "2 . 7432E6 . 

% 

tSHOCK 

RSHP-2.  16135. 
SIP-0  09067. 
SI -O  0061. 

% 

sboov 
REXP-2  . 

CIO  n. 

C2-  -  I 

MORE-  T..  : 

i 

SBOOV 

RJUNC-1  905. 
ZJUNC-1  905. 
REXP-1  . 

C1"0  00. 

MORE ■ . T . . 

f 

SBOOV 

ZJUNC-3  50. 
RJUNC-I  905 
REAP"  I  . . 
2EAP-1 . . 
C0-0.07. 

MORE  -  F  . 

% 

SCONTROL 

NBSIT-3. 

1 PRNT "5 .  i 
OELK"  05. 

AO"  05. 
METM00-10 
0ELN1N" I  OF  * 6 
PSIRAT-0  «0. 


;tENTERLINF  VEIOCITV  CM/SEC 
;  ALTITUDE  IN  CM 


;CM  RADIUS  01  CURVATURE  OF  SHOCK  AT  A»|S 


; EXPUNEN1  OF  |R  R JUNC ) 

CM  1MICF  T I  IF  NOSE  RADIUS 
;  SINCE  IMIS  IS  A  SPHERE 
.ANOTHER  p.ouv  stgment  r 01  lows 


CM  RAD1AI  COORDINATE  SPHERE  CONE  JUNCTION 
CM  AA1AI  COORD  I NA I E  SPHERE-CONE  JUNCTION 


TANGENT  0T  CONE  HALF  ANGLE 


.BEGINNING  Of  FI  ARF 


iC-PONfNI  OT  17  *  r.MJNCl 
.  TANGENT  Or  FI  ARE  ANGLE 
.THIS  IS  LASI  nODV  SFGMENT 

-NOT  MORE  THAN  J  SHOCK  SHAPE  ITERATIONS 
STATUS  RF POO  I  TVER*  5TM  STEP  FROM  SHOCK 
; AVERAGE  RAV  SPACING.  DEFAULT  IS  04 
:  POSITION  or  11051  thin  axis  RAv 
PPOBAOL  »  IINT4FC  ESSARY 

i  M'NIMUM  SirrSI7E--INITIAL  STEP  “  I  .  OC  -  *1  •  RSHP 
DEFAULT  MAP*ll|N<;  STOPS  WITH  107  Of  f  l  UX  LEFT 


PS  I  RAT  "0  TO ,  .  lit 1  AUL  i  mipt.  a.u.j  .  . 

EBBOB'I  OE-S.  CUNVEBOlNlE  CBITEBION  EOB  MUSI  UFBM10NS 
USE  15  POINTS  TO  START  M.O.C.  CAICIIIAIION 


NOUT- T5. 

t 

% SUPER 

ZMAX-6  25.  I  BOP  •  IS  «  75  M  I  ON' l 

2E  T  AM»  -O  0.  ;M  O  r  RAV  SPACING  IS  O.B'C  F  I  Al  I  OMANI  F 
NMAX-35.  ;MAXimiM  NUMBER  OT  STRFAMUNFS  IN  Mf  SH 
MOORAV-tO  : PRINT  EVERY  HUH  RAV 

NCORN12I-2.  ;  SIGNALS  THAT  11115  CORNER  IS  A  COMPRESSION 


Figure  F-3.  Sample  input  for  equilibrium  airflow  over  flared  body. 
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CAL  BLUNT  BODY  PROGRAM 


PUN_NC  EQUILI  _ 

09/08/68 

SAMPLE  INPUT  KOH  EQUILIBRIUM  A  IB  FLOP  OVER  A  SPHERE -CVl INUER  CONE 


ALTITUDE  •  9. OOOOOOOE *04  _T  T  VELOCITY  *  I  60000001*04  FT/SEC 


“PRESSURE  OYNE/SO/ CM  TEMPERATURE  DEG.K  PENS  I  TV  GM/CU. CM  VILOCtTr  LM/5EC  MOLECULAR  WEIGHT 

I -7609575E*04  7  23964111*02  27388232105  _ 4  876BOOOE+OS  2  89614831*01 


INIliAL  SHOCK  PARAME1FR5  -  RSHP  •  1  1613  SIP  ■9.B6700C  02  Si  ■«  IOOOOE  03 

PEEK  -0.0500  XO.  *5^0500 


_ IBDY  YBDV  801 _  B02  _ BDC(l)  _BOC.(  »!  ...  B0C(3)_  B0C(4J _ BOyS) _ 

o  oooooooe *oo  o  ooooooor >oo  2  00000001*00  o  oooooooE«oq  0.00000001*00  3  8 ioooooe *00" 1  0000000000  0  oooooooE»oo  o  ggoggooE*gg _ 

T  9050000EiOO  1  90S10LX**00  i  OOOOOOOE *00  0  00000001*00  O  OOOOOOOE  *00  0.00000001*00  6  00000001*00  0  00000001*00  0  00000001*00 

3  50000001*00  I  90000001  *  DO  I  00000001*00  I  OOOOOOOE *00  2  00000001-02  O. OOOOOOOE *00  0  OOOOOOOE  *00  0.00000001*00  0  00000001*00  _ 

EM1SSIVJTY  •  0  00000001*00  INITIAL  MALL  TEMP  •  2  2222222E *02010. K _  _  _ 


I  UNIT  1EJBLB  1  NOVI  I  NAT  T  INOIM  I PRN 1  NRDY  MOD  F  l  NSPEC  NR  NC  _  NVIB  NBAO HKOW  M1TH00 — 

- o - “o  6  I  '1  's'  3  2  8  0  O  O  O  20  10 

- slaS -  skap  errob“  rmstst  tesim  msf  gamtsT  deImJeT  “ 

-S  OOOOOOOE -01  3.5000000E-0I  I  .OOOOOOOE  “05  5  00000001-04  I.  10000001*00  I  500000QE0I  I  00000001^05  1  oooooooe_-o6 


ERIE  STREAM  MACH  NO  -16  2088 

KLMROA  *  3  69900571*02  HINT  *  I  DB4T597E»02  BIGAM  *  I  4079C86E *00 


SHOCK  GEUME1RY 


M 

XCM» 

S2IM) 

5B(M) 

SIGMA 

curvimi 

0.  12496E  02 

0  49979E  01 

1.5208 

0  99920 

- - 

2 

6. loooo 

0  49931E-02 

O  99B34E  01 

1.4709 

-0  9966B 

3 

O  15000 

0  1 1215E  01 

0  14944 

l  4212 

0  99322 

4 

0  20000 

0  I98B9EOI 

0  19868 

1.3717 

-0  98548 

0  25000 

(1  30977E-O1 

O  24743 

1.3227 

-0  97607 

— 

6 

0  30000 

0.44429E  01 

0  29558 

i  .2742 

0  96367 

7 

0  35000 

6  40000 

0  60IB5E-01 

0  34303 

•  .  22C4 

0  94798 

— 

ft 

0  78  170E-01 

O  38967 

t  I7Q4 

-6  92886 

9 

0  45000 

0  98301E-01 

0.43544 

1  1335 

1  OftBB 

•0  90614 

10 

0  50000 

O  1 70*18 

0  40024 

-0  88057 

1  1 

0  55000 

O  14462 

0.52403 

1  0455 

-0  85109 

12 

0  60000 

0  17059 

U  56675 

1 .003; 

-6  82075 

13 

0  65000 

O  19830 

0.60837 

0  9C348 

0  7B76B 

1A 

0  70000 

O  73762 

<)  64886 

0  92495 

-0.75376 

15 

0  75000 

0.35B45 

0.68822 

0  86816 

-071 806 

- - - 

16 

6  80000 

0  29067 

0.72645 

6  85315 

-0  68262 

17 

0  85000 

0.324  17 

O. 76356 

0  ft  1990 

-0.64740 

a.  First  page  of  output  from  run  ‘EQUILI* 
Figure  F-4.  Output  from  Example  2. 


ITERATION  NO  3  SHOCK  RADIUS  •  2  18276  SIP  *  0.  10333  SI  -  6  110639 

RMS  DEVIATION  ■  2  1 1866781  -04  MAA  DEV  *  I  . 07932B9E  02  /STAG  =  4.7660I6IE  02 


PROGRAM 

DID  NOI 

CONVERGE 

VAIUFS  AT  INTERSECTION  Of  RAT 

WITH  BOOT 

RAT 

?B 

PB 

PH  IB 

flEtTA 

FLUX 

5 

4  895B5E  02 

4  759231  02 

1  516?5E*00 

0  'J0000E*00 

0  OOOOOE *00 

r- 

5  ?8 40 BE -02 

9  50364 E -02 

1  46 1 7  IE  +  00 

o  oooooe  «oo 

0  OOOOOE *00 

7 

5.93I76E  02 

1  42IBSE-01 

1  407I9E+00 

0  OOOOOE « CO 

0  OOOOOE  *00 

8 

6  BJ422E-02 

1  H8B95E  01 

1  35268000 

o  ooooor  *co 

0  OOOOOE *00 

9 

7.PB922E02 

2  .1502 1C  01 

»  29831E+00 

o  ooooor*oo 

O  000006*00 

lO 

9  79286E  02 

2  B0427E01  " 

1  24376E lOO 

0. OOOOOF  *0O 

O  OOOOOE *00 

1  1 

1  10406E01 

3  24977C  01 

1  10934F  +00 

0 . OOOOOE  *oo 

0  OOOOOE  *oo 

12 

1  29770E  01 

j  4B542E  Ol 

1  134«UF*00 

0.000006*00 

O  OOOOOE 'OO 

13 

1  50464E  01 

4  I0994E01 

1  08059E  +00 

0  'JO000fc*00 

O  OOOOOE *00 

14 

1  7392  IE  01 

4.53?  I2E  01 

1  03625E*00 

O  OOOOOF  *ob 

0  OOOOOE 10 O 

15 

1  9957  IE -01 

4  9208  IE-01 

9  7 1938E-0I 

o  OOOO0E«OO 

0  OOOOOE *00 

16 

2  77340E  01 

5  304RBE-01 

9. 17639E -01 

O  OOOOOF  *O0 

0. OOOOOE *oo 

17 

3  57150E  01 

5.6737BE01 

8  63345E-0! 

O . OOOOOE  *00 

0  OOOOOE  *QO 

1ft 

2  flfl9  1 9F  *01 

6.02500E  Oi 

B.09045E  01 

0  OOOOOE •«> 

0. OOOOOF *0O 

— 

b. 

Results  of  the  third  marching  pass 

Figure  F-4.  Continued. 
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STEP  size 

•  f . 09 1 1 3f  *02 

MUMRF  R 

OF  SUrCESSFULL 

STEP*.  *  13 

D1S1ANCF 

FROM  SHOCK  •  9. 

012731  02  CM 

RAY 

i  Chi 

t?  i  M 

r  D.Nt.cM-.? 

1  PIG  i 

mo  GM/CC 

*.MW  GM-MUl 

V10I  CM/ SIC 

FLOW  AND 

— 

— 

HTOl  IRC/ CM 

S  cai  /*:m  k 

AVER  CHI 

El  DfNS.A/CC 

MACH  ND 

FLUA  RAI 10 

RAD  ERG/CM  SEC 

PSI  GM/51C  _ 

— 

— 

4 

9  OI27331FPJ 

0  l>VX>N>OE  ♦OO 

6  725S604E  *t)0 

'.»  7627539E  *03 

3.U3J9272E  04 

7  J34 10B3F  *0 1 

6  10 1 BOB9E  *03 

O  OOOOOOOE *00 

— 

— 

—  ■ 

t  20276671*11 

4  6?u*»Bir.r*03 

0  OOOOOOOE  *oo 

O  '  XJOOOE  *  OO 

3  9591 136E -07 

1  00000001*00 

0  OOOOOOOE ♦ OO 

0  OOOOOOOE *oo 

— 

5  9  274170lf-n? 

1, '4564071  ui 

6  20 1 00? 4  E  *0-4 

T.  75B4693E ♦  03 

3 .02443691 -04 

2. 33448  7  7E  *01 

2.4514^371*04 

7 . 20796 1 4E  *  6 1 

— 

— 

1  2043612E *11 

4  5707 764 F* 03 

0  OOOOOOOf  *oo 

0  OOOOOOOE *00 

1  5706729E  Ol 

1.3?BI0I2E*00 

0  OOOOOOOE  *oo 

6.6293433E  02 

— 

6 

1  0057523E-0I 

7  OSR^3741-O1 

fi  1279  IB5E*Oti 

5  7469474E  +03 

2. 99627  191  -04 

}  33664 79E *01 

4 . 7  ?G984 1 E  *04 

7 . 593867 16*01 

— 

— 

t  2090573E  *  t • 

4  5t90in9F»03 

0  00000001*00 

O  OOOOOOOE *00 

3.0C4I202E  Ol 

1  38?  1 923E  *00 

0  OOOOOOOE *00 

2.7520583EO1 

— 

7 

1  1 35988  IE -0 1 

3  i2fi7nC7f01 

6.008  1252F  ♦HO 

5  72  /373 IE  *03 

2  95026961  *>4 

i  3J-G  1791*01 

7. 1 12551 1E«04 

7 . 5062886E  *01 

— 

1  2  166365F* 1  1 

4 . 5  tn4<?46C«03 

0. OOOOOOOE *00 

0  'tOOOOOOF  *oo 

4  5710740F-0I 

t  47  131 R3F  *00 

O  00000001*00 

6  5660002E  01 

— 

B 

1  31762986  Ol 

4  Ot 

6  B447700F  *rw* 

5  «9«KM05F*03 

?  Bft/MOIF  04 

7. 34050671 *01 

9  43027761*04 

7.30i9<»t2E*01 

— 

— 

_JK2267584F*  1 1 

4  51621031*03 

0  OOOOOOOF  +00 

0  OOOOOOOE  *00 

6  O 7 75598 E -01 

1  5935927f*DO 

0  OOOOOOOE *00 

1 . 2569909E *00 

-  - - 

9 

1  5499696E-OI 

3  |T*«4-»7E  0\ 

5  8424  1296  >0*j 

5  66116701*03 

2*11 1475E  04 

2  344  J776T  »0 1 

1  1 7 16715E  *06 

7  04  R  999  2  E  ♦  0 1 

— 

— 

1  23BB91  f • i . 

4 , 5 I29049F  *03 

0 . OOOOOOOE *00 

0  OOOOOOOF  *00 

7  579368  IF -01 

1  745B302F*00 

0  OOOOOOOE  *oo 

2  1357776E  *00 

— 

lo- 

1  B320T64E  C\ 

6.  ‘8G7P951.  oi 

5  406BG44F ♦HC 

7  G 1732876  *03 

2  72247GI6  04 

2  3497 127F  *01 

1  39604 39E  *05 

6  775 184 1E*0! 

— 

— 

1  2526628E* 1 1 

4.^0840401*03 

O. OOOOOOOF.  U*0 

o.‘jooooooe*oo 

9  O’  35005E  01 

1  92*46661*00 

0. OOOOOOOE »O0 

3  359816  IF  *00 

— 

n  “ 

2  I627910E -01 

7  IB04//7E  Ol 

5  1450t0Bf06 

5  5544962F  *03 

2  6748975F  04 

7  3554437E  *01 

1  8  ISO  1026  *05 

6  49292626*0 1 

— 

t  2874244E+I1 

4 . 5074 !?4E«03 

O. OOOOOOOE *oo 

0 . OOOOOOOE  *00 

1 . Q557659E  *CO 

2.  I740398C  *00 

0  OOOOOOOE • 00 

A  9943300E*00 

— 

12 

2  54O7265E-0I 

4  16fl»94BF  Ol 

4  .e*>4o->«i^t*o'; 

6  4B  J54  HX  *03 

2  52 1601  IE -04 

i  31. 786391  »OI 

i  8274/796  <  05 

6  20902986*01 

— 

— 

t  2B25I26E*  » 1 

4  49407O«E«O3 

0  OOOOOOOE *oo 

O  OOOOOOOE *00 

1  70300?  3E  *00 

7  3394400E  *00 

0  OOOOOOOE *oo 

7  0**P6598E *00 

13 

2  96427B0F-01 

9  1  186  1  I3E -01 

4  5708885E*O6 

5  400 1064 E  *03 

2  4150370E  04 

2  37152326*01 

2  0324  *686*05 

5  92774636*01 

— 

— 

1  .  2975540E ♦ 1 1 

4  48GOOGGF  *03 

0. OOOOOOOE *00 

O. OOOOOOOE *00 

1  34 B8098E  *0O 

2  56509531  *00 

0  OOOOOOOE *oo 

97 1 89004 E *00 

— 

14 

3  43  I640BF  -0 « 

t  OOr*9  1.*3f  *oo 

4.271184  /E**:B 

r.  304  t  nos E* 03 

2  3069IB2E  04 

2  JSI?R?RE*Ot 

7  27B97B9E*»>5 

5  65 13U14E *01 

— 

— 

1  3  1 23056F ♦ t 1 

4  4736397E  *03 

0. OOOOOOOF *O0 

0. OOOOOOOE *00 

1  4928973E*0n 

2. 7965399E *00 

0. OOOOOOOE* 00 

1  2B80557E  *0 1 

— 

16 

7  94083S2E -0 1 

t  rmsfil'tr  *oo 

3  9G«»7l30l*0e 

5  193  t 12 IE*03 

2  ->OOt37?E  04 

2  3973754F  *01 

7.41075746*05 

5  30O7  22OE  *0 1 

1  3257337E ♦ 1 1 

4 . *15003  .IL  *03 

0.  OOOOOOOE  *oo 

0 . OOOOOOOE *oo 

1  61539G5E  *00 

3  0104 24  YE  *00 

0. OOOOOOOE *00 

1  6629347E*0I 

— 

Is 

4  4B97555E-OI 

1  M7HUMCI00 

3  670693BI *06 

5.06535431*03 

2  0964  706E  04 

2  4048528E  *01 

2  50302996 *05 

5  ueiooSeioT 

— 

— 

1  3372457C*  t 1 

4 . 444O0O3E  *03 

O  .  OOOOOOOE  *'X> 

0  nuO*XXX)f  *O0 

1  .  77723I6F  *0«* 

1  26621541*00 

0  OOOOOOOF* CO 

2  o^r  45-»3r*o» 

1 . 27547 121*00 

3  3785497E  *06 

4.92322931*03 

1  99623 1BE04 

2.41/86801*01 

2.761 4550E*05 

4  85904981*01 

t . 347 1750E* 1 1 

4  42651561*03 

0  OOOOOOOE *O0 

0^  OOOOOOOE  *oo 

1  91958316*00 

3  5035604 £♦ 00 

0  00000001*00 

2.591 197«f+01 

— 

c.  Output  from  the  13th  data  surface 
Figure  F-4.  Continued. 
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2  2664 1 70C *00  2 . 77555771 *00  4 . 2398b26f *06  2  314938  I  1  1485604E+06  1 . 543 1604E -04  2  5596172E+03  2  9028 

1 .21 18949E4  1 1 _ 2  .  8584868E»O  I _ 4^  4  497830F  »00 _  _ 


_ 2 ■ 26641 7QE*OQ  2 . 8 1 73468E+OQ  4  30128 18E«OS  2  3B33239E«01  1  2782196E+06  »  832B30BE  04  2 . 4007» 15E*03  3  3306085 E *02 

1.21 16292E+ 1  I  2  865B443E+0«  46J76H45F+00  “  -  - — - 


ELAPSEO  TIME  » _ 11  90  SEC  DRAG  -3.731E+Q7  DfNCS 


NRCORD 

• 

50 

NRCORD 

m 

51 

NRCOHD 

* 

52 

NRCORD 

• 

53 

MRcQRO 

■ 

54 

NRCORD 

. 

55 

NRCORD 

* 

56 

2 

S 

1 

. 

57 

END  OF  SEGMENT  3  ENCOUNTEPEO  3CHEO(  ._  I«3«33P»»~ - T«f73 


d.  Method-of-Characteristics  (MOC)  printout  at  junction  between  body  segments 

Figure  F-4.  Continued. 
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RAY  RUN 

60  60UILI 


sfoM*  i*i  _  delta  sraq  z  nooi  RTgoov  PHI  BOOY 

7  26088846*01  O  00000001*00  O .  OOOOOOOE  *  OO  O  OOOOOOCE  *00  _  3  87  10  1 786  *O0  I  8B9742«E*00  O  0OOOOOOE*OO 
SURFACE  0151  ANC E  *  0  1492  FEET  ENTROPY*  O  4524886*04  CAL/GM/K 


NPT 

Z  1  CM ) 

H  IERG/GM) 

r  rear 

(OM/MOEE) 

U  ( CM/ SEC  1 

MACH  NO 

PHI  (6ed)  P  <DVNE/C«-*Jt 

RHO  (QM/CCl 

f  1060  K) 

PS1  ( GM/CM* "2SEC ) 

0 

3 . 57 1 6 1 78E  *00 

1 .80974246*00 

3  470004  76  *03 

0. OOOOOOOE *00 

1 . 84273046*05 

1  6  19828BE  05 

3.49890201*03 

0.00000006*00 

1.27069766*11 

2  53009916*01 

3  00234946  *00 

1.7432646?  OB 

2 

3.57161786*00 

1  39  15990F*  1  t 

7  T 1 2R5 166  *00 

2  5726545F401 

3 . 9247076E  *05 

3  49 1B577E100 

1  250I963E*00 

1  89104376*05 

3 . 3B7B2  14E403 

1 . 75766606  *01 

3 

3  57 16 1 7BE  *00 

2-  16641526*00 

2  58472576  *01 

3. 920253  IE  *05 

3  5  142407F  *00 

1.55358266*00 

1  91  150216*05 

1  78543306  05  3.32922686*03 

2.25593116*01 

1  37677356*  1  1 

"4 

3.57 16l7BE*0O 

2  2147 1416  *00 

2  596021 16*01 

3  91471341*05 

3  5333 160f ♦OO 

1.862  14956*00 

1  93484106*05 

1  B292354E  05 

3  30355406*03 

3 . 57  161786*  OO 

2  -  2 758044 E*00 

2  *1079*1 26*01 

3 . 90504  3 1 E  *05 

2.30794086*00 

1.97257726*05 

1.09419756-05 

3  27097706*03 

3 . 35  *44896*01 

1 .343439  IE*  11 

3.55573356*00 

T.  32 370406 *00 

2 . 622989  16  *0 1 

7  7021266E*00 

3  O09B473E-OB 

1.95442986  05 

3  2446724E*03 

3  97276046*01 

6 

3.57'6178E«00 

3 . R968044F  *05 

3  57238956*00 

3 . 57 16 1 78E  *00 

2 . 39863316*00 

2 . 64244926  *01 

2 . 44575351^00 

2 . 6556639E  *01 

3  87943616*05 

3.44524016*00 

2.08844326*05 

2.0722242E  05 

3  20398736*03 

4  7419FB4F-0I 

1  3009274E* 1 1 

3, 5945624F 1 00 

2.  1 4  86  344E  *05 

2.  16121286  05 

3  17640S0f'05 

5732873096*01 

8 

3.57161786*00 

3  8665 16OE  *05 

3  608 1023E  *00 

3  96529941*00 

_9 _ 

3 . 57 16 17BE  *00 

1 . 29253276*1 1 

2.49240506*00 

2 . 66934556 ♦61 

4.54262516*00 

2  22 15234E  *05 

2  2B66264E  05 

3  14755976*03 

5.94544946*01 

— 

3  62042226  *00 

5  «9t6202E*0O 

2.5BT«29lEi6S 

274529474^05 

3  1034280?  *03 

67927805  IE  *01 

to 

3  57 16 1786*00 

1  2317 130E* »  1 

7  5609570E400 
2.69004 15E10I 

3  83062836*05 

3  8367499F  *00 

3. 67 161786*00 

1  .20933«7€~*11 

3  8228 19BL  *05 

6  5BB92 151 *00 

2  5  1 75326E  *05 

2  6858*246  05 

3 . 0568702E  *03 

8.07463276*01 

7. 71030436+01  ' 

3  672249BE*CO 

3.57161786*00  2 . 7 1 17626F *00 

1  2081745F*11  2  73078756*01 

3  87495656  *05 

1 . 7626 1 ?4f  *00 

7  90nS7OE<00 

2  74197656*05 

2  9942651605 

3  00857166*03 

3.57 161786*00 

2  80125026*00 

3 . 943749  IE  *05 

9  44957296*00 

3  05 10622E  *05 

3  44137836-05 

2 . 9403 !50E*03 

1  14436096*02 

1  J07S7B7E* 1 1 

2  7508  14  16*01 

3.884779BE»66 

3  57161786*00 

2  B803B8O? *00 

2  70349O4E  *01 

4.01962946*05 

4 . 02482286*00 

».09J4383E«01 

3  4C«6690E*O5 

3  9934431605 

2  6S6744IE*03 

3.57161 78E  *00 

2  9733860E  *00 

4.  109527BE  *05 

1  23057 146*01 

3.8H8269E*05 

4 . 7 1232536 -05 

2  73793066*03 

1.64343806*02 

1  JOBSP3BE-MI 

2.81335206*01 

4  204056BE  *0O 

<i."379796»E*6s 

B  63BO366E-0B 

3  5716178F*O0 

1  20923556* 1 1 

3.05669026*00 
2. B38375U*01 

4  19496116*05 

4  39314006*00 

1  3595818C-01 

2  59816396*03 

17 

3 . 57 1617BE*00 

3.  M22820E *00 

2 ; 8595594 E* 01 

4  27867116*05 

1  4854366F  *0 1 

4.85469016*05 

6  83493926  -05 

3  44101046*03 

2.40149236*02 

1  S1J5035EMI 

4  59624186*00 

If 

3 . 57 16 1 78 E  *00 

1 . 2130494E* 1 1 

3  230B893E  *00 
2  86888946*01 

4 . 35C9166E  *05 
4  B335274E  *00 

1  608156.  1 

5  58987986*05 

8  57879126  06- 

2  24804046*03 

2.97175  *02 

19 

3.5716178E*O0 

3. 32408206*00 

4  43364926*05 

1  7 3QQ694 E *01 

6 . 599  1 2586*05 

1 . 13990296-04 

2  00239836*03 

_ 3. 78 107296*02 _ 

1 . 2  1 22B2BE  +  11 

2  87  49670E  *0 1 

5  160739 16*00 

1 . 5954 I05E  04 

"20 

3.57 101786*66 

1 . 2 1 196B5E*  1  1 

3  42319646*00 
2  8732957E  *01 

4  507664  .T  *05 

5.57786956*00 

1 .85457386*01 

8  094 1 397  E  *05 

21 

3.57  16  1 78E*00 

3  43 1 55896  *00 

4  5  1 288 1BE *05 

1 ,8633039F*01 

8  19830856*05 

1.64 101936  04 

1 .72937946*03 

4  94 10B08E*02  _ 

1  2  M6389E*  1 1 

2  87  J9062F  *0 ! 

5  6 1 170286*00 

ELAPSEO  TIME  ■ 

13  25  SEC 

DRAG  -3  7396*07  DYNES. 

.  , ...  . 

- - - - -  — 

e.  Method-of-Characteristics  (MOC)  printout,  ray  normal  to  cylinder  at  beginning  of  flare 

Figure  F-4.  Continued. 
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18 

3  571617BE*00 

1  2 I38494E4 1 t 

3  23O0893E*OO 

3  .  8686894E *0 1 

4 

4 

356916GC*0fi 

8  335274E ♦ OO 

1.808156  1 

5 . 589H79BE *05  B.5787912E  05 

2  24884O4E»03 

2  97171  02 

_ 19 _ 

3 . 57 16 1 78E ♦ OU 
1.2?22828E*  ii  ' 

3 . 32 408 JOE  *00 

2 . 8749670E  +Q 1 

5 

4336492E«05 
160719  IE *00 

1  7300694E40I 

6 . 599 1258 E *05  1  1399O29E-04 

2 . 0023963E  *03 

_3. 78 107 26E  *02 _ 

20 

3  5716178E«00 
1.21 16685E  *11 

3 . 4  233964E  *00 

2 . 8732957C  *0 1 

4 

5 

507664 7E  *05 
5778695F  <00 

1  B54573BF*01 

fl  054 1397E*05  1.5054IO5E  04 

1  7451 197E.03 

5 . 0097734C  *02 

21 

3.5716178E»00 

343 15589E ♦OO 

4 

51288 1 8  E ♦ 05 

1  86330391401 

8  1 96 3 08 5 E  *05  I.64  101D3E  04 

1  7263794E*03 

4.94 I0808E  -02 

1  2HG389E*  11 

2  8729B62E  *OI 

5 

6li7028FtOO 

ElAPSEO  TIME  ■ 

13.28  SEC 

DRAG  *J  739£*07  OfNtS 

4  004  DEG  FLARE  INlilAL  ANGLE  OF  IM0; Ollf 0  SHOCK  IS  72 .  *02  DEGREES 
MAFCM  ERROR  •  O  510981  09 


f.  Printout  from  Subroutine  ‘FLARE’ 
Figure  F-4.  Continued. 


RAV  R&N  "  SIGMA  PI  DELTA  S8AR  7  BODY  ~  S  BODY  PHI  BOOY 

61  EQUlLl  2  260BB84C  *0 1  0 - OOOOOOOE >00  O  OOOOOOOr  »0O  0 .  OOO  JOOOE  ♦  00  3  57I6I7SE400  I  8897424E»00  O  OOOOOOOE^OO 

SURFACE  DISTANCE-  6  M93  FiET  ENTROPY'  0  452619F*04  CAL/GM/K 


. 

■Rrt 

M) 

H  (ERG/GM) 

fl  ICMI 

XMW  (GN/MOLE) 

U  (CM.'SFCI 

MACH  NO 

PHI  iOEGj 

P  t OYNE/Cli.  *2  J 

RHO  1  GM/CC  1 

T  (DEG  Ki 

PS!  |dH/cS»»55fCj 

O 

3.5716178E*00 

1  6897424C  *00 

3. 383 12521*00 

4.004I729€*01> 

2  3B2I007E1O5 

2  0  12830  *E  05 

3.5789782E>03 

o.oooooooe*oo 

1  2 7 0696 3 E ♦ 1 1 

2  r3t2rc41E+Ol 

7  71f*2  I9BE  «00 

-  o 

3  57 1G 178C  *00 

1  389 /424E i 00 

3  4 7 0804 2 F  -05 

6  OOOOOOOF  *oo 

1  8427304E  *05 

1  619B2RBF  05 

3  46390201 *03 

O.OOOOOOOE4QO 

1 . 2706976C* 1 1  2  530899  I E*0 1  _3  0023494E-00 


IMRCOOf D  SHOCK  ANGLE  *  32  107  DECREES 


NRCORD  _*  _S8  .  _ 

g.  Flow  conditions  at  discontinuity  (first  block  is  behind  the  shock) 
Figure  F-4.  Continued. 
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RAY 

62 

RUN 

COUILI 

SIGMA  PI  DELtA  SOAR  7  BOO* 

2  1 D39355E*01  O  OOOOOOOE *00  0  OOOC000F*O0  O. OOOOOOOE *00  3  B5329I3C* 

R  BODY  PHI  BODY 

00  1  9094595E*00  40041729E+00 

SURFACE  01SIANCE*  0  1585 

FEET  ENTROPY- 

0.452619E-04  CAL/GM/K 

NPT 

Z  ( CM  I 

H  (ERG/GMI 

P  <CM> 

7  MW  (RM/MOLEt 

U  1  CM/SEC  1 

MACH  NO. 

PHI  |DEO) 

P  IDYNE/CM--2) 

oi  in  (  S5/CC  l 

T  IDEG  Kl  PSI  IGM/CM--2SEC) 

0 

3. 85329*3E*00 

1  9094595E  *00 

3 . 3882339E  *03 

4  0O41729E*O0 

2 . 34738 15E  *05 

1.987826  IE -05 

3.570494  IE *03 

0  OOOOOOOE *00 

1 . 2704 13  IE 1 11 

2.51318  7BE  *01 

2. 864701 7F*0O 

1  3910075fToi 

t  aaOaTjsF-oi 

lOOO 

3.B532913E*00 

B . 26S2335E  *  10 

2  004 1 299E  *00 

2 . 8753026E  *01 

3  4B10096F  >05 

4  0684 7 89E  *00 

B  79947  1  t  E  *09 

1  9371982E*03 

B  7B8330 IE*00 

lOOO 

3 . B5329 1 3E  ♦  OO 

2.004 I299E+00 

□ . 7050730E >05 

5  272B956E  01 

1  7966706E  *05 

5.0270998E  05 

1 . 23435 17C *03 

8  788330 IE +00 

8  1 950046E ♦ 10 

2 . 86  •  5346E  *01 

5  3858427E  *0O 

1 .  0773OME*66 

T.  69 14754?  05 

2 

3  85329 1 3E *00 

1  39 1 6387F ♦ 1 1 

7  •  *85~ 36F*00 

2 . 5753059E  *0 1 

3  934347GE >05 

3 . 5  I03895E  >00 

1  8260359E>O5 

3  344B021E«03 

3.7376783E*01 

3 

3  85329 1 3E  *00 

2. 1 7  35639E  *00 

3 . 93060 1 7E  >05 

1  354O528C*00 

1  840158  IE >05 

1.7277420C  05 

3 . 3  15732 1 E  *03 

3  226  1 344E  *01 

1  376B035E *11 

2.5876392E*01 

3  53419771 >00 

i  0S44 IB3E .OS 

4 

3  85329 1 3E *00 

1  3625609E  *  1 1 

2 . 22326471 • 00 

2  5992989E  *01 

3.926 1377E >05 

3 . 5553998E  >00 

1  6 1539 18E  *00 

1  7632853E  OB 

P) 

? 

* 

" 

3  69677 IIE*Ot 

5 

3  8532913E*00 

7.286 I5B9E*00 

3 . 9 1RG593E  >05 

1  9838249E  *00 

1  870 1686E  *05 

1  8153966F-05 

3  2543692E *03 

4 . 32 192 16E*0I 

1  3434568E * 1 1 

2  6 146039F  *0 1 

3  5R09657F»Ot) 

*  3  22654  13? *65“ 

a 

3  85329 1 3E  *00 

1  32732  1 3E.  ♦  1  1 

2  336O5B0E  *OC 

2  G272423E  *0 1 

3  91112554*05 

3 . 6002 1 24 E  *00 

2  32C8388C  *00 

1  902058  7  E ♦ 05 

1  R64U9I IE  05 

4  B435003E  *01 

7 

3 . 85329  1 3E  *CO 

2. 4  142646E*00 

J  89bB190€>05 

2  9073731E<00 

1  952476«E*05 

1  9544628E  05 

3.  102171 7E+03 

5.71278771*01 

1  1006  4  8 1 E ♦ 1 1 

2. 6476919E*01 

3 . 62R4459E  >130 

a 

3  85329  1 3E+00 

1  2B223*lE*11 

2.463789OE*0O 

2  66 1 2  205E*01 

3 . BB50791 E  >05 
3  64436501*00 

3 . 3617421 E  *00 

1  998046«C*0S 

2  02«83»9E0B 

3.  I530997E  *03 

8  30OI2B0C70I 

9 

3  85329 13E*00 

2  5 1 3045 1E*00 

3  8  722  305  L ♦ 05 

3 . 83B94  25E  *00 

2  04 7394 2E *05 

2  1 1057 72 E  *05 

3. 12 20904 €  >03 

6  9 1722 1 7E  *0 1 

1  262B26BE+1 1 

2  6751573E*01 

3.66031634*00 

j  139Bi3~tt  *05  " 

2  2B9?9i4t-OB 

10 

3  85329 13E *00 

1  23  182  13E*  1  1 

7  5859347E  *00 

2  6966564E«01 

3 . B500R47E *05 

3  682B684F  >00 

4  642 1741 E*00 

3 . 0T25O67E  *03  7 . B99B390E+01 

11 

3 . 85329 1 3E  *O0 

2 . 6G296 17E*00 

3  8485299E  *05 

5  6 1R6 152E  *00 

2  264 8503 E  *05 

2  4  49258 1 E -  05 

3  02306 96E*03 

9 . 0468463E  *0 1 

1  2094839E *11 

2  7173454E«01 

3  7234599F  MX) 

3  704343  IE  OB 

12 

3  85329 1 3E *O0 

1  20B3S33EM1 

2 . 748 1095E  *0Q 
2.73812 15E*01 

3  902362  IE *05 
3,81886861*00 

6  B035SJ7E.0O 

2  4  390634E  *05 

2  9710B8BE+03 

1  04854 10E  *02 

13 

3  85329 1 3E  *00 

2. 945208  IE  *00 

3.97I2444E105 

8  2885907F-0O 

2  6964 706 E* 05 

3  0934857E-05 

2  898686  4  E *03 

1.24 1 8955E  *02 

1 . 2073802E ♦ 1 1 

2. 7641 147E*0i 

394406 / 1 E *0O 

9  1652979E  *00 

3  5B56093E - 06 

2  8 12406  3EM>3 

14 

3  .  B5329  1 3E  ♦•00 

1 .2076472EM  1 

2  93O0O99E*OO 
2 . 79Q4898E+01 

4  0460604 F *  05 

4  086 2 4 37 E  *00 

3. 0037454 E *05 

1 ,4696I03E»03 

15 

3  85329 1 3E *00 

3  0319395E-00 

4  1 34 1022 E  *05 

1  1  178439E*0 1 

3  3599496E*05 

4  2342167E05 

2  59 1 370BE  *03 

1 .74196 135*02 

1  7GR4B52E  *11 

2  B  19 1 370C  *01 

4 . 26576031 *00 

16 

3  853291 3E *00 

1  209347 1 E*  1  1 

3  12 1989 IE  *00 
?  B475794E401 

4  2174902E*05 
4  455G564C  *00 

1  2502935E«01 

3  77465B5E  >05 

5  O639650C-05 

2  549 I609E  *03 

2.07341 17E*02 

17 

3  853291 3E  *00 

3.214 1965E*00 

4 . 2990374E  *05 

1 . 3787697E  *01 

4  2880689E  *05 

6.  1714340C-05 

2  390207 8 E *03 

2  50 1 8093F *02 

1  2I27«7U»M 
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Ak 

Aki 

AU 

Bk 


Bki 

bi 


Ch 


CP 

CPj 

Cf 

Cki 

Dj 

Dik 

Eif 

Ej 

AFi 

F1-F5 


NOMENCLATURE 
Constant  in  radiation  rate  of  ktb  process 
Constant  in  expression  for  rate  of  ith  reaction,  Eq.  (21) 

Matrix  specifying  reactions  affected  by  j*  vibrator 

“Frozen”  sound  speed;  aj2  =  Cp  <=/<?  E-yj  T/(Cp  -  E  yj) 

Constant  in  radiation  rate  of  k*  process 

Constant  in  expression  for  rate  of  ith  reaction,  Eq.  (21) 

Constant  in  harmonic  oscillator  formula  for  chemical  potential 

Stanton  number 

Specific  heat  of  mixture 

Specific  heat  at  constant  pressure  of  jth  species 

Coefficient  of  friction 

Constant  in  expression  for  rate  of  ith  reaction,  Eq.  (21) 
Constant  specifying  reaction  type,  Eq.  (22) 

Constant  in  expression  for  rate  of  i*11  reaction,  Eq.  (21) 

Binary  diffusion  coefficient 

Energy  of  1th  electronic  level  of  jth  species,  cal  per  molecule 
Molar  free  energy  of  jth  species,  ergs/mole 

Change  of  free  energy  in  ilh  reaction  in  Eq.  (19) 

Protile  integrals  used  in  boundary-layer  calculation 
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G  Boundary-layer  enthalpy  ratio 

<S  Rate  of  change  of  static  enthalpy  due  to  chemical  and  vibrational  relaxation, 

Eq>  (A-24) 

gj t  Degeneracy  of  1th  electronic  level  of  j,h  species 

H  Total  enthalpy 

hj  Static  enthalpy  of  jth  species,  ergs/mole 

h?  Heat  of  formation  of  jlh  species,  ergs/mole 

Ij  Rotational  contribution  to  species  enthalpy,  ergs/ mole 

J  Jacobian  of  transformation  from  z,  r  to  x,  y  coordinates 

K;  Equilibrium  constant  ilh  reaction 

k'bi  Backward  reaction  rate  constant  ith  reaction 

k'fj  Forward  reaction  rate  constant  iIh  reaction 

<M  Molecular  weight,  gm/mole 

me  Mass  of  electron,  gm 

M  Ray  index,  also  any  species 

Nc  Number  of  elements 

Nr  Number  of  reactions 

Nrad  Number  of  radiative  processes 

N§  Number  of  species 

Nvib  Number  of  nonequilibrium  vibrators 

nj  Number  of  atoms  in  jth  species 
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Ne 

Nj 

o 

p 

p 

Qii 

R 


Rad 


s 

s 

T 


‘Vj 


U 


u 


Vj 


W; 


Electron  number  density,  cm-3 
Number  density  of  jth  species,  cm-3 
Number  of  vibrational  levels  in  cut-off  oscillator 
Pressure,  dynes/cm2 
Dimensionless  pressure, 

Source  term  in  species  conservation,  Eq.  (22) 

Radius  of  curvature 

Universal  gas  constant  =  8.313405  x  107  ergs/mole  °K 

Coordinate  normal  to  axis 

Total  energy  radiated,  ergs  gm-1  sec-1 

Streamline  coordinate 

Boundary-layer  history  parameter 

Translational  temperature 

Vibrational  temperature  of  jth  oscillator 

Resultant  velocity  =  U/Ua 

Free-stream  velocity  at  centerline 

Velocity  component  normal  to  ray  =  u/U 

Modifies  reaction  to  account  for  nonequilibrium  of  jlh  vibrator 

Velocity  component  along  ray  =  v/U^ 

Source  function  for  nonequilibrium  vibration  -  defined,  Eq.  (16) 
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Wi 

Xi 

X 

y 

z 

Zi 


<*jk 


0 

r 


r, 


•yj 

5 

€ 


T 

Ta 

7b 

7c 

7d 


Constant  specifying  reaction  type 
Mass  fraction  of  i'h  species 
Coordinate  along  shock 
Coordinate  normal  to  shock 
Coordinate  along  axis  from  shock  vertex 
Constant  specifying  reaction  type 
Number  of  atoms  of  type  k  in  species  j 
Mach  angle,  also  cos2  <t> b 
Averaged  specific  heat  ratio  (frozen) 

Concentration  of  j111  species,  moles/gm 
Concentration  of  j,h  species,  moles/free-stream  mole 
Boundary-layer  velocity  thickness,  cm 
Vibrational  energy,  erg/gm 
Viscosity,  gm  cm-1  sec-1 

Chemical  potential  of  jth  species,  ergs/mole,  Eq.  (4) 
Stoichiometric  coefficient  of  species  j  in  reaction  i 

Vibrational  relaxation  time,  cm,  Eq.  (27) 


j.'lju  rl  1 L '! n  muitdifinpVi'Oif  id!  mows  r>:  i.r.iv 

Constants  in  expression  |pr  ^Eq.  ,(27),^,  0 

(oi)  .p3  J/'Kiftoi:  -  noit&id'v’  ntjmdiliup'."i.-i;  v fioitar;i't  'jyiHr?. 
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V 

0ii 
0 

ir(s) 

7T 

X 

X 

X 
A 

a 

<7 

U) 

a 

e 

Xi 

Subscripts 

oo 

o 


Characteristic  vibrational  temperature,  °K 

Boundary-layer  coordinate  normal  to  body;  also  coordinate  normal  to  a 
streamline 

Viscosity  cross-section  (Az) 

Flow  angle  with  respect  to  z-axis;  azimuthal  angle 
Mass  flux 

Boundary-layer  pressure  gradient  parameter 

(1/s)  J*  T(s)ds 

Local  shock  curvature 

Constants  used  in  determining  new  shock  shape 

Uoo2  constant  in  dimensionless  state  equation 

Shock  angle  with  respect  to  Z-axis 

Stefan-Boltzmann  Constant  =  5.6687  x  10“ 5  erg  cm-2  (°K)"4  sec 

0  two-dimensional  flow 
1  axisymmetric  flow 

Solid  angle 

Density,  gm  cm-3 

Degree  of  nonequilibrium  of  ith  reaction,  Eq.  (24) 

Evaluated  in  free  stream  at  centerline 
Evaluated  at  stagnation  point 
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b  Evaluated  at  body  surface 

sh  Evaluated  at  shock 

e  Evaluated  at  boundary-layer  edge 

N  Normal  to  shock 

t  Tangent  to  shock 

am  Local  free-stream  value  at  shock  location 

Superscripts 


_  Indicates  dimensional  quantity 

o  Evaluated  at  0DK  and  1  atm  pressure 

*  Pertaining  to  products  of  a  reaction 

Note:  The  nondimensional  quantities  used  throughout  this  report  are  referenced  to  free- 
stream  values  according  to  the  following  system: 


Types  of  Quantity 

Length 

Velocity 

Density 

Pressure 

Energy  of  mixture 
Energy  of  particular  species 
Molecular  weight 


Reference 


R&h 


IF 


G  oo 

G„Ui 


T 

ca  a 


T0 


NS 

Concentration  (yj)  E  =  cM~ 

j=i 
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